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瞬 变 电磁 法 以 其 独特 的 优点 广泛 应 用 于 资源 勘查 、 工 程 勘察 、 地 质 灾害 调 
查 等 领域 ， 并 越 来 越 受 到 广大 地 质 科 技工 作者 的 青睐 。 随 着 把 地 震 学 中 已 有 的 
和 正在 发 展 的 方法 拉 术 应 用 到 肯 变 电磁 法 中 的 研究 不 断 进步 ， 瞬 变 电 磁 测 深 偏 
移 成 像 方法 技术 的 研究 也 成 为 当今 勤 控 电磁 领域 研究 的 热点 。 人 们 越 来 越 关注 
像 拟 地 震 解释 法 、 电 磁 俩 移 成 像 以 及 偶 移 与 反 演 相 结合 的 研究 方法 。 这 些 方 法 
从 波 场 的 角度 ， 拓 展 和 丰富 了 有 瞬 变 电 磁 场 的 内 浮 ， 使 得 从 实测 资料 中 提取 到 和 常 
规 瞬 变 电 磁 法 提取 不 到 的 、 对 地 下 目标 体 成 像 更 有 利 的 信息 。 

面 对 资 源 勘查 深部 矿 、 讶 矿 ) 和 工程 勘察 (精细 勘察 ) 方面 对 地 球 物理 
工作 要 求 的 提高 ， 人 研究 开发 瞬 变 电磁 法 拟 地 震 的 多 次 覆盖 技术 是 十 分 必要 而 且 
有 重要 实际 意义 的 。 男 外 ， 多 分 量 、 阵 列 式 的 观测 技术 可 以 获得 丰富 的 地 电信 
恩 ， 其 方法 技术 研究 也 已 经 引起 人 们 的 重视 。 研 制 抗 干扰 能 力 强 的 多 分 量 、 多 
通道 适合 开展 阵列 式 观 测 的 瞬 变 电磁 探测 系统 ， 大 功率 、 多 功能 、 智 能 化 及 适 
时 数据 处 理 功 能 等 是 现代 仪器 的 发 展 方向 ， 另 外 ， 研 制 高 灵 人 敏 度 、 高 性 能 、 三 
分 量 观 测 探头 也 是 高 分 辨 退 变 电 磁 仪 的 重要 环节 。 

在 前 人 研究 的 基础 上 ， 笔 者 先后 和 二 十 多 位 研究 生 一 起 研究 瞬 变 电磁 拟 地 
尾 资 料 解 释 方 法 技术 。 提 出 从 瞬 变 电磁 场 到 波 场 的 优化 算法 ， 将 正则 化 算法 用 
于 逆 变 换 过 程 的 计算 中 ， 成 功 地 提取 出 虚拟 地 震 子 波 。 通 过 实验 发 现 ， 有 来 用 多 
孔径 辐射 场 源 能 得 到 比 单 孔 径 辐 射 场 源 更 强 的 磁场 分 布 这 一 现象 的 存在 ， 从 一 
次 场 观测 结果 看 出 多 孔径 场 明显 地 改善 了 磁场 的 方向 性 ， 证 明了 多 孔径 瞬 变 电 
磁场 存在 相关 人 加 现象 ;完成 了 虚拟 子 波 宽 度 压缩 技术 研究 ， 用 基 尔 上 霍 夫 积分 
法 进行 电磁 波 场 偏 移 成 像 处 理 ， 实 现 了 瞬 变 电磁 场 的 三 维 曲面 延 拓 计算 ; 从 瞬 
变 电 磁场 本 号 的 特性 出 发 ， 基 于 等 效 导电 平面 的 理论 ， 提 出 瞬 变 电磁 虚拟 波 场 
连续 速度 分 析 方 法 ， 初 步 取 得 了 瞬 变 电磁 拟 地 震 资 料 解释 的 阶段 性 成 果 。 

我 们 出 版 本 研究 专著 ， 以 便 对 该 方法 技术 感 兴趣 的 读者 阅读 和 参考 。 全 书 
共 分 7 章 ， 第 1 草 介 绍 瞬 变 电磁 法 的 发 展 概述 和 研究 展望 ;第 2 章 为 波 场 转换 
基本 原理 ; 第 3 章 人 研究 瞬 变 电磁 虚拟 波 的 波 场 特性 ; 第 4 章 介 绍 合成 孔径 技术 
及 虚拟 子 波 锅 度 压 缩 技 术 ; 第 5 草 介绍 瞬 变 电磁 虚拟 波 场 三 维 曲 面 延 拓 成 像 技 
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术 ; 第 6 章 介 绍 虚 拟 波 场 的 速度 分 析 ; 第 7 章 给 出 理论 模型 数值 模拟 结 末 以 及 
我 们 近年 来 的 工作 实例 。 

全 书 主要 由 李 儿 负责 撰写 ， 醇 国 强 负责 第 1 音 、 第 4 章 的 撰写 和 全 书 的 编 
排 和 校对 工作 。 往 届 毕 业 的 研究 生存 志 鹏 、 郭 文 波 、 马 宇 、 全 红 娟 、 武 军 杰 、 
范 涛 、 刘 银 爱 、 吴 琼 、 朱 宏伟 、 张 旭 等 ， 以 及 在 读 研究 生 孙 怀 风 、 钱 建 兵 、 
张 焉 芝 等 者 对 方法 的 研究 以 及 本 书 的 出 版 做 出 了 贡献 。 博 士 生 周 楠 楠 对 眠 变 电 
磁 物 理 量 信息 进行 了 整理 , 博士 生 陈 卫 营 对 公式 进行 了 核对 。 硕 士 生 李 海 对 书 
中 参考 文献 进行 了 校对 。 

本 书 是 在 作者 近年 来 所 承担 的 国家 上 自然 科学 基金 项 目 和 企业 委托 项 目的 基 
础 上 完成 的 。 但 由 于 瞬 变 电磁 法 在 我 国 开 展 研 究 和 应 用 的 时 间 相 对 较 晚 ， 特 别 
是 瞬 变 电磁 拟 震 解 释 理论 与 应 用 研究 还 处 于 初步 阶段 ， 还 需要 继续 这 人 研究 。 
加 上 作者 水 平 有 限 ， 不 妥 之 处 在 所 难免 ， 敬 请 专家 和 读者 批评 指正 。 


李 fk 
2012 年 11 月 8 日 于 长 安 大 学 
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1.1 EBRA AREE 


瞬 变 电磁 法 最 早 在 20 世纪 30 年 代 由 原 苏 联 科学 家 提出， 当时 采用 的 是 远 
区 工作 模式 。 但 是 ， 利 用 电流 脉冲 激发 供电 侦 极 形成 时 域 电磁 场 是 由 美国 科学 
家 L. W. Blau 于 1933 年 最 先 提 出 的 ， 当 时 利用 不 同 电 导 率 地 层 界面 电磁 波 的 
反射 与 地 震 反 射 波 信和 号 的 相似 性 ， 进 行 了 大 量 的 实验 和 比较 ,但 由 于 脉冲 激发 
的 瞬 变 啊 应 频率 较 低 ， 难 以 得 到 识别 各 个 反射 波 所 需要 的 分 辩 率 ， 瞬 变 电 磁 法 
一 直 受 到 冷落 。 到 了 20 世纪 五 、 六 十 年 代 ， 苏 联 科学 家 成 功 地 完成 了 瞬 变 电 
磁 法 的 一 维 正 、 反 演 ， 建 立 了 瞬 变 电磁 法 的 解释 理论 和 野外 工作 方法 之 后 ， 瞬 
变 电 磁 法 才 开 始 进 入 实用 阶段 。20 世纪 60 年 代 以 后 ， 当 意识 到 时 间 域 电磁 测 
深 法 可 以 利用 远 远 小 于 期 望 探 测 诬 度 的 收发 距 时 ， 该 方法 有 了 一 个 快速 发 展 。 
在 办 联 ， 组 建 了 数 十 个 瞬 变 电磁 法 探测 队伍 ， 在 油气 勘探 领域 开展 工作 ， 并 取 
得 了 齐 有 成 效 的 勘查 效果 。 随 之 ,“ 短 偏 移 “晚期 ”、“ 近 区 ”等 技术 研究 迅 
速 发 展 起 来 。20 世纪 七 、 八 十 年 代 ， 在 美国 等 西方 国家 短 偶 移 法 一 直 处 于 研 
究 和 实验 阶段 ， 未 被 广泛 运用 ， 而 长 偏 移 法 已 得 到 了 应 用 ， 特 别 是 在 地 热 调 查 
和 地 元 绪 构 的 调查 中 。 随 后 一 些 专家 对 瞬 变 电磁 法 的 一 维 正 、 反 泪 方 法 技术 进 
行 了 大 量 研究 。20 世纪 80 年 代 后 随 着 计算 机 技术 的 发 展 ， 欧 美学 者 在 二 、 三 
维 正 演 模 拟 技 术 方面 ， 发 表 了 大 量 用 有 限 差 分 法 、 有 限 元 法 、 积 分 方程 法 等 进 
行 电磁 法 数值 模拟 研究 的 论文 。 理 论 和 认识 的 提高 推动 了 瞬 变 电磁 法 探测 应 用 
的 发 展 。 电 磁 法 理论 方面 的 主要 代表 著作 有 考 夫 曼 和 凯 勒 (1987) 的 专著 《 频 
域 和 时 域 电磁 测 深 》， 米 桂 殉 N， 纳 比 吉 安 (1992) 主编 的 《勘查 地 球 物理 。 
电磁 法 》 第 一 卷 。 与 此 同时 ， 苏 联 学 者 别 尔 季 切 夫 斯 基 和 日 丹 诺 夫 (1985) E 
出 电磁 波 拟 地 震波 的 偏 移 方法 ， 他 们 吸取 了 “ 仿 移 成 像 ” 的 广义 概念 ， 在 电磁 
法 中 确定 了 正则 偏 称 和 解析 法 偏 移 两 种 方法 。 随 后 ， 一 些 成 像 技 术 研 究 的 成 果 
相继 出 现 ，Macnae (1987). Eaton 和 Hohmamn (1987, 1989) 等 完善 了 电 
导 率 成 像 方法 ，Ghristensen (2002) 对 瞬 变 电磁 数据 一 维 成 像 进行 了 研究 。20 
世纪 90 年 代 初 ，Lee 和 Xie (1993) 推导 了 二 维 成 像 解释 方法 。de Hoop 
(1996) 认为 瞬 变 电磁 法 与 反射 地 震 勤 控 具 有 相似 性 原理 。Gershenson (1993) 
提出 用 波 的 扩散 传播 特性 解释 时 间 域 电磁 测 深 资料 。Kunetz (1972) 研究 了 
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大 地 电磁 与 地 震 弹 性 波 场 之 间 的 内 在 联系 ，Lee 和 Memechan (1987) 根据 
电磁 波 与 地 震波 的 相似 性 利用 差分 法 的 二 维 偏 移 实现 了 电磁 数据 成 像 。20 
世纪 80 年 代 末 ， 别 尔 季 切 夫 斯 基 和 日 丹 诺 夫 (1985) 又 从 激发 极 化 现象 理 
论 出 发 ， 研 究 了 时 间 域 瞬 变 电磁 法 的 激 电 效应 特征 及 影响 ,解释 了 瞬 变 电磁 
法 晚期 电磁 响应 的 变 号 现象 ， 并 对 三 维 极 化 体 的 瞬 变 电磁 响应 特征 进行 了 数 
值 计 算 。 

在 我 国 ， 对 瞬 变 电磁 法 的 研究 始 于 20 世纪 70 年 代 ， 较 早 开展 工作 这 项 工 
EHAZE (1986, 1992). #HE2S (1990) 、 曾 孝 敌 等 ， 他 们 将 脉冲 式 航 电 
仪 用 于 地 质 填 图 和 找 矿 中 ; 蒋 邦 远 (1998) 等 将 瞬 变 电磁 法 用 于 普查 勘探 良 导 
金属 矿 实践 中 ， 并 研制 出 了 电磁 脉冲 瞬 变 系统 ; Fes (1986, 1992) 将 瞬 变 
电磁 法 用 于 金属 矿 勘 探 上 ， 取 得 了 显著 的 效果 ， 并 与 长 沙 智通 新 技术 研究 所 合 
作 人 研制 并 生产 了 智能 化 瞬 变 电磁 仪 ; MER (1990) 用 G- SMREKE 
现 了 电 性 源 瞬 变 电 磁 测 深 的 正 演 计算 ; SARE (1993) 用 线性 数字 滤波 技术 
实现 了 大 回 线 源 瞬 变 电 磁 测 深 的 正 演 计算 ， 并 将 大 回 线 源 瞬 变 电磁 测 深 法 广泛 
用 于 地 热 和 地 下 水 调查 、 工 程 调查 和 地 质 灾 害 调查 ， 又 将 瞬 变 电磁 法 用 于 大 地 
电磁 测 深 曲线 的 静 校 正 ， 取 得 了 良好 的 效果 。 在 一 维 、 二 维 数值 模拟 方面 ， 王 
华军 和 罗 延 钟 (2003) 从 事 人 研究 并 发 表 了 相关 论文 ， 王 华军 (2003). AN 
等 (1998)、 陈 明生 (1999)、 惧 长春 和 刘斌 (1994) 等 对 二 维 、 三 维 瞬 变 电磁 
正 、 反 演进 行 了 研究 。 王 家 映 等 (1985)、 王 家 映 和 方 胜 (1986)、 王 家 映 
(1990) 通过 分 析 平 面 电磁 波 与 地 震波 的 相似 性 ， 对 大 地 电磁 的 拟 地 震 解 释 进 
行 了 系统 的 研究 ， 于 鹏 和 王家 林 (2001), FHE (2003) 在 对 电磁 场 偏 移 研 
究 的 基础 上 ， 提 出 了 改进 的 有 限 差 分 法 大 地 电磁 场 偏 移 成 像 技 术 ; 李 杀 等 
(2005b) 对 瞬 变 电磁 测 深 开 展 了 微分 电导 成 像 的 研究 工作 ， 并 提出 从 瞬 变 电 
磁 到 波 场 的 优化 算法 ， 通 过 波 场 转换 ， 对 瞬 变 电磁 虚拟 波 场 的 三 维 曲 面 延 拓 成 
像 进行 了 研究 ; 陈 本 池 等 (1999) 对 波 场 变换 进行 了 研究 ， 用 模拟 数据 成 功 地 
获得 了 单一 界面 的 二 维 地 质 结 构 的 图 像 ， 王 华军 和 罗 延 钟 (2003) 对 瞬 变 电磁 
测 深 波 场 变换 中 的 波形 展 宽 原 因 探 讨 ; SEAE (1998) 研究 了 二 维 瞬 变 电 磁 道 
时 偏 移 ; BE (2005) 提出 并 开展 了 瞬 变 电磁 对 地 成 像 理 论 和 成 像 数值 计算 
方法 的 研究 ; 由 述 (2002) 用 时 间 域 有 限 差分 法 对 地 下 瞬 变 电磁 啊 应 进行 了 数 
值 计算 ， 并 在 瞬 变 电磁 测 深 的 联合 时 - 频 解 释 方 面 进行 了 研究 ; 郭 文 波 (2001) 
通过 分 析 均 匀 介 质 和 分 层 媒质 中 电磁 波 与 地 震波 的 传播 特征 ， 认 为 回 线 源 中 心 
激发 的 瞬 变 电磁 场 在 地 下 扩散 可 近似 看 成 球面 波形 式 传播 ， 在 晚期 球面 波 接近 
于 平面 波 ， 与 地 震波 有 相似 的 传播 特征 ， 在 时 - 频 转 换 的 基础 上 ， 通 过 求 取 发 
射 函数 序列 实现 了 瞬 变 电磁 资料 拟 地 震 解释 法 。 

国内 主要 的 著作 有 《电磁 测 深 法 原理 》( 朴 化 荣 ，1990)、《 脉 冲 瞬 变 电 磁 
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法 及 应 用 》 (和 牛 之 莲 ，1986) 、《 时 间 域 电磁 法 原理 》( 牛 之 和 莲 ，1992) 、《 瞬 变 电 
磁 测 深 原 理 》( 方 文 党 等 ，1993)、《 实 用 近 区 磁 源 瞬 变 电磁 法 勘探 》( 蒋 邦 远 ， 
1998)、《 瞬 变 电 磁 测 深 的 理论 与 应 用 》( 李 狄 ，2002) 等 。 

在 仪器 的 研制 方面 ， 国 内 多 家 科研 院 所 从 20 世纪 80 年 代 末 开始 相继 研 
制 成 功 了 各 具 特 色 的 瞬 变 电磁 仪 。 例 如 ， 中 国 地 质 科 学 院 物化 探 研 究 所 研制 
开发 的 WDC- 1. WDC - 2 与 IGGETEM - 20 瞬 变 电磁 仪 ; 西安 强 场 源 物化 
探 研 究 所 刘 显 跃 等 研制 的 LC- 1、EMRS - 3 等 小 线圈 、 大 电流 瞬 变 电磁 系 
统 ; 中 南 工 业 大 学 牛 之 焉 与 长 沙 智通 新 技术 研究 所 联合 研制 生产 的 SD- 1、 
SD-2 瞬 变 电磁 仪 ; 长 沙 白云 仪器 有 限 公 司 生产 的 MSD - 1 轻便 型 瞬 变 电磁 
仪器 ; 北京 矿产 地 质 研究 所 王 庆 乙 (1999) 研制 的 TEMS -3S 瞬 变 电 磁 仪 ; 
吉林 大 学 林 君 等 (2004〉 研 制 的 ATEM - 2 瞬 变 电磁 仪 等 。 很 多 勘查 单位 购 
买 了 国际 上 先进 的 装备 ， 如 美国 ZONG 研制 生产 的 GDP - 12、GDP - 16. 
GDP- 32 多 功能 电 法 系统 ， 加 拿 大 凤凰 公司 生产 的 V5、V6、V8 系统 ， 加 
拿 大 Geonics 公司 生产 的 PROTEM 瞬 变 电磁 系统 接收 机 (EM37、EM47、 
EM57、EM67 等 )， 澳 大 利 亚 生 产 的 SIROTEM -I + HD, Terra TEM 仪器 
等 ， 从 而 推动 了 瞬 变 电磁 法 的 工作 ， 并 获得 了 大 量 有 价值 的 研究 成 果 及 成 功 
的 应 用 实例 。 

1992 年 以 来 ， 随 着 仪器 的 智能 化 与 数字 化 ， 瞬 变 电 磁 法 开始 迅速 步 和 人工 
程 、 环 境 、 灾 害 地 质 调查 中 ， 如 探测 地 下 采 空 区 、 陷 落 柱 等 煤田 灾害 、 划 分 地 
下 断层 、 寻 找 地 下 水 、 金 属 矿产 勘探 、 石 油 和 煤炭 等 非 金属 矿产 调查 、 工 程 场 
地 地 质 勘 察 、 隧 道 超 前 地 质 预报 等 领域 。 目 前 ， 瞬 变 电 磁 法 已 经 几乎 涉足 了 其 
探 地 球 物 理 的 所 有 领域 ， 取 得 了 良好 的 效果 。 


1.2 了 瞬 变 电磁 法 的 发 展 方向 


瞬 变 电磁 法 作为 一 种 重要 的 电磁 勘探 方法 ， 以 其 独特 的 优点 (经 济 、 无 
损 、 快 速 、 信 息 丰 富 、 分 辩 力 强 等 ) 广泛 用 于 资源 勘探 和 工程 勘察 中 。 但 针对 
一 些 具体 的 精细 探测 问题 ， 如 高 速 公路 和 铁路 建设 中 的 隧道 超前 地 质 预报 精细 
探测 、 大 型 重要 证 墓 的 墓室 结构 精细 探测 、 大 型 水 坝 隐患 精细 探测 等 ， 背 用 的 
瞬 变 电磁 法 的 分 辨 率 受 到 限制 ， 探 测 效果 受到 影响 ， 这 是 因为 前 规 方法 要 想 提 
高 信 品 比 ， 就 要 加 大 发 射 回 线 的 面积 ， 加 大 发 射 磁 矩 ， 但 随 着 发 射 回 线 面积 的 
加 大 ， 体 积 效应 的 影响 随 之 加 强 ， 也 就 影响 了 分 辩 率 的 提高 。 从 目前 看 ， 在 理 
论 研究 、 仪 器 研制 方面 处 于 初级 阶段 。 理 论 研 究 虽 已 解决 了 一 维 正 、 反 演 问 
题 ， 但 在 二 、 三 维 的 研究 还 未 达到 应 用 程度 ， 而 仪器 大 多 是 在 国外 仪器 基础 上 
的 开发 和 改进 ， 应 用 推广 并 不 普遍 。 


瞬 变 电磁 法 拟 地 震 偏 移 成 像 研究 


随 着 探测 分 辩 率 和 精度 要 求 的 提高 ， 迫 切 需 要 从 两 个 方面 对 极 有 发 展 前 途 
的 瞬 变 电磁 法 进行 深入 的 研究 ， 一 方面 是 在 探测 的 方法 技术 方面 研究 多 孔径 阵 
列 式 观测 方式 ， 目 的 是 在 保证 一 定 探测 深度 及 不 增加 体积 效应 的 前 提 下 ， 提 高 
对 目标 体 的 分 辩 率 ; 另 一 方面 是 在 原 有 解释 方法 理论 基础 上 对 正 、 反 演 问 题 做 
更 深入 的 研究 ， 探 索 新 的 反 演 理论 并 对 其 进行 系统 化 ， 建 立 系统 的 瞬 变 电磁 解 
释 正 、 反 演 理 论 ， 同 时 使 瞬 变 电磁 法 的 解释 向 三 维 方向 迈进 ， 使 这 一 方法 能 更 
好 地 解决 一 些 高 难度 的 精细 探测 问题 。 

B 20 世纪 70 年 代 以 来 ， 瞬 变 电 磁 法 在 我 国 研究 、 应 用 已 有 30 余年 的 历 
史 ， 在 理论 〈 包 括 解释 理论 ) 人 研究、 方法 技术 方面 取得 了 长 足 进步 。 虽 然 先 后 
发 表 出 版 了 大 量 有 关 瞬 变 电 磁 理论 与 应 用 方面 的 论文 与 专著 ， 解 释 理 论 、 方 法 
技术 、 勘 察 效果 等 方面 均 有 报道 ， 但 从 目前 看 ， 在 理论 研究 、 仪 器 研制 方面 仍 
处 于 初级 阶段 。 

有 专家 认为 ， 正 是 由 于 电磁 法 理论 本 身 的 复杂 性 ， 对 用 户 在 仪器 设备 和 基 
础 理论 方面 的 知识 要 求 更 高 ， 因 此 方法 的 进展 相对 缓慢 。 国 内 仪器 生产 的 规模 
性 和 持续 性 较 差 ， 加 之 应 用 推广 并 不 普遍 ， 受 市 场 经 济 的 影响 ， 国 内 越 来 越 多 
的 生产 单位 和 科研 院 所 都 购置 了 国外 生产 的 多 功能 、 多 方法 仪器 ， 使 国内 上 自 成 
体系 的 方法 研究 、 装 备 研 发 以 及 技术 发 展 受 到 制约 。 与 发 达 国 家 相 比 ， 瞬 变 电 
磁 法 在 我 国 的 应 用 方面 与 西方 国家 大 体 相当 ， 其 他 方面 仅 相 当 于 国外 20 世纪 
90 年 代 的 水 平 ， 而 复杂 地 电 条 件 下 的 正 、 反 演 解 释 国 内 可 以 说 尚 处 于 入 门 阶 
段 。 何 继 善 (1997) 院士 指出 ， 电 磁 法 二 维 、 三 维 反 演 方 法 和 成 像 技术 将 是 理 
论 研 究 中 的 重要 方向 。 就 目前 瞬 变 电磁 法 的 理论 与 应 用 ， 仍 面临 以 下 几 方 面 的 
问题 。 l 

(1) 在 信号 采集 方面 : RAAN., BERE h RERI E — K aH E Hb 
下 介质 感应 产生 的 随时 间 变 化 的 二 次 场 。 其 特点 是 : 信号 动态 范围 大 ， 观 测 
信和 号 的 幅 值 及 其 衰减 快慢 取决 于 探测 的 地 下 介质 〈 地 质 体 ) 的 时 间 常 数值 ， 
对 于 良 导 体 ， 信 号 衰减 较 慢 ， 时 间 范 围 一 般 从 几 十 微 秒 到 几 百 毫秒 ， 信 和 号 的 幅 
值 一 般 从 数 十 或 上 百 毫 伏 变 到 零点 几 微 伏 ; 名 感应 信和 号 的 频带 宽 ， 瞬 变 电 磁 探 
测 信和 号 的 频谱 成 分 十 分 丰富 ， 频 带 范 围 一 般 为 "一 ?2X10':Hz; 图 探测 的 信和 号 受 
干扰 噪声 影响 大 ， 瞬 变 电 磁 探测 一 般 在 晚期 观测 ， 晚 期 探测 信号 一 般 为 微 伏 量 
级 ， 这 时 观测 频带 内 的 各 种 天 然 及 人 文 干 扰 品 声 相 对 较 强 ， 实 际 观测 到 的 是 各 
种 信号 的 全 加 ， 要 提高 信 品 比 ， 就 要 提高 仪器 的 灵敏 度 或 增 大 发 射 功 率 。 因 
此 ， 如 何 采 集 到 并 分 离 出 丰富 的 能 反映 地 下 介质 地 电 特 性 的 有 用 信号 ， 是 瞬 变 
电磁 勘查 能 否 成 功 的 关键 。 

目前 的 观测 技术 ， 多 数 瞬 变 电 磁 仪 紫 是 在 连续 的 等 对 数 时 间 间 隔 窗 口内 采 
集 来 获得 感应 二 次 场 的 衰减 曲线 的 ， 一 般 为 20—40 道 ， 信 息 量 十 分 有 限 。 
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GDP - 32 多 功能 电 法 系统 实现 了 算术 等 时 间 间 隔 的 采样 功能 ， 大 大 提高 了 采 
样 密度 ， 可 观测 几 百 至 上 千 道 的 数据 ， 比 常规 对 数 间 隅 采样 获得 更 丰富 的 地 电 
信息 ， 但 晚期 道 很 多 比 有 用 信息 强 的 干扰 噪声 同时 进来 了 。 如 何 获取 并 分 离 出 
含 丰 富 地 电信 息 的 有 用 信号， 仍 是 地 球 物理 工作 者 需 解决 的 问题 。 

对 于 瞬 变 电磁 探测 ， 不 论 长 偏 移 距 还 是 近 区 观测 ， 目 前 的 观测 参数 多 为 磁 
场 的 垂直 分 量 。 有 学 者 提出 应 同时 观测 水 平分 量 ， 但 鲜 有 成 功 应 用 的 文献 报 
道 。 因 此 ， 研 究 开发 多 分 量 、 阵 列 式 的 观测 技术 对 促进 瞬 变 电磁 探测 的 应 用 发 
展 具 有 重要 的 现实 意义 。 

(2) 解释 理论 方面 : 国内 外 有 关 瞬 变 电 磁 法 解释 理论 方面 的 文献 报道 很 
多 ,主要 有 浮动 薄板 解释 法 ( 视 纵 向 电导 解释 法 )、 烟 圈 理 论 解 释 法 、 瞬 变 电 
磁 时 - 频 转换 反 演 ， 一 维 、 二 维 正 、 反 演 人 机 联合 解释 法 等 。 

以 往 解释 所 用 信号 和 参数 单一 ， 尚 未 见 有 多 分 量 、 矢 量 合成 、 多 参数 联合 
解释 方面 的 报道 ， 复 杂 条 件 下 的 二 维 、 三 维 正 反 演 解 释 还 不 成 熟 ， 尚 处 于 研究 
探索 阶段 。 因 此 ， 探 索 研究 瞬 变 电磁 多 分 量 快速 解释 方法 和 合成 矢量 解释 方 
法 ， 研 究 开发 阵列 式 观测 资料 的 处 理 方 法 和 快速 适用 的 二 维 、 三 维 正 反 演 解释 
理论 的 研发 ， 将 极 大 地 提高 瞬 变 电磁 探测 的 分 辨 能 力 及 解释 精度 ， 丰 富 瞬 变 电 
磁 探 测 的 解释 理论 。 

(3) 方法 实用 性 方面 : 瞬 变 电磁 法 应 用 领域 广泛 ， 受 人 们 重视 、 关 注 的 程 
度 越 来 越 大 。 尽 管 该 方法 有 独特 的 优点 ， 但 在 某 些 特殊 地 质地 貌 开 展 工作 仍 会 
遇 到 各 种 困难 ， 例 如 ， 地 形 切 割 严 重 的 秦岭 山区 开展 金属 矿 瞬 变 电 磁 探测 ， 工 
作 布 置 困 难 ， 采 样 难度 大 ， 由 于 地 形 恶 劣 、 植 被 发 育 ， 工 作 甚 至 不 能 开展 ; Wg 
洋 瞬 变 电 磁 工作 ， 虽 有 学 者 在 从 事实 验 研 究 ， 但 仍 沿用 传统 的 装置 方式 、 观 测 
技术 ， 效 率 低 、 工 作 难 度 大 。 另 外 ， 随 着 找 矿 难度 的 增 大 ， 需 要 在 干扰 大 的 矿 
区 及 外 围 开 展 大 深度 的 勘查 工作 ; 工程 勘察 方面 往往 需要 解决 精细 构造 、 复 杂 
断面 等 方面 的 问题 ， 这 就 给 地 球 物理 工作 者 提出 了 更 高 的 要 求 ， 迫 切 需 要 研发 
抗 干扰 能 力 强 的 、 勘 探 深 度 大 的 、 分 辨 能 力 强 的 、 解 释 精 度 高 的 方法 技术 。 因 
此 ， 开 发 经 济 有 效 的 瞬 变 电磁 探测 工作 装置 、 快 速 观测 技术 、 精 度 高 的 解释 理 
论 是 瞬 变 电磁 法 探测 应 用 研究 的 一 个 方向 。 

地 震 勘 探 共 深度 点 多 次 覆盖 技术 为 我 们 提供 了 研究 思路 ， 开 发 应 用 阵列 式 
多 分 量 观 测 方法 与 技术 将 快速 获得 更 丰富 的 地 电信 息 ， 模 仿 地 震 勘 探 共 深度 点 
合 加 法 ， 人 研究 阵列 式 多 分 量 观测 资料 的 快速 适用 的 解释 方法 将 大 大 提高 瞬 变 电 
磁 探 测 的 解释 精度 ， 丰 富 瞬 变 电 磁 探测 的 解释 理论 ， 从 而 对 瞬 变 电磁 (航空 和 
海洋 瞬 变 电磁 ) 的 应 用 发 展 带 来 更 加 广阔 的 前 景 。 

因此 ， 可 以 将 瞬 变 电磁 法 的 发 展 方向 概括 为 以 下 几 个 方面 。 

(1) 理论 研究 方面 : 在 注意 发 展 研 究 复 杂 地 电 条 件 下 二 、 三 维 问题 正 反 演 
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的 同时 ， 更 应 注意 实际 应 用 效果 ， 瞬 变 电 磁 拟 地 震 的 偏 移 成 像 报 术 、 瞬 变 电 磁 
法 各 种 干扰 信号 的 分 离 技 术 人 研究 和 解释 方法 也 是 一 个 值得 注意 的 研究 方向 。 当 
前 国内 、 外 学 者 正在 寻找 一 个 数学 上 的 处 理 办 法 ， 将 瞬 变 电磁 场 变 换 成 波 场 ， 
即将 电磁 啊 应 中 与 传播 有 关 的 特征 提取 出 来 ， 而 把 电磁 波 传播 过 程 中 与 频 散 、 
诅 减 有 关 的 特征 进行 压制 或 去 除 。 如 果 能 够 实现 ， 就 能 将 瞬 变 电磁 场 的 求解 
问题 苇 化 为 大 家 所 熟悉 的 波动 方程 的 求解 问题 ， 就 能 将 地 震 偶 移 成 像 技术 、 
Born 近似 反 演技 术 、 脉 冲 谱 反 演技 术 、 层 析 成 像 反 演技 术 等 用 于 瞬 变 电磁 
场 的 反 演 解释 中 。 基 于 电磁 波 在 导电 介质 中 的 传播 和 地 震波 在 弹性 介质 中 传 
播 的 运动 的 相似 性 建立 起 来 的 反 演 解释 方法 这 些 年 来 有 了 很 大 发 展 ， 取 得 了 
里 有 成 效 的 研究 成 末 ， 在 一 定 程度 上 缩小 了 解 的 非 唯一 问题 ， 扩 大 了 电磁 法 
的 应 用 范围 ， 并 可 希望 从 观测 资料 中 提取 常规 反 演 方 法 得 不 到 的 更 多 、 更 准 
确 的 地 下 信息 。 

万 外 ， 经 典 瞬 变 电 磁 理论 将 磁性 源 和 电 性 源 看 做 偶 极 子 ， 利 用 稳 恒 电流 场 
的 磁 偶 极 子 或 静电 场 的 电 偶 极 子 公 式 ， 比 拟 到 谐 变 场 的 频 域 表达 式 ， 再 经 频 域 
到 时 域 的 变换 得 到 解析 公式 。 这 种 基于 偶 极 于 的 理论 公式 在 瞬 变 电厂 法 的 发 展 
中 起 到 了 重要 作用 ， 但 还 不 能 完全 反映 场 的 特性 、 不 能 适应 精确 勘探 的 需要 。 
可 以 用 时 变 挟 电 三 载 流 微 线 元 代 蔡 侦 极 子 第 元 ， 可 以 不 再 经 过 Fourier 变换 或 
Laplace 变换 ， 直 接 在 时 间 域 中 推导 层 状 介质 表面 上 大 回 线 和 长 接地 导线 源 的 
瞬 变 电磁 场 解 析 解 。 通 过 研究 时 域 辅助 函数 的 选取 、 时 变 点 电荷 载 流 微 元 的 处 
理 、 仿 特殊 冰 数 积分 算法 ， 分 析 瞬 变 场 在 典型 地 层 和 不 同 场 区 情况 下 的 啊 应 特 
征 等 ， 建 立新 的 电磁 场 理 论 体系 。 这 一 人 研究 方向 将 突破 长 期 沿用 的 偶 极 子 理 
论 ， 用 真正 的 微 元 代替 偶 极 子 微 元 ， 减 小 非 偶 极 子 和 时 域 频 域 转换 误差 ， 可 以 
很 好 地 反映 全 场 区 的 电磁 特性 ， 为 瞬 变 电磁 法 的 进一步 发 展 和 实际 勘探 提供 研 
究 思 路 。 

(2) 方法 拉 术 方面 : 瞬 变 电磁 法 拟 地 震 的 多 次 覆 兹 技术 人 研究 是 十 分 重要 
的 ， 为 外 ， 多 分 量 、 阵 列 式 观测 技术 研究 也 不 应 忽视 。 

针对 目前 较 广泛 使 用 的 瞬 变 电磁 回 线 源 汶 置 ， 开 展 高 分 辨 瞬 变 电 磁 法 全 域 
探测 技术 研究， 或 者 开展 对 以 往 方 法 技术 的 改进 人 研究， 将 会 大 大 提高 瞬 变 电磁 
探测 的 纵向 和 横 疝 分 辨 能 力 ， 从 而 提高 解决 地 质问 题 的 精度 ， 为 精细 工程 地 质 
勘察 提供 方法 、 技 术 上 的 支撑 。 

(3) NADE: 主要 是 发 展 多 通道 、 大 功率 、 多 功能 、 智 能 化 电 测 系统 ， 
高 灵敏 度 、 高 性 能 探头 的 研制 是 现代 电磁 仪 研制 开发 的 难点 。 受 市 场 和 国外 仪 
货 的 冲击 ， 国 内 仪 硕 需要 问 高 性 能 、 高 分 辩 方 向 发 展 ， 万 是 当务之急 。 
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1.3 了 瞬 变 电磁 法 拟 地 震 解 释 研 究 的 前 沿 课题 


1.3.1 基于 时 - 频 等 效 转换 的 瞬 变 电磁 成 像 技术 


近年 来 ， 利 用 瞬 变 电 磁 场 勘探 石油 、 地 热源 和 各 种 矿产 资源 的 理论 和 应 用 
研究 工作 在 不 断 发 展 ， 对 于 探测 埋 在 地 下 的 低 阻 异常 体 ， 瞬 变 电 磁 法 已 被 证 明 
是 一 种 有 效 的 方法 。 但 是 ， 由 于 在 分 层 、 有 耗 媒 质 中 电磁 现象 的 复杂 性 ， 目 前 
对 实测 数据 的 解释 水 平 仍 很 低 。 近 年 来 人 们 正 深入 研究 二 维 、 三 维 的 复杂 模 
型 ， 力 图 更 准确 地 描述 大 地 中 的 电磁 瞬 变 现象 ， 并 设法 从 电磁 响应 中 获得 地 下 
结构 的 局 部 形状 和 尺寸 等 高 分 辨 信息 。 

中 心 回 线 瞬 变 电 磁 法 是 一 种 近 区 观测 的 电磁 探测 方法 。 由 于 场 源 的 特殊 
性 ， 在 介质 中 传播 的 电磁 场 是 扩散 场 。 在 研究 回 线 源 瞬 变 电 磁 法 对 地 成 像 时 ， 
这 一 问题 是 不 容 回避 的 问题 ， 这 样 ， 由 扩散 场 向 平面 波 场 数据 转换 就 成 了 问题 
研究 的 关键 。 

郭 文 波 等 (2005) 从 大 量 模型 计算 入 手 ， 通 过 对 两 种 场 源 测 深 正 演 数 据 的 
分 析 、 对 比 ， 以 及 对 两 种 场 在 地 下 介质 中 传播 的 特性 分 析 ， 建 立 了 一 种 从 瞬 变 
电磁 测 深 数据 向 平面 波 场 转换 的 时 间 - 频 率 对 应 关系 。 

通过 大 量 的 理论 模型 正 演 计算 、 曲 线 对 比 、 误 差分 析 ， 结 合 场 的 特性 分 
析 ， 从 两 种 场 的 穿 透 深度 及 反映 地 下 电 性 结构 一 致 性 角度 ， 经 过 详细 推导 ， 得 
出 瞬 变 电磁 测 深 视 电阻 率 数 据 可 以 转换 成 平面 波 场 视 电阻 率 数据 的 结论 ， 构 造 
了 由 时 间 到 频率 的 转换 关系 式 

210/f =t (1. 1) 

AP, f HHR, t 为 观测 时 间 延 迟 ， 这 一 关系 可 以 从 两 种 场 在 地 下 介质 中 传 

播 的 时 间 、 速 度 、 深 度 等 方面 加 以 说 明 。 电 磁场 在 地 下 介质 中 的 扩散 速度 与 地 

下 介质 的 电阻 率 及 时 间 有 关系 。 在 阶 跃 波 断 开 后 的 某 一 时 间 延 迟 内 ， 扩 散 深 度 
与 速度 的 关系 如 下 

DO) = | Valotar (1.2) 


式 中 ，D G) 是 扩散 深度 ，Va (o, t) 是 扩散 速度 ，p 是 介质 的 电阻 率 ,，z 是 
断 电 后 的 时 间 延 迟 。 
为 了 推导 瞬 变 电磁 场 在 地 下 的 传播 速度 及 次 度 ， 我 们 分 析 肯 变 电 磁 场 在 地 
下 的 建 场 机 制 及 过 程 。 在 场 源 激励 下 ， 地 下 介质 中 产生 涡流 ， 在 阶 帐 波 断 开 
后 ， 疝 流 不 会 立即 消失 ， 而 是 有 一 个 过 滤 过 程 ， 在 这 个 过 程 中 ， 由 地 下 介质 所 
产生 的 二 次 感应 场 经 过 了 一 个 由 无 到 有 ， 由 小 到 大 、 到 极 大 、 到 衰减 、 再 到 无 
eTa 
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的 过 程 。 对 于 地 下 某 一 深度 ， 磁 场 微分 量 初始 值 为 零 ， 经 过 建 场 ， 到 达 极 大 
值 ， 最 后 误 减 到 零 。 某 一 频率 或 者 某 一 时 间 的 电磁 场 分布 在 地 下 的 任何 深度 
内 。 从 地 球 物 理 勘 探 的 角度 看 ， 电 磁 探 测 深 度 与 仪 右 的 检测 灵敏 度 、 地 电 情 
况 、 围 岩 情况 、 品 声 电 平等 因素 有 关 。 在 比较 理想 的 情况 下 ， 也 可 能 探测 到 埋 
深 相 当 于 几 倍 的 趋 肤 深度 (或 者 扩散 深度 ) 的 地 下 地 质 体 ， 在 复杂 地 质 情 况 
下 ， 也 可 能 探测 不 到 埋藏 在 趋 肤 深度 〈 或 者 扩散 深度 ) 以 内 的 地 质 体 。 但 是 ， 
总 的 来 说 ， 对 于 同一 介质 的 同一 深度 ， 扩 散场 的 视 电 阻 率 与 平面 波 场 的 视 电 阻 
率 对 此 深度 地 电 性 结构 应 该 有 相同 的 反映 。 

在 一 维 近似 的 情况 下 ， 把 趋 肤 深度 和 扩散 深度 等 效 认 为 是 电磁 场 的 探测 次 
度 。 这 样 ， 令 两 个 深度 相等 ， 从 而 可 以 找到 频率 与 时 间 的 类 似 式 (1. 1) 的 关 
系 式 。 

薛 国 强 等 〈(2004) 提出 了 瞬 变 电磁 法 对 地 成 像 理 论 和 成 像 数 值 计算 方法 ， 
其 原理 是 根据 所 构造 的 扩散 场 向 平面 波 场 转换 时 的 时 间 - 频 率 关 系 式 ， 把 时 间 
转换 成 频率 ， 这 样 就 完成 了 扩散 场 数据 癌 平 面 波 场 的 转换 。 由 转换 后 的 视 电 阻 
率 值 求 出 平面 波 场 的 波 阻 抗 ， 再 由 计算 出 来 的 波 阻 抗 ， 建 立 关 于 电 性 层 反 射 系 
数 序列 的 方程 组 ， 再 通过 线性 规划 法 ， 求 出 反射 系数 序列 ， 最 终 以 反射 系数 为 
参数 绘制 成 像 剖 面 。 

该 方法 给 出 了 瞬 变 电磁 法 对 地 成 像 数值 计算 步骤 : 中 给 出 各 种 不 同 地 电 模 
型 ， 分 别 进行 正 演 计算 ， 计 算出 视 电 阻 率 值 ， 根 据 不同 的 地 表 电 阻 率 ， 对 早期 
道 数据 进行 校正 ; 名 通过 转换 关系 ， 把 时 间 延 迟 变 成 频率 ， 对 理论 模型 正 演 数 
据 或 者 实测 视 电 阻 率 值 进行 域 的 变换 ， 把 时 间 域 扩散 场 视 电 阻 率 值 变 成 平面 波 
场 视 电 阻 率 值 ，@@ 由 平面 波 场 视 电 阻 率 值 在 频率 域 求 出 波 阻 抗 ; 由 以 波 阻 抗 为 
参数 ， 构 建 方程 组 ; @ 用 线性 规划 法 求 出 反射 系数 序列 ;，@@ 最 终 以 反射 系数 为 
参数 进行 成 像 。 


1.3.2 基于 波 场 转换 的 阳 变 电磁 成 像 技 术 


由 于 瞬 变 电磁 场 满 足 的 微分 方程 事实 上 是 一 个 扩散 方程 ， 因 而 不 能 采用 目 
前 大 家 熟悉 的 波动 方程 求解 方法 。 所 谓 瞬 变 电 磁 场 的 波 场 变换 是 指 通过 数学 积 
分 变换 ， 将 满足 扩散 方程 的 时 域 瞬 变 电磁 场 转换 为 满足 波动 方程 的 波 场 ， 然 后 
借助 于 地 震中 发 展 起 来 的 一 些 比较 成 熟 的 成 像 方法 技术 ， 求 解 被 探 目标 体 的 物 
性 和 几何 参数 。 

Weidelt (1972). Kunetz (1972). Levy 等 (1988)、Lee 和 Memechan 
(1987) 等 于 20 世纪 70 年 代 和 80 年 代 的 一 些 研究 成 果 ， 都 分 别 揭示 了 在 层 状 
大 地 介质 中 ， 电 磁 扩散 方程 与 地 震波 动 方程 间 存 在 有 趣 的 数学 对 应 形式 ， 但 他 
们 研究 问题 的 着 眼 点 都 是 将 对 应 地 电 模 型 的 波 场 模拟 结果 变换 成 时 域 电 磁 啊 
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应 。 但 是 ， 更 能 激 起 学 者 们 研究 兴趣 的 是 波 场 逆 变 换 ， 将 已 知 时 域 场 转换 为 波 
场 ， 这 将 有 利于 偏 移 以 及 更 加 复 洒 的 成 像 技 术 的 应 用 。 
从 波 场 到 时 域 场 的 波 场 正 变换 式 
H.(t) = es], se td (1. 3) 
AP., H, G) 表示 瞬 变 电磁 场 量 ,，U (z) 表示 虚拟 波动 参量 ，r 表示 虚拟 
时 间 。 

这 一 变换 过 程 称 为 正 问题 (direct problems)。 如 果 反 过 来 ， 已 知 时 域 场 
求 流 场 ， 则 称 为 反问 题 (inverse problems), 

众所周知 ， 反 问题 往往 和 不 适 定性 (ill-posed) 紧密 相 联 ， 各 种 各 样 的 反 
问题 ， 不 仅 出 现 于 地 球 物理 问题 中 ， 而 且 出 现 于 数学 本 身 。 波 场 正 变换 式 是 典 
型 的 第 一 类 Fredholm 型 算 子 方程 ， 而 它 有 一 个 重要 的 特征 ， 就 是 所 谓 的 “不 
WETE”. XF “Ww” (well-posed) 和 “不 适 定 ” (ill-posed) 的 概念 是 
Hadamard 为 了 描述 数学 物理 问题 与 定 解 条 件 的 合理 搭配 ， 于 20 世纪 初 引 
人 的 。 

瞬 变 电磁 场所 满足 的 扩散 方程 主要 刻画 电磁 涡流 场 的 感应 扩散 特征 ， 而 基 
于 扩散 方程 的 偏 移 成 像 及 反 演 方法 ， 一 般 对 电 性 界面 的 分 辨 能 力 较 差 ， 这 就 需 
要 寻找 到 一 个 数学 上 的 处 理 办 法 ， 将 瞬 变 电磁 场 变换 成 “ 波 场 "， 即 提取 出 电 
磁 啊 应 中 与 传播 有 关 的 特征 ， 压 制 或 去 除 电磁 波 传播 过 程 中 与 频 做 、 误 减 有 天 
的 特征 。 将 瞬 变 电磁 场 的 求解 问题 转化 为 波动 方程 的 求解 问题 ， 就 能 将 地 震 仿 
移 成 像 技 术 、Born 近似 反 演 技术 、 层 析 成 像 反 演技 术 等 用 于 瞬 变 电磁 场 的 反 
演 解释 中 。 

EME 〈2005b) 提出 了 从 瞬 变 电磁 场 到 波 场 的 优化 算法 。 在 波 场 的 正 变 
换 中 ， 在 保证 瞬 变 电磁 场 的 计算 精度 的 同时 ， 应 用 两 步 最 优化 算法 ， 成 功 地 控 
制 和 减少 了 积分 系数 的 个 数 和 高 散 数 字 积 分 的 采样 点 个 数 ， 解 决 了 在 波 场 反 变 
换 式 中 ， 由 于 积分 系数 过 多 而 产生 的 欠 定 方程 组 的 问题 ， 同 时 改善 了 第 一 类 算 
子 方程 的 不 适 定性 。 在 波 场 反 变换 中 ， 将 正则 化 算法 用 于 逆 变 换 过 程 的 计算 
中 ， 通 过 采用 偏差 原理 和 Newton 迭代 格式 选 出 最 优 的 正则 化 参数 ， 使 得 反 变 
换 所 得 到 的 波 场 稳定 、 可 靠 。 通 过 对 数值 计算 结果 与 已 知 波 场 限 数 对 比 ， 证 明 
了 该 方法 的 有 效 性 和 实用 性 。 大 量 的 理论 计算 表明 ， 这 种 变换 得 到 的 虚拟 波 
场 ， 不 仅 满 足 波动 方程 ， 而 且 还 类 似 于 地 震 子 波 一 样 ， 具 有 传播 、 反 射 、 透 射 

通过 理论 研究 已 经 解决 了 瞬 变 电磁 场 的 波 场 变换 问题 ， 这 就 为 瞬 变 电磁 场 
的 偏 移 成 像 创 造 了 条 件 。 但 瞬 变 电磁 场 偏 移 成 像 问题 与 地 震 勘 探 中 的 弹性 波 偶 
移 成 像 问 题 有 相似 之 处 ， 但 也 不 完全 相同 。 需 要 指出 ，TEM 虚拟 波 场 是 经 过 

. Ü. 


瞬 变 电磁 法 拟 地 震 偏 移 成 像 研究 


1.3.3 地 时 仿 移 成 像 方法 


1. Stratton-chu 积分 仿 移 成 像 方法 

Zhdanov 等 (1988)、Zhdanov 和 Booker (1993)、Zhdanov 和 Li (1997). 
Zhdanov 和 Portniaguine (1997) 信和 鉴 了 地 震 勘 探 中 的 逆 时 偏 移 概 念 ， 经 过 十 
余年 的 系统 研究 ， 对 时 间 域 的 瞬 变 电磁 场 进 行 了 逆 时 偏 移 成 像 的 深入 研究 ， 提 
出 了 偶 移 电厂 场 的 概念 ， 并 且 在 逆 时 俩 移 电 磁场 基础 上 对 二 、 三 维 反 演 问 题 也 
展开 了 最 新 的 研究 。 瞬 变 电 磁 场 的 偏 移 应 用 的 是 Stratton-Chu 积分 











] oG ` "t 
-HfS (c. ` E) VG” + (n X E) X VG” + (n X HD Z- | .dsd 
= H(r yÉ) (1. 4) 
AP, 
1/2 : 
° == ( o0) oz|r —r | n (1. 5' 
G zna Lp 537z EXP | B i TE u (t t) Cs Sy 
这 里 ， 
1 t F 
ult — t) = | wasi (1. 6) 
0 >t 


AP, G° 是 引入 的 格林 上 函数 ， 和 HH Hh P EBJHR4082, n PE T 的 外 
法 线 方 癌 的 单位 天 量 。 式 〈1.4) RHA, aE R k WW B Fa g£ 1⁄0 (8 [e] 
下 延 拓 ， 得 到 地 下 电磁 场 的 分 布 ， 延 拓 后 的 电磁 场 值 在 界面 处 发 生 突 变 ， 由 此 
可 以 确定 地 下 介质 的 地 电 构 造 。Zhdanov 等 主 要 对 孤立 和 多 个 电 性 异常 体 的 模 
型 做 了 实验 研究 ， 成 像 效 果 较 好 。 

2. 基 尔 霍 夫 (Kirchhoff) 积分 偏 移 成 像 方法 

利用 优化 算法 ， 将 瞬 变 电磁 场 进行 波 场 变换 ， 得 到 “ 波 场 "， 将 波 场 分 析 
的 原理 即 地 震波 场 从 地 面向 地 下 反 向 外 推进 行 偏 移 成 像 的 方法 用 于 对 瞬 变 电磁 
场 的 解释 ， 形 成 瞬 变 电磁 场 的 偏 移 方法 ， 也 就 是 对 观测 电磁 数据 向 下 延 拓 ， 构 
建 产生 啊 应 的 地 下 电阻 率 图 像 。 

23k (2005) 给 出 了 基 尔 霍 夫 积分 法 进行 电磁 波 场 偶 移 成 像 处 理 的 数值 计 
算 方 法 。 人 其 方程 为 








uB”. F. r= F 
上 述 方程 的 基 尔 堆 夫 积分 解 为 
Uye) =— 0 Orl a alare 
(1.7) 
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对 于 & (zx，y，z，t)， 当 把 t 变 为 一 :时 ， 既 是 w (z, y, z, t) =u (x, 
y, z, —t) BF, w (z, y, z, t) 仍 可 满足 同样 的 波动 方程 。 对 于 w(x，y， 
z, t) 是 时 间 向 前 的 问题 ， 对 于 wu (z, y, z, t) 就 是 时 间 “ 倒 退 ” 的 问题 。 
可 以 把 反射 界面 的 各 点 作为 同时 激发 上 行 波 的 源 点 。 这 样 ， 我 们 可 以 把 地 面 上 
的 接收 点 作为 二 次 发 射 源 ， 将 这 些 信息 值 时 间 “ 倒 退 ” 到 原来 状态 ， 寻 找 反 射 
界面 的 波 场 困 数 ， 以 确定 反射 界面 。 

设 自 激 自 收 的 上 行 波 为 G(xzx，y，zo，t)， 它 是 地 下 反射 界面 作为 源 点 发 
WANEH g (w, y, z, t) 在 地 面 上 z= 二 xz。 上 的 值 ， 由 式 (17) 可 得 





IL 1. l 2 L pr 2 r 
中 (工区 yt s 34 Pa = Ú F Ja Sleeper ti )dO 
+ (1. 8) 
ro 


AP, G (ë, p {os HZ) 中 的 时 间 变 量 取 是 由 于 现在 考虑 的 是 波 


动 的 道 过 程 ; 事实 上 是 由 8g (z, y, z, t) 在 地 表面 上 的 值 g (6, p Ú t) 
求 出 地 下 的 波 场 值 ， 从 而 确定 反射 界面 ， 这 即 是 波 场 的 向 下 延 拓 。 

将 三 维 边界 元 技术 引入 到 基 尔 霍 夫 积 分 计算 中 来 ， 其 最 大 的 好 处 是 能 够 实 
现 三 维 曲 面 延 拓 计 算 ， 这 在 地 震 勘 探 中 是 很 难 实现 的 ， 这 是 由 于 式 (1.8) 中 


存在 人 项， 地 震 勘探 中 无 法 进行 地 表 垂 向 梯度 测量 ， 相 反 瞬 变 电 磁 法 利用 梯 


度 探 头 可 以 很 容易 地 实现 地 表 垂 回 梯 度 测 量 。 另 外 ， 在 速度 分 析 中 ， 可 将 瞬 变 
电磁 法 中 的 反 演 技术 融 人 其 中 ， 一 方面 ， 实 现 连 续 速 度 分 析 ; 另 一 方面 ， 也 可 
将 反 演 与 偏 移 成 像 相 结合 ， 提 高 解释 的 可 靠 程度 。 

通过 对 时 间 域 的 瞬 变 电磁 场 进行 逆 时 仿 移 成 像 的 深入 人 研究， 将 瞬 变 电磁 虚 
拟 波 场 从 地 面向 地 下 反 向 外 推进 行 偏 移 成 像 ， 建 立 了 有 有 瞬 变 电磁 法 特点 的 偏 移 
成 像 方法 。 


第 2 革 BEERE 
变换 原理 


所 谓 瞬 变 电磁 场 的 波 场 变换 是 指 ， 通 过 数学 积分 变换 ， 将 满足 扩散 方程 的 
时 域 瞬 变 电磁 场 转换 为 满足 波动 方程 的 波 场 ， 然 后 借助 于 地 震中 发 展 起 来 的 一 
些 比较 成 熟 的 成 像 方法 技术 ， 求 解 被 探查 目标 体 的 物性 和 几何 参数 。 


2. 1 有 瞬 变 电 人 大。 场 与 虚拟 波 场 关系 式 的 建立 


Lee 于 1989 年 在 电 性 源 基础 上 ， 建 立 了 满足 时 域 扩散 方程 的 电场 强度 与 
虚拟 流 场 的 关系 式 ， 依 据 这 一 思路 ， 本 书 从 麦克 斯 韦 方 程 出 发 ， 建 立 以 大 回 线 
源 为 基础 的 时 域 瞬 变 啊 应 与 虚拟 波 场 的 对 应 关系 。 

在 均匀 各 回 同 性 寻 电 介质 中 ， 大 回 线 源 以 外 的 空间 中 ， 电 磁场 都 满足 麦克 
斯 韦 方程 


VX E. (r,t) 一 一 -SBa (rst) (2.1) 

VX H. (r,t) = aE. (r,t) + D.C.) (2. 2) 
Wa S (poly == 0 (2. 3) 

V. E. (r.t) = 0 (2.4) 


AP, E, (r, t) 表示 电场 强度 ，H,，(r，1) ERKE, Da (r, t) 表 
示 电 位 移 矢量 ， 且 有 
Balrst) = pHnlrst) (2. 5) 
Du (r,t) = e(r)E,,(r,t) (2. 6) 
AP, u 表示 磁 导 率 ，e 表示 介 电 常数 。 
对 式 (2.2) 两 端 求 旋 度 ， 并 将 式 (2.5), A (2.6) RA, 4 
Ə š 
Ə 





VX VX H,.(r,t) =—olr)p Z Harst) + pe (r) 
AP, p5p 为 真空 中 的 磁 导 率 。 
在 导电 介质 中 ， 如 果 忽略 位 移 电 流 ， 则 方程 C.D 可 化 简 为 
VX VX Ha (r,t) +o (r) -Hn (r,t) = 0 (2. 8) 


相应 的 初始 和 边界 条 件 可 写成 


H.G D (AT) 
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H. (r,0) = 0,H,|r 一再 (ri t>œ>0 
这 里 厂 是 体积 V fEr=r,BFBJ3 R, S[APR2K U (r, +z) 如 下 


2 
VX VX U(r,r) + o (7) AU) =0 (2.9) 


U(r,0) = Z UCr,r) l = h 
q 
U| = U(ro,z),r > 0 
对 比方 程 (2.8) 和 方程 (2.9) ， 就 可 清楚 地 看 到 自 变量 r 是 时 间 平 方 根 
EW, KHAU (r, r) 是 以 波束 一 Ç 方 传播 的 波 场 。 分 别 对 扩散 方程 (2.8) 
中 的 Hn ir T) M t zl] s 和 波动 方程 (2. 9) 中 的 U (r, r) 从 rt 到 zp 进行 La- 
place 变换 ， 得 





VX VX H(r,s) +s (r)sH(r,s) = 0 


. . 2, 10 
Hl = Hiros), — < atg(s) < > ; ) 
式 中 ， 
rin =e [Hnr Dde 
Ü 
和 
VX VX U(r,p) +o (r) p2U(r,p) = 0 
. I | 2.11 
U|r = U(ro, p); — < arg(p) < + ' ; 
其 中 ， 
U(r,p) = [Urdd 
如 果 令 s=p, MHE (2.10) 变 为 
VX VX H(n,p2) +o (r) p Hír, p’) = 0 
、 2,12) 
H|r = Hp’), — + < arg) < + . i 
将 方程 (2. 11) 与 方程 (2. 12) 两 式 相 减 ， 并 取 
D(r,p) = H(r,p2) —U(r, p) 
HD (r, p) 满足 
VX VX D(r,p) +r) Dr, p) = 0 
(2. 13) 


D|r =0,—+ < arg(p) DI 


ED (r, p) BJ3t$SEPR3KD"` (r, p) RUR (2.13) 并 在 体积 V 上 做 积 
分 ， 得 到 
. 13 ° 
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|. | VX VX Dr, p) | aV + p | or) | D(r, p) | 2 dV = Ü 


— x ELS 
1 < arg( p) < 1 


H(r,p2) = U(r, p) 
成 立 ， 去 掉 空 间 变 量 +， 并 写成 标量 形式 
[Hoedt = | Ud 


由 于 一 9 <arg (s) EET 将 s=p RAER, 48 


[HOt m | Ude 
将 此 式 从 s 到 rz 进行 Laplace 反 变 换 ， 得 到 如 下 波 场 变换 。 


] -u l 

H.G) =- N re uU C) dr (2.14) 
这 便 是 时 域 扩散 场 H (r, t) 与 虚拟 波 场 U (r, t) 之 间 的 积分 关系 表 

达 式 。 


2.2 波 场 反 变 换 式 的 不 适 定性 





2.2.1 反问 题 与 第 一 类 算 子 方程 的 不 适 定性 


在 上 节 中 ， 式 (2.14) 给 出 了 从 波 场 到 时 域 场 的 波 场 正 变 换 式 ， 这 一 变换 
过 程 称 为 正 问 题 。 如 果 反 过 来 ,已 知 时 域 场 求 流 场 ， 则 称 为 反问 题 。 波 场 正 变 
换 式 (2.14) 的 道 变换 是 典型 的 第 一 类 Fredholm 型 算 子 方程 ， 其 反问 题 有 一 
个 重要 的 特征 ， 就 是 所 谓 的 “不 适 定 性 ”。“ 适 定 ” 和 “不 适 定 ” 的 概念 是 20 
世纪 初 引 入 的 。 

设 or 和 pr 分 别 是 空间 F 和 U 的 度量 ， 算 子 A:，，F-~U 是 从 空间 下 线性 或 
非 线 性 映射 到 U。 则 式 (2.14) 中 的 反问 题 可 写成 如 下 的 第 一 类 算 子 的 形式 

Az =u (x€ F,u € U) ra. 15) 

目 然 ， 其 中 的 A 可 为 积分 算 子 。 为 此 我 们 做 如 下 定义 。 

定义 : 称 问 题 方 程 (2.15) 为 适 定 的 ， 如 果 它 同时 满足 下 述 三 个 条 件 

(1) YuEV， 都 存在 zEF 满足 方程 (2. 15) ( 解 的 存在 性 ); 

(2) B us, u CU, # > 和 zz 分 别 是 方程 (2.15〉 对 应 于 uu, 的 解 ， 
则 z1 2 ( 解 的 唯一 性 )， 
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(3) 解 相 对 于 空间 域 (F, U) 而 言 是 稳定 的 〈( 解 的 稳定 性 )， 即 Ye 守 0， 


Pumm Oe) (rt € U) (2.16) 





EA 
Prap TE (Ax, = t År: = uz) (2. 17) 
反之 ， 若 上 述 三 个 条 件 中 ， 至 少 有 一 个 不 能 满足 ， 则 称 其 为 不 适 定 性 。 

由 此 ， 我 们 知道 解 的 存在 性 与 唯一 性 取决 于 空间 (F, U) 与 算 子 A 的 
代数 特征 ， 即 算 子 是 否 为 满 射 ， 或 是 和 理 为 一 对 一 映射 〈 单 射 )， 而 稳定 性 则 
取决 于 空间 的 拓 朴 性 质 ， 即 逆 算 子 A :是否 连续 。 由 此 可 见 ， 问 题 的 适 定 与 
不 适 定 ， 不 仅 与 算 子 A 及 其 定义 域 与 值 域 有 关 ， 而 且 与 相应 空间 的 度量 

另外 ， 上 面 的 三 个 适 定性 条 件 具 有 深刻 的 实际 意义 。 首 先 ， 对 于 实际 问题 
而 言 ， 我 们 自然 期 望 其 解 是 存在 且 唯 一 的 。 更 重要 的 是 ， 实 际 获取 的 数据 资料 
都 是 近似 的 ， 即 我 们 实际 处 理 的 是 近似 数据 ， 而 不 是 “精确 ”数据 。 因 为， 
“精确 ”数据 往往 是 未 知 的 ， 若 原始 数据 的 小 的 误差 将 导致 近似 解 对 于 真 解 的 
严重 偏离 ， 则 计算 所 得 的 数值 结果 将 毫 无 意义 。 

可 见 ， 在 处 理 地球 物 理 问 题 时 ， 算 子 A 的 适 定 与 和 否 是 非常 重要 的 。 如 何 
判断 一 个 算 子 的 适 定性 ， 如 下 定理 可 以 解决 这 一 问题 。 

EM: RA: FU 为 全 连续 算 子 ( 紧 算 子 )，F,U 均 为 Banach 空间 。 
EAC FE, MF, U 之 中 一 个 是 无 穷 维 时 ， 它 必 为 无 界 算 子 ， 即 它 在 U 上 
是 不 连续 的 。 

由 上 述 定 理 可 知 ， 具 有 全 连续 算 子 的 第 一 类 积分 方程 即使 有 解 ， 也 是 不 稳 
定 的 ， 再 者 ， 当 它 的 值 域 非 闭 时 〈 由 于 观测 误差 的 存在 常常 导致 数据 资料 越 出 
该 算 子 的 值 域 )， 该 方程 无 经 典 解 。 从 而 ， 第 一 类 积分 方程 (包括 Fredholm 
方程 和 Volterra 方程 ) 一 般 是 不 适 定 的 。 

2.2.2 波 场 变换 方程 的 不 适 定性 分 析 

根据 电磁 场 理 论 中 的 对 偶 原 理 ， 瞪 变 电 磁 场 中 的 各 个 场 量 均 满 足 式 
(2.14) 的 波 场 变换 关系 ， 为 了 不 失 一 般 性 ， 取 f (xz, y, z, t) 代表 瞬 变 电 
磁场 的 场 分 量 ，w (z, yx, z, +) 代表 虚拟 波 场 U (r, z+) 革 一 分 量 ， 故 式 
(2.14) 的 离散 数值 积分 形式 可 写 为 

f(z,y,z,t,) = Pjur ys zt) a ri), (2.18) 
AP, h; 为 积分 系数 ， 并 且 

















w. 
altist) = te (2. 19) 
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为 核 函数 ， 是 随 着 r 的 增加 达到 某 一 极 值 后 快速 衰减 到 零 的 图 数 〈 图 2.1， 见 
彩 图 )， 对 于 不 同 的 时 刻 而 言 ， 核 函 数 曲线 的 形态 基本 一 致 ， 但 不 同时 间 核 消 
数 的 幅度 和 虚拟 时 间 > 的 动态 范围 很 大 ， 这 就 导致 由 核 消 数 构 成 的 系数 答 阵 的 
条 件数 会 很 大 ， 不 适 定性 非常 严重 ， 不 利于 波 场 反 变换 的 计算 。 





00 007 — — 0.04 0.06 0.08 





虚拟 时 间 /s 
(b) 
图 2.1 A Eat pR ç ñ P #p El 
(a) ! 一 1 一 100ms 核 函 数 展 布 图 ; (b) ¿=0.1—10s Ë AUE # E 


我 们 对 比 一 些 常 规 的 离散 方式 ， 如 等 间距 离散 、 对 数 等 间距 离散 、 等 面积 
离散 和 高 度 等 间距 离散 四 种 方式 ， 不 同 离散 方式 的 虚拟 时 间 分 布 示意 图 如 图 
2. 2 所 示 。 根 据 不 同 离散 方式 选取 虚拟 时 间 ， 采 用 数值 积分 方法 进行 积分 ， 比 
较 各 数值 积分 精度 和 离散 后 各 矩阵 条 件数 如 表 2. 1 所 示 。 通 过 比较 选取 均 方 误 
差 最 小 ， 形 成 系数 和 矩阵 条 件数 最 好 的 是 等 间距 离散 。 


第 2 章 瞬 变 电磁 场 的 全 域 波 场 变换 原理 





0.35 0.4 0.45 0.5 


0.25 0.3 


0.2 


0.15 


0.05 0.1 


0 


[a]/sz 


虚拟 时 


(a) 


0.07 


um: um a — —. s um mm =—=-— — — — J .— uz ma: a a wan” am aan am man tm “m a 8 m° ma —— — =— = = 





0.06 
0.05 
0.04 
0.03 
0.02 
0.01 f 


0.45 


虚拟 时 间 /s7 


(b) 
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(a) 等 面积 离散 ; (b) 对 数 等 间距 离散 ，; 
(c) 高 度 等 间距 离散 ;(d) 线性 等 间距 离散 
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表 2.1 不 同 离散 方式 比较 


离散 方式 均 方 误差 /2% 系数 矩阵 条 件数 cond (A) 
线性 等 间距 .53 7. 58 X10 
xC S [a] EB 4. 92 1. 05x10” 
等 面积 离散 6. 94 3. 25X10” 
等 高 度 离散 8.74 1. 53X10" 


2.3 波 场 变 换 的 数值 化 万 法 


由 2. 2. 1 小 节 知 ， 式 (2.14) 的 一 般 表达 式 可 写 为 
I J” ë | 
= re U(rT,y TT) dr (2. 20) 
如 果 已 知 瞬 变 电磁 场 值 上 (xz, x, z, t) 通过 对 式 (2.20) TEER, 
就 可 求 得 “ 波 场 ” 值 w(x，y，z，7r)。 但 这 一 逆 变 换 过 程 在 数学 上 归结 为 第 
一 类 算 子 方程 问题 ， 它 的 数学 特征 是 不 适 定 的 ， 将 其 离散 化 后 得 到 的 线性 代数 
方程 组 常常 是 病态 的 ， 而 且 随 着 方程 组 阶 数 的 增加 ， 病 态 情 况 更 加 严重 。 因 
此 ， 要 求 得 稳定 、 可 靠 的 “ 波 场 ” 值 必须 在 采用 正则 化 算法 的 同时 ， 还 要 控制 
方程 组 的 阶 数 。 
2.3.1 数值 化 方法 
波 场 变换 式 (2.20) 的 离散 数值 积分 形式 可 写 为 


了 (yy = 


flt, yszsti) = 2 ulz yszsrj)altist;)hj (%. RI) 
j=1 
式 中 ， 
2 
altist) = ——;e % (2. 22) 
| 2 y mt? i 


称 为 核 函 数 ，h; HENDA HPE PR a (ti gy) 是 随 r 的 增加 达到 某 一 
极 值 后 快速 衰减 的 函数 ， 而 波 场 值 w (z, y, z, +) 一 般 为 变化 幅度 不 大 的 连 
续 函 数 ， 对 数值 积分 计算 结果 影响 不 大 ， 因 此 ， 精 确 计算 式 (2.21) 的 关键 是 
求 得 一 组 hj, T; (Q =l, Da y n), 使 其 最 大 限度 地 满足 式 (2. 21), 为 此 ， 
考虑 一 种 特殊 情况 。 

邻 (m=, y, z, r) =1, WK (2.20) 变 为 

1l _ < 
dh re “dr (2. 23) 
¿ JÜ + 


faz ,y, zs t) = 
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利用 如 下 特殊 积分 








~ _ l 
| re dr = De 
A (2.23) 可 变 为 
ee OE E 
Jf (a y, zt) 一 一 Ja 2 l. TE dr (2. 24) 
对 式 (2.24) 进行 数值 积分 ， 其 离散 表达 式 为 
Jm, ,ysz,t,) = =" = 2, aG rD)h, (2. 25) 


对 式 此 用 最 优化 法 可 以 选择 出 一 组 ， Tjo h; (j =l, 2, "t, n) 
值 ， 然 后 将 这 组 值 用 到 式 (2. 20) P, FARKI t G=1, 2, +, m), 1E 
可 求 出 一 系列 的 瞬 变 电磁 场 值 (r, y, z, ti), i=l, 2, +, m, 

2.3.2 最 优化 方法 

取 瞬 变 电 磁 场 的 采样 时 间 门 数 为 m， 波 场 变 换 的 积分 系数 个 数 为 n， 为 了 
保证 波 场 变换 式 的 系数 矩阵 满 秩 ， 取 mEn. 

我 们 的 目的 就 是 要 利用 最 优化 法 选择 出 一 组 t Ma G=1, 2, |, n) 值 ， 
最 好 地 满足 式 (2. 25)， 对 于 不 同 的 i (i 二 1]，2，…，m)， 取 mm 宇 n， 可 得 








pt (2. 26) 
式 中 ， 
a =L E A. (2. 27) 
2 A Sat: | 
与 成 矩阵 形式 
A * H=F (2. 28) 


AP, AS [aj jxw， 为 系数 矩阵 ; F= [六 ]。 为 已 知 瞬 变 场 值 向 量 ; H= 
hilh 为 积分 系数 癌 量 。 
选取 和 hj; 使 目标 函数 
o= |F—A-H ||? (2. 29) 
取 极 小 。 
rt 和 有 hh; 的 选取 可 分 为 两 步 : 
第 一 步 ， 给 定 一 组 m G=1, 2, =, n), TiEa; 为 已 知 ， 由 p 的 极 小 决 
E — H h; Sle Bs "> nia H T h; 是 线性 地 满足 目标 函数 o 的 ， 故 式 
(2.29) 取 极 小 的 充 要 条 件 是 H 满足 方程 组 
A'A-H=A! .F (2. 30) 
解 该 方程 组 ， 便 可 求 出 在 给 定 z 条 件 下 ， 使 目标 函数 o 取 极 小 的 一 组 h; 
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(J=1, 2, =, m) 值 。 
第 二 步 ， 研 究 在 求 得 一 组 hh Gj=1, 2, =, n) 值 后 ， 如 何 获得 一 组 = 
(=l, 2, =, n) 使 目标 函 数 o 取 极 小 。 设 
V, = ajh; (2. 31) 
写成 矩阵 形式 
V=A.H CE S) 
为 此 ， 将 Y 在 一 组 初始 的 rm (三 1，2，…，n) ARR Taylor 级 数 ， 并 
取 ór; 的 一 次 项 


V =Vo + 全-6r (2. 33) 
AP» 
— y») 2V _ [ 9V, = a 
V, (V: Jaos JT | Ə z; | ar Lër; J, 


VO 是 一 组 初始 rt G=1, 2, =, n) 时 的 V; 值 ，sr E c 的 微小 增 量 ， 将 
式 (2.33) 代入 目标 函数 p 中 


1F 一 Yo 一 Yir | (2. 34) 
此 时 的 目标 函数 o 是 关于 57 的 函数 ， 由 o 的 极 小 可 决定 Sr， 故 有 
(B+al) ° ór = C (2. 35) 
成 立 ， 式 中 
ƏVyl/Əəv kaw (0) 
B= (ZYV (38) c= (HW) ev 


a 为 阻尼 系数 ，T 为 单位 矩阵 。 从 式 (2.35) 中 解 出 ór, TEE HY 
的 rz 
T” =t” + ór (2. 36) 
再 以 t 作为 新 的 初 值 ， 将 站 在 * 附近 展开 成 Taylor 级 数 ， 重复 上 述 过 
程 ， 可 得 tf 的 二 次 估计 值 r'*?， 如 此 往复 直至 得 到 最 佳 的 + 值 。 


2.4 波 场 逆 变 损 的 预 条 件 共 窑 榜 层 正则 化 算法 


在 前 述 的 波 场 正 变换 中 ， 应 用 两 次 最 优化 法 ， 可 以 在 保证 计算 精度 的 前 提 
下 最 大 限度 地 减少 积分 系数 和 离散 采样 点 个 数 ， 也 就 减少 了 方程 组 的 阶 数 ， 改 
善 其 不 适 定性 。 考 虑 到 波 场 变换 离散 形式 方程 (2. 21) 中 的 系数 矩阵 条 件数 较 
大 ， 单 纯 的 正则 化 共 思 梯度 法 或 者 预 条 件 共 罗 梯 度 法 ， 都 不 能 很 好 地 解决 问 
题 。 因 此 ， 将 两 种 方法 相 结合 形成 预 条 件 正 则 化 共 斩 梯 度 法 。 将 式 〈2. 21) 5 
成 矩阵 形式 为 
«Hj 
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AU = F (2, 37) 
AP, A=[a;h; J, USL u jwxi 为 虚拟 子 波 ，F 二 [fi jmxi 为 接收 的 瞬 变 场 时 
间 信 号 。 
为 了 利用 共 轿 梯度 迁 代 ， 将 式 (2. 37) 转化 为 
ATAU = ATF (2. 38) 


只 要 4 是 列 满 秩 和 矩阵 ，4 A 就 是 对 称 正定 和 矩阵， 因此 可 以 利用 共 示 梯 度 
法 。 但 是 A A 的 条 件数 较 A 的 条 件数 更 大 ， 使 方程 的 病态 更 加 严重 。 为 了 进 
一 步 降 低 息 阵 的 条 件数 ， 改 善 方程 的 病态 程度 ， 在 进行 正则 化 共 斩 梯 度 之 前 ， 
对 系数 矩阵 进行 预 条 件 。 对 于 预 条 件 和 矩阵 的 构造 采用 超 松弛 预 条 件 法 ， 这 是 因 
为 对 称 超 松弛 预 条 件 是 一 种 较为 有 将 的 预 条 件 方法 ， 不 仅 预 条 件 因 于 容易 求 
得 ， 而 且 能 有 效 降 低 矩 阵 的 条 件数 。 


设 和 矩阵 

S 王 工 工 十 中 ID (2. 39) 

可 分 解 为 
S=M—N (2. 40) 

式 中 ， 
M= — [OK AK (K La. (2.41) 
w2 — w) 

N= S A w)K HaC)’ K C(1—e)K+aC,)] (2.42) 


AF, K. C, 和 C, 分 别 为 $ 的 对 角 元 、 下 三 角 元 和 上 三 角 元 ,ww 为 (0，2) 
内 的 参数 。 于 是 我 们 可 以 选择 如 下 的 预 条 件 矩 阵 | 
P = (KiaC) !K'! (K -+¿C,) (2. 43) 
数学 上 已 经 证 明 经 过 超 松 预 条 件 后 ， 和 矩阵 条 件数 降 为 原来 的 平方 根 。 
假设 根据 矩阵 A (v) 构 造 的 预 条 件 因子 为 M (o), 构造 新 的 方程 如 式 
(2.41) IR, EEM (o) 4(z) 接 近 单位 阵 ， 因 此 迭代 很 快 收敛 。 具 体 计 算 
过 程 如 图 2. 3 所 示 ， 其 中 kw; 为 外 层 循 环 最 大 迭代 次 数 ，e A EEH HE Ez 
法 迭代 终止 条 件 ，Lwww 为 内 层 循环 最 大 的 迭代 次 数 ，& S W Etu SEERA 
止 条 件 ，v 为 选 定 的 正则 化 参数 。 
M (vY A(vx = M (v)! (wt + F) (2. 44) 
AP, x' 为 第 上 次 迭代 的 x 值 ; x 的 初 值 x“" 选 为 单位 向 量 ; Aloa 
+ATA, 
正则 化 参数 v 的 选择 非常 重要 ， 正 则 化 参数 w(6) 使 U 在 近似 性 与 稳定 性 
之 间 进 行 优 化 选择 。Zhdanov 等 〈1988) 提出 正则 化 因子 不 断 递 减 的 自 适 应 算 
法 。 并 与 工 曲线 法 做 了 比较 ， 认 为 采用 自 适 应 算法 所 得 反 演 结果 不 逊 于 工 曲 
. 2? 。 
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线 法 所 得 ， 由 于 工 曲线 法 需要 通过 多 次 反 演 计算 来 确定 最 佳 的 正则 化 因子 ， 
计算 量 将 成 倍增 长 ， 而 自 适应 算法 只 需 在 每 次 迭代 前 确定 参与 当 次 反 演 的 正则 
化 因子 ， 因 此 可 以 大 大 减少 计算 时 间 。 王 彦 飞 (2007) 依据 偏差 原理 提出 了 重 
FRH (RSCG) 计算 最 佳 正则 化 因子 ， 并 与 CGNR 和 CGER 方法 
进行 比较 ， 认 为 RSCG 方法 更 加 稳定 。 






G9<slbl, 或 k> kun aun > VE I<L. 


~ JES BEF . 
图 2.3 MREMA LH E 


依照 王 彦 飞 〈2007) 3UTR83ESSEEE (RSCG) 计算 方法 ， 根 据 部 分 先 验 
信息 来 选择 最 优 的 正则 化 参数 值 。 正 则 化 因子 o 为 一 渐变 的 量 ， 其 初始 值 v 
为 数据 拟 合 泛 函 与 稳定 泛 函 的 比值 ， 在 前 期 大 量 模拟 计算 的 基础 上 基本 可 以 确 
定 正则 化 参数 的 初 值 ， 根 据 经 验 可 以 取 o 等 于 0. 00005 为 分 段 最 大 的 正则 化 
因子 。 在 此 后 的 迭代 过 程 中 ， 如 果 数 据 拟 合 残 差 随和 迭代 次 数 逐 闻 变 小 ， 正 则 化 
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因子 可 保持 不 变 ， 否 则 按照 式 (2.45) 进行 选择 

U= wwe, 天 一 0 (2. 45) 
RP, 1 5228063828, k ERARE IE (PRCG) 的 迭代 过 程 
中 第 次 迭代 。 

在 完成 积分 离散 以 及 正则 化 参数 选择 后 ， 对 方程 (2.21) 运用 预 条 件 正 则 

化 共 斩 梯 度 法 进行 求解 计算 。 取 波 场 理 论 值 为 & (z，y，z，r) 王 1， 这 时 方 
程 为 一 个 简单 积分 ， 求 得 函数 值 C, y, z, t) = 二 1/Vxt。 正 则 化 参数 结果 
如 表 2. 2 所 示 ， 为 适用 时 间 区 间 CO. 000078，0. 006280] 的 正则 化 参数 。 利 用 
(PRCG) 法 求 得 反 变 换 结果 如 图 2. 4 所 示 ， 最 大 误差 小 于 2%， 可 知 文中 所 述 
方法 满足 要 求 。 

表 2.2 正则 化 参数 估计 结果 


正则 参数 〈v) AREN SJA 
0. 0011285 308226 2. 006% 
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图 2.4 波 场 反 变换 结果 
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我 们 知道 ， 瞬 变 场 满足 的 电磁 扩散 方程 主要 刻画 电磁 涡流 场 的 感应 扩散 特 
征 ， 体 积 效应 较 强 ， 因 此 ， 基 于 扩散 方程 的 偏 移 成 像 方法 一 般 对 电 性 乔 面 的 分 
辨 能 力 较 差 。 理 论 研 究 表 明 ， 将 瞬 变 场 转换 成 虚拟 波 场 ， 就 是 把 电磁 啊 应 中 与 
传播 有 关 的 特征 提取 出 来 ， 而 把 电磁 波 传 播 过 程 中 与 频 散 、 衰 减 有 关 的 特性 进 
行 压 制 或 去 除 。 大 量 的 理论 计算 和 实验 表明 ， 这 种 变换 得 到 的 虚拟 波 场 ， 不 仅 
满足 波动 方程 ， 而 且 还 类 似 于 地 震 子 波 一 样 ， 具 有 传播 、 反 射 、 透 射 特征 。 

利用 波 场 变 换 方法 分 别 对 两 层 模型 和 三 层 模型 啊 应 进行 波 场 变换 分 析 ， 给 
定 不 同 的 虚拟 时 间 和 子 波 宽度 ， 总 结 虚拟 波 场 的 反射 、 透 射 、 吸 收 衰减 等 传播 
规律 。 这 对 研究 虚拟 波 场 对 电 性 界面 和 电阻 率 差异 等 的 最 小 分 辩 率 具有 重要 
意义 。 


3. 1 MERENDA 


根据 特殊 子 波 u Ga, y, z, z) =1 可 以 确定 正则 化 参数 的 合理 范围 。 选 择 
正则 化 参数 后 ， 以 高 斯 脉冲 为 子 波 合成 虚拟 波 场 记录 ， 对 地 电 模 型 进行 模拟 。 


3.1.1 单 正 峰 波 场记 录 模 型 


以 高 斯 脉冲 为 子 波 合成 波 场记 录 ， 如 图 3.1 (a) 所 示 。 正 变换 后 模型 衰 
减 曲 线 如 图 3. 1 (b) 所 示 ， 图 中 虚线 为 特殊 子 波 u (z) = 1 时 的 衰减 曲线 。 
丢 加 子 波 的 积分 函数 如 图 3.1 (c) 所 示 。 误 减 曲 线形 态 与 D 型 模型 正 演 的 豪 
减 曲线 变化 规律 一 致 ， 可 假设 单 峰 记 录 为 D 型 地 电 模 型 波 场 图 数 。 分 析 可 知 ， 
脉冲 峰值 时 刻 对 应 着 电 性 变化 的 分 界面 ， 高 斯 脉冲 的 宽度 对 应 着 电阻 率 差异 ， 
幅 值 符号 对 应 地 电 模 型 的 高 低 电 阻 率 相 对 变化 。 按 照 第 2 章 所 述 方法 选择 最 优 
正则 化 参数 后 ， 进 行 波 场 反 变换 。 

给 出 系列 的 单 峰 虚拟 合成 记录 。 固 定子 波 宽度 (W) 和 幅 值 (A), WE 
峰值 出 现时 刻 (T,)， 范 围 为 0.02 一 0. 1si1 。 对 给 定 模型 记录 进行 反 演 拟 合 ， 
采样 时 间 为 10ks 一 10ms， 虚 拟 时 间 范 围 为 0.005~0. 5s? ， 离 散 点 数 为 M==NN 
二 101， 最 优 正则 化 参数 为 0. 003。 反 演 结果 如 图 3.2 所 示 ， 图 中 实 线 为 合成 
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模型 记录 U(r)， 空 虚线 为 反 演 拟 合 结果 U (z). 






0 0.1 0.2 0.3 0.4 0.5 
虚拟 时 间 /sz 
(a) 


menso 半空 间 模型 训 减 曲线 
D 型 模型 衰减 曲线 
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RT [B]/s 
(b) 
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HEB B]/s2 
(c) 
图 3.1 利用 子 波 合成 DD 型 模型 波 场 记录 
(a) 虚拟 波 场 记录 ; (b) 等 效 模型 衰减 曲线 ;(c) 不 同 采样 时 刻 亚 加 子 波 的 核 函 数 


. 26 ° 


第 3 章 瞬 变 电磁 虚拟 波 场 的 波 场 特性 分 析 





“00 0.1 0.2 0.3 04 05 
| 
虚拟 时 间 /sz 
(a) 








“0.0 0.1 0.2 0.3 04 05 
| 
虚拟 时 间 /s= 
(b) 
L.7 
1.6 
1.5 | 
à U' (可 
1.4 ° U (z) 
1i T,=0.04 
-FE | 
12 | 
1.1 
10 wS, 
0.9 L° es 
0.0 0.1 0.2 0.3 0.4 0.5 
l 
虚拟 时 间 /sz 
(c) 


图 3.2 不 同 虚拟 时 间 D 型 模型 波 场 变 换 结 果 图 
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图 3.2 不 同 虚 拟 时 间 D 型 模型 波 场 变换 结果 图 ( 续 ) 
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图 3.2 不 同 虚拟 时 间 D 型 模型 波 场 变换 结果 图 (SE) 
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由 图 3. 2 可 知 ， 随 着 峰值 时 间 的 推迟 ， 反 演 子 波幅 值 逐 渐 降 低 ， 且 子 波 宽 
度 略 有 展 宽 。 统 计 计 算 结 果 ， 如 表 3. 1 所 示 ， 随 着 虚拟 时 间 向 后 延续 ， 反 演 误 
差 逐 渐 增 大 ， 误 差 变化 如 图 3. 3 所 示 。 这 与 电磁 波 在 有 耗 介 质 中 的 传播 规律 基 
本 相符 ， 也 是 由 此 原因 ， 使 得 瞬 变 电磁 的 分 辨 率 有 所 降低 。 

表 3.1 不 同 虚拟 时 间 模型 迭代 反 演 次 数 和 方差 表 


序号 gimis? THRE 子 波 宽度 迭代 次 数 。 ” 均 方差 /% 
a 0. 02 0. 6 0. 04 66941 0. 83 
b 0. 03 0. 6 0. 04 66203 0. 54 
C 0. 04 0.6 0. 04 60065 0. 96 
d 0. 05 0.6 0. 04 61031 L. Aa 
e 0. 06 0.6 0. 04 63133 1. 38 
Í 0. 0/ 0. 6 0. 04 63692 1.7 
g 0. 08 0.6 0. 04 60975 1.91 
h 0. 09 0. 6 0. 04 57203 2. 03 


1 0. 10 0.6 0. 04 61172 2. 41 


FJ Jy 25 /% 
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图 3.3 均 方 误差 随 峰值 出 现时 间 变 化 曲线 


固定 子 波 峰值 出 现时 间 ， 改 变 子 波 宽度 ， 对 模型 进行 反 演 拟 合 ， 反 演 结 采 
如 图 3.4 所 示 ， 图 中 实 线 为 合成 的 波 场记 录 ， 小 圆圈 为 反 演 拟 合 结 采 ， 和 采样 间 
隅 与 虚拟 时 间 间 隔 不 变 。 
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图 3.4 固定 模型 峰值 时 间 改 变 子 波 宽度 的 反 演 拟 合 图 
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图 3.4 固定 模型 峰值 时 间 改 变 子 波 宽度 的 反 演 拟 合 图 (2) 
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图 3.4 固定 模型 峰值 时 间 改 变 子 波 宽 度 的 反 演 拟 合 图 ( 续 ) 
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图 3.4 固定 模型 峰值 时 间 改 变 子 波 宽度 的 反 演 拟 合 图 ( 续 ) 


分 析 模 型 的 反 演 结果 ， 统 计 参 数 如 表 3. 2 所 示 。 由 图 3. 4 可 知 在 相同 深度 
下 ， 随 着 虚拟 子 波 宽度 逐渐 变 窗 ， 反 演 拟 合 的 振幅 与 模型 振幅 差距 逐渐 增 大 ， 
反 演 误差 增 大 ， 当 子 波 宽度 很 罕 时 〈 电 阻 率 差异 不 大 ) 反 演 拟 合 数据 与 模型 差 
异 加 大 ， 因 此 ， 不 能 分 辨 该 宽度 的 子 波 。 拟 合 误差 与 子 波 宽度 之 间 的 关系 如 图 
3. 5 所 示 ， 当 子 波 宽度 很 小 时 ， 拟 合 误差 反而 变 小 但 是 拟 合 曲线 旁 瓣 与 主 幅 值 
差距 不 大 ， 因 此 分 辨 效果 不 佳 。 
表 3.2 D 型 模型 固定 深度 参数 反 演 误差 


序号 峰值 时 间 /s + IB IB TERE 迭代 步 数 | 均 方差 /% 
a 0.2 0.6 0.20 62863 1.49 
b 0. 2 0. 6 0. 18 62373 1. 78 
c 0.2 0.6 0.16 60892 2.12 
d Ú, Š 0.6 0.14 61301 2. 57 
e 0.2 0.6 0.12 56076 3. 09 
f 0. 2 0.6 0.10 58696 3.58 
g 0. 2 0.6 0. 08 61728 3. 95 
h 0. 2 0.6 0.06 61882 4. 02 
i 0.2 0.6 0.04 60665 3.53 


] 0. 2 0.6 0. 02 62295 2. 34 
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子 波 宽度 
图 3.5 均 方 误差 随 子 波 宽度 的 变化 曲线 


3.1.2 单 负 峰 模型 分 析 


按照 对 单 正 峰 模型 进行 模拟 的 方式 对 单 负 峰 模型 进行 分 析 ， 以 负 高 斯 脉冲 
为 子 波 合成 的 虚拟 波 场记 录 如 图 3.6 (a) 所 示 ， 积 分 后 衰减 曲线 如 图 3.6 (b) 
所 示 。 友 加 子 波 的 积分 晒 数 如 图 3.6 (c) 所 示 ， 对 其 进行 离散 积分 得 到 的 总 
减 曲 线 与 正常 正 演 G 型 地 电 模 型 衰减 曲线 的 变化 规律 相符 。 因 此 ， 可 将 单 负 
峰 高 斯 脉冲 合成 记录 视 为 G 型 模型 ， 按 照 上 述 方法 进行 反 变 换 。 

首先 ， 固 定子 波 宽度 和 幅 值 ， 改 变 子 波峰 值 出 现时 间 ， 并 利用 预 条 件 正则 
CHIB BEA HITEN, MARUR 3. 7 所 示 。 
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图 3.6 利用 子 波 合成 G 型 模型 的 波 场记 录 
(a) 虚拟 波 场 记录 ; (b) 等 效 模型 衰减 曲线 ;(c) 不 同 采 样 时 刻 秋 加 子 波 的 核 沙 数 
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3.6 ”利用 子 波 合成 G 型 模型 的 波 场 记录 ( 续 ) 
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图 3.7 固定 模型 子 波 宽度 改变 峰值 时 间 参 数 的 反 演 拟 合 图 
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图 3.7 “固定 模型 子 波 宽度 改变 峰值 时 间 参数 的 反 演 拟 合 图 ( 续 ) 
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图 3.7 固定 模型 子 波 宽度 改变 峰值 时 间 参 数 的 反 演 拟 合 图 (SE) 


由 图 3. 7 可 知 ， 随 着 负 峰 值 的 延迟 ， 反 演 子 波 振幅 逐渐 衰减 ， 子 波 宽度 略 
有 增加 。 这 与 波 在 有 耗 介质 中 的 传播 规律 相符 ， 从 而 证 明了 算法 的 正确 性 。 对 
模型 反 演 参数 做 统计 ， 如 表 3. 3 所 示 。 随 着 子 波峰 值 出 现时 间 的 延迟 ， 拟 合 误 
差 逐 渐 增 大 ， 如 图 3. 8 所 示 。 模 型 的 反 演 结果 充分 说 明了 瞬 变 电磁 方法 对 浅 部 
异常 分 辩 率 较 高 。 随 着 深度 的 增加 ， 电 磁 波 高 频 部 分 损失 严重 ， 啊 应 主要 为 低 
频 响 应 ， 因 此 分 辩 率 降低 。 

固定 子 波 峰值 时 间 ， 改 变 子 波 宽 度 ， 按 照 文中 所 述 方法 进行 波 场 变换 ， 以 
找到 波 场 对 子 波 宽度 的 最 小 分 辨 。 波 场 变换 结果 如 图 3. 9 所 示 。 
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表 3.3 G 型 模型 固定 电阻 率 参数 的 反 演 误差 
子 波幅 值 


序号 ”峰值 时 间 /st 
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图 3.9 ”固定 模型 第 一 层 厚 度 改 变 电 阻 率 差 异 的 反 演 拟 合 图 
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图 3.9 固定 模型 第 一 层 厚 度 改变 电阻 率 差 异 的 反 演 拟 合 图 (22) 
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根据 反 变 换 结果 可 知 ， 在 相同 峰值 时 间 下 ， 子 波 宽度 越 大 越 容 易 分 辨 。 随 
着 子 波 宽度 的 变 小 ， 反 变换 结果 对 界面 的 显示 越 来 越 不 明显 。 对 反 演 模型 结果 
进行 列表 分 析 ， 如 表 3. 4 与 图 3. 10 所 示 ， 除 了 前 两 点 外 ， 随 着 子 波 宽度 的 增 
加 ， 拟 合 误差 逐渐 减少 ， 当 子 波 宽度 为 0.04s 时 ， 拟 合 误差 较 小 ， 而 且 虚 拟 
子 波 幅 值 比 最 大 旁 瘀 幅 值 约 大 一 个 量 级 。 当 子 波 宽度 为 0.02s? 时 ， 虽 然 拟 合 
误差 较 小 ， 但 是 虚拟 子 波幅 值 比 最 大 旁 辩 幅 值 仅 大 几 倍 ， 因 此 ， 随 着 子 波 宽度 
变 窜 ， 波 场 不 能 很 好 地 分 辨 电 性 界面 。 

表 3.4 G 型 模型 固定 模型 第 一 层 厚 度 参 数 反 演 误 差 表 





序号 峰值 时 间 /s? 子 波幅 值 子 波 宽 度 迭代 步 数 SPREA 
a 0.2 — 0. Ü 0. 20 62519 1.44 
b 0. 2 — 0. 6 0. 18 63641 L: TE 
C 0. 2 —0. 6 0.16 60894 m. l 
d 0. 2 一 0.6 0. 14 63855 2. 59 
0.2 —0. 6 0. 12 64006 3.14 
f 0.2 —0. ô 0.10 62725 3. 67 
g 0.2 —0. 6 0. 08 62404 4. 06 
h 0. 2 —0. 6 0. 06 61728 4.15 
1 0. 2 —0. 6 0. 04 62643 3. 65 
j 0. 2 —0. 6 0. 02 64018 2. 44 
4.5 
4.0} 
3.5 
g 30 
$ 
g 25 
2.0 
1.5 
1.0 
0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.22 





子 波 宽度 
图 3. 10 均 方 误差 随 了 于 波 宽 度 变 化 曲线 
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分 析 两 层 模型 的 波 场 变换 结果 ， 和 总 体 上 来 说 G 型 模型 与 D 型 模型 变换 结 
果 基 本 相似 。 在 相同 深度 下 ， 随 着 层 间 电阻 率 的 差异 减 小 ， 虚 拟 波 场子 波幅 值 
逐渐 减 小 ， 并 且 子 波 锅 度 小 到 一 定 程 度 时 ， 波 场 不 能 分 辨 介质 电 性 分 界面 。 在 
相同 电阻 率 差 异 条 件 下 ， 随 着 峰值 虚拟 时 间 的 推迟 ， 虚 拟 波 场子 波 振幅 逐渐 训 
减 ， 而 且 子 波 具 有 展 宽 现 象 。 


3.2 三 层 模型 分 析 


在 导电 妹 质 中 ， 电 磁 流 的 传播 速度 和 衰减 速度 都 较 快 ， 因 此 对 多 层 界面 的 
识别 比较 困难 。 为 了 分 析 反 变换 对 多 层 模 型 的 计算 效果 ， 研 究 人 员 设 计 了 具有 
两 个 子 流 的 模型 ， 并 进行 反 演 拟 合 分 析 。 


3.2.1 双 正 峰 虚 拟 波 场记 录 


波 场 记录 如 图 3.11 (a) 所 示 ， 正 变换 后 所 得 扩散 场 时 间 衰 减 信号 如 图 
3.11 (b) 中 实 线 所 示 ， 为 了 方便 对 比 ， 图 中 虚线 给 出 了 无 反射 界面 情况 下 豪 
减 信 号 。Q 型 模型 电阻 率 逐 层 减 小 ， 表 现在 瞬 变 电磁 衰减 曲线 上 为 衰减 信号 
降低 速度 逐渐 变 缓 ， 与 所 给 模型 扩散 信号 相似 ， 因 此 可 将 该 记录 视 为 Q 型 模 
型 虚拟 波 场 合成 记录 。 


1.8 





一 Q 型 模型 衰减 曲线 





0 0.1 0.2 0.3 0.4 I 0.5 
虚拟 时 间 /s 时 间 /s 
(a) (b) 
图 3.11 利用 子 波 合成 Q 型 模型 波 场 的 记录 
(a) 虚拟 波 场 合成 记录 ; (b) 正 变换 总 减 曲 线 


根据 文中 所 述 方 法 进行 反 演 。 为 了 便于 分 析 首 先 确定 子 波 宽度 ， 给 定 第 一 
峰值 时 刻 与 第 二 峰值 时 刻 的 间距 ， 改 变 第 一 峰值 出 现时 间 〈T)。 结 果 如 图 
3.12 (a) 一 图 3.12 (c) 所 示 。 然 后 ， 给 定 第 一 峰值 时 间 (Ti)， 了 逐渐 提前 第 
二 峰值 出 现时 间 〈T: ) ， 结 果 如 图 3.12 (c) 一 图 3. 12 (f) 所 示 。 
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图 3. 12 Q 型 模型 波 场 反 变 换 结果 
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图 3.12 QQ 型 模型 波 场 反 变 换 结果 (SE) 
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由 Q 型 模型 反 变换 结 末 可 和 项， 对 于 多 层 界 面 模型 ， 随 看 模型 峰值 出 现时 
刻 的 推迟 ， 虚 拟 子 波 的 幅 值 逐 渐 误 减 ， 而 且 第 二 个 子 波 振幅 衰减 尤为 严重 。 随 
着 两 个 峰值 出 现 间隔 的 减 小 ， 虚 拟 波 场 反射 子 波 逐渐 相互 重 丢 ， 使 得 对 薄 层 的 
分 辩 率 降低 。 由 表 3. 5 可 知 ， 当 已 知 子 波 距离 为 0. 04s 时 ， 波 场 变换 子 波 相 
互 重 到， 对 薄 层 基本 不 能 分 辨 。 给 出 模型 电阻 率 参 数 ， 由 此 可 以 确定 波 场 对 不 
同 深 度 薄 层 的 最 小 分 辨 厚度 。 

RIS Q 型 模型 波 场 变换 误差 统计 表 
第 一 峰 ”第 一 子 ”第 一 子 第 二 峰 第 二 子 第 二 子 选 代 


ii EVE WREE WEE EWE WAE WREE 步 数 all 
a 0 1 0. 06 0. 6 0. 20 0. 06 0.5 126213 4.7582 
b 0.08 0.06 0. 6 0. 18 0. 06 0.5 1277407 3.8505 
c 0. 06 0.06 0.6 0.16 0. 06 0.5 1271561 3.7838 
d 0. 06 0. 06 0.6 0.14 0. 06 05 1272517 3.4252 
Ë 0. 06 0. 06 0. 6 0. 12 0. 06 0.5 1272594 3.7434 


f 0. 06 0. 06 0.6 0. 10 0. 06 0.5 1267121 2. 8502 


3.2.2 正 负 双 峰 虚拟 波 场 记录 


图 3.13 (a) 为 虚拟 波 场 合成 记录 ， 其 在 时 间 域 扩散 场 的 衰减 信号 如 图 
3.13 (b) 所 示 。 由 于 H 型 模型 电阻 率先 减 小 后 增 大 ， 表 现在 瞬 变 电磁 衰减 曲 
线 上 为 衰减 信号 幅度 降低 速度 先 变 缓 后 逐渐 加 快 ， 与 模型 在 时 间 域 衰减 信号 的 


2 10? 
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图 3.13 合成 H 型 模型 波 场子 波 示意 图 
(a) 虚拟 波 场合 成 记录 ; (b) 正 变换 衰减 曲线 
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变化 一 致 。 因 此 ， 合 成 瞬 变 电磁 虚拟 波 场 记录 为 类 H 型 虚拟 波 场记 录 。 根 据 
文章 所 述 算法 ， 对 衰减 信号 进行 反 变换 ， 并 与 已 知 波 场 对 比 ， 反 变换 结果 如 图 
3. 14 所 示 。 

依照 文中 所 述 方法 对 H 型 模型 的 衰减 信号 进行 波 场 变换 。 首 先 固定 模型 
子 波 宽度 ， 给 定 第 一 峰值 时 间 ， 第 二 峰值 时 间 逐 步 延迟 ， 对 模型 进行 变换 ， 波 
场 变换 如 图 3. 14 所 示 。 

由 图 3. 14 (a) 一 图 3. 14 (g) 可 知 由 于 第 一 峰值 出 现时 间 不 变 ， 反 变换 
波形 中 第 一 个 子 波 与 理论 子 波 吻合 较 好 ， 随 着 两 峰值 间隔 的 增加 ， 变 换 的 第 二 
个 子 波 振幅 逐渐 训 减 ， 由 表 3.6 可 知 拟 合 误差 逐渐 增 大 。 
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(b) 
图 3.14 HH 型 模型 波 场 变 换 图 
(a) 第 二 峰值 时 间 T;=0.06; (b) 第 二 峰值 时 间 T; =0.08; (c) 第 二 峰值 时 间 T; =0.10; 
(d) 第 二 峰值 时 间 T; =0. 12; (e) 第 二 峰值 时 间 T; =0.14; 
(D 第 二 峰值 时 间 T; =0.16; (g) 第 二 峰值 时 间 T; =0. 18 
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图 3.14 H 型 模型 波 场 变换 图 (SE) 
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图 3.14 H 型 模型 波 场 变换 图 (SE) 


表 3.6 H 型 模型 波 场 变换 误差 统计 表 


第 一 子 第 一 子 第 二 峰 BT BT BR 
波 宽度 WWE EAE ” 波 宽 度 ” 波 振幅 FA 


0.06 0.6 0. 06 0.06 — Ü. 9 73491 


0. 06 0. 6 0. 08 0. 06 — 0. Š 71107 


0. 06 0. 6 0.10 0. 06 —0.5 1276072 


. 5] ° 


方差 /入 


1. 3743 
1.7725 


3. 032 


BE S R024030 8 835 pR: 


续 表 
序号 第 一 峰 ”第 一 子 第 一 子 Fe 第 二 子 ”第 二 子 ÆR 方差 /% 
值 位 置 WAE ”波幅 值 AE WAE 波 振 幅 。 步 数 
d 0.04 0.06 0.6 012 0.06 一 0.5 1273393 2.9119 
e 00 — 0.06 0.6 0.14 006 一 05 1275730 2.8213 
f 0.06 0.06 0. 6 016 0.06 一 0.5 1269651 4.3382 
g 0.08 0.06 0.6 018 006 一 0.5 1274553 5.3975 


3.2.3 负 正 双 峰 波 场记 录 


负 正 双 峰 波 场 记录 如 图 3.15 (a) 所 示 ， 其 对 应 的 时 间 域 彭 减 信号 如 图 
3.15 (b) 所 示 。 误 减 曲线 变换 趋势 与 K 型 模型 相似 ， 表 现在 瞬 变 电磁 衰减 曲 
线 上 均 为 衰减 信号 ， 首 支 衰减 较 快 尾 支 衰减 变 缓 。 因 此 ， 此 波 场记 录 可 视 为 开 
型 模型 波 场记 录 。 





-KK 型 模型 衰减 曲线 





0 0.1 0.2 0.3 0.4 0.5 
虚拟 时 间 /s7 时 间 /s 
(a) (b) 
图 3.15 合成 K 型 模型 波 场子 波 与 衰减 曲线 示意 图 
(a) 虚拟 波 场 合成 记录 ; (b) 正 变换 衰减 曲线 


依照 第 2 章 所 述 方法 对 K 型 模型 的 衰减 信号 进行 波 场 变换 。 固 定 模型 子 
波 宽 度 参 数 和 峰值 时 间 间 隔 。 第 一 峰值 时 间 逐 渐 推 惨 ， 波 场 变 换 结 来 如 图 
3.16 所 示 。 由 图 可 知 虽然 峰值 时 间 间 隔 不 变 ， 但 是 随 着 第 一 峰值 时 间 的 推迟 ， 
第 二 峰值 逐渐 延迟 。 在 反 变 换 中 ， 第 一 个 子 波 振幅 逐渐 减低 ， 第 二 个 子 波 振幅 
衰减 较 第 一 个 子 波 更 快 。 由 表 3.7 可 知 拟 合 误差 也 逐渐 增 大 。 
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图 3.16 K 型 模型 波 场 反 变换 结果 


. 53 o 


0.5 


0.5 





0.5 


瞬 变 电磁 法 拟 地 震 偏 移 成 像 研究 








“00 0.1 02 03 0.4 0.5 


虚拟 时 间 /s 
(d) 





U (z) 


| ° U (r) 
1.4 
12 Josa T=0.18 


0.0 0.1 02 0.3 0.4 0.5 
l 
虚拟 时 间 A 
(e) 





幅度 


%0 al a 0 Qi ùs 
虚拟 时 间 /sz 
(Í) 
图 3.16 KK 型 模型 波 场 反 变 换 结果 (SE) 
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图 3. 16 K 型 模型 波 场 反 变换 结果 (Sk) 
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表 3.7 KK 型 模型 波 场 变换 误差 统计 表 
n 第 一 峰 第 - 子 第 一 子 第 二 峰 第 二 子 第 二 子 ÆR | 
TS ata nae wma maw wre wrm sa JE 


a 0.04 0.06  —0.6 014 0.06 0.5 1272259 3.373 
b 0.05 006  —0.6 0.15 0.06 0.5 1273564 3.9016 
c 0.06 0.06  —0.6 0.16 0.06 0.5 1262827 4.2297 
d 007 0.06  —0.6 01 0.06 0.5 1270457 4.284 
e 0.08 006  —0.6 018 0.06 0.5 1270215 4.6335 
f 0.09 0.06  —0.6 0.19 0.06 0.5 127149 3.8511 
g 0.10 0.06  —0.6 0.20 0.06 0.5 127621 4.7129 
h 0.11 0.06  —0.6 202 0.06 0.5 1267150 4.9841 
l 


012 0.06 —06 02 0.06 0.5 1273142 5.0905 


3.2.4 双 负 峰 波 场记 录 

双人 负 峰 波 场 记录 如 图 3.17 (a) 所 示 ， 其 对 应 的 时 间 域 衰减 信号 如 图 3.17 (b) 
所 示 。 误 减 曲 线 变 换 趋势 与 A 型 模型 衰减 信号 相似 ，A 型 模型 电阻 率 逐 渐 增 
大 ， 表 现在 瞬 变 电磁 芭 减 曲线 上 为 衰减 信号 逐步 变 快 ， 与 所 给 波 场记 录 减 信号 


相似 。 因 此 ， 此 波 场 记录 可 视 为 A 型 模型 波 场 记录 。 


10 







somaar -ka hE 
— A 型 模型 衰减 曲线 














10? 
其 
yE 
10' 
10° — WI 
0 Ol 02 03 04 05 [Ga yp <= Ñ 
虚拟 时 间 /s7 时 间 /s 
(a) (b) 


图 3.17 A 型 模型 合成 波 场 记录 与 衰减 曲线 示意 图 
(a) 虚拟 波 场 合成 记录 ; (b) EWH 
依照 第 2 章 所 述 方法 对 A 型 模型 的 衰减 信号 进行 波 场 变换 。 固 定子 波 宽 
度 和 第 一 峰值 时 间 。 第 二 峰值 时 间 逐 渐 提 前 ， 波 场 变换 结果 如 图 3. 18 所 示 ， 
并 给 出 误差 统计 ， 如 表 3. 8 所 示 。 由 图 3. 18 可 知 由 于 第 一 峰值 时 间 不 变 ， 波 
场 变换 的 首 个 子 波 振幅 不 变 。 随 着 第 二 峰值 时 间 的 提前 ， 第 二 个 子 波 振幅 逐渐 
增加 。 两 个 子 波 相互 靠近 ， 当 子 波 距离 为 0.04s? 时 ， 两 个 子 波 释 加 在 一 起 。 
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因此 ， 波 场 反 变换 对 于 小 于 0. 04s? 的 薄 层 ， 不 能 分 辨 。 
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图 3.18 A 型 模型 波 场 反 变换 结果 
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图 3.18 A 型 模型 波 场 反 变换 结果 (2E) 
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表 3.8 A 型 模型 波 场 变化 误差 统计 表 


第 一 峰 SF “第 一 子 ”第 二 峰 ”第 二 子 ”第 二 子 


值 位 置 ” 波 宽 度 WRA EME WKE 波 振幅 


序号 
a 0. 04 
b 0. 04 
c 0. 04 
d 0. 04 
e 0. 04 
f 0. 04 


0. 06 —0. 6 0. 13 0. 06 — 0. 5 
0. 06 —0. 6 0. 12 0. 06 — 河 
0. 06 —0. 6 0.11 0. 06 一 局, 5 
0. 06 —0. 6 0. 10 0. 06 — 0. 5 
0. 06 —0. 6 0. 09 0. 06 — Ü. 9 
0. 06 一 0.6 0. 08 0. 06 — 0. 5 


3.3 HERRIAN AIRE EIR AR E 


野外 测量 数据 受 各 种 因素 影响 难免 会 有 误差 ， 下面 先 以 单位 向 量 为 例 ， 首 
先进 行 波 场 正 变 换 ， 得 到 时 间 域 响应 ， 然 后 在 时 间 域 响应 中 分 别 加 入 信 噪 比 为 
100:1、50:1、30:1 和 20:1 的 高 斯 日 噪声 。 其 含 噪 时 间 域 啊 应 和 相应 的 


反 变 换 结果 如 图 3. 19( 见 彩 图 ) 和 图 3. 20 ( 见 彩 图 ) 所 示 。 


102 


10' 


磁感应 强度 杂 


信和 上 曲 比 =100 : 1 


磁感应 强度 /T 
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图 3.19 U (+) =1 的 含 噪 时 间 域 响应 
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ENR 
步 数 
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1276592 
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1272992 
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1269852 


方差 /% 


0. 029219 
0. 026293 
0. 024571 
0. 023016 
0. 025678 
0. 025171 


Bi B SET E RE kR SI: 
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图 3.19 U e) =1 的 含 噪 时 间 域 响应 〈 续 ) 


“| 信 品 比 =100 : 1 
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图 3.20 U (z) =1 的 含 噪 时 间 域 响应 反 变 换 
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图 3.20 U (+) =1 的 含 品 时 间 域 响应 反 变 换 (SE) 


由 图 3. 19 (d) 可 以 明显 看 出 噪声 使 得 时 间 域 响应 曲线 变 得 很 不 平滑 。 由 
图 3. 20 的 波 场 反 变换 结果 来 看 ， 当 信 了 噪 比 为 50 : 1 时 ， 反 演 值 与 期 望 值 吻 合 
得 还 可 以 ， 当 信和 了 噪 比 降低 到 30 : 1 和 20 : 1 时 ， 反 演 值 波动 剧烈 ， 与 期 望 值 已 
相差 甚 远 。 

下 面 让 虚拟 波 依 次 取 高 斯 正 向 单 脉 冲 、 负 向 单 脉冲 、 正 向 双 脉 冲 、 负 向 双 
脉冲 、 正 负 向 双 脉 冲 和 负 正 向 双 脉 冲 六 种 不 同 的 波形 ， 含 噪 实验 结果 如 图 
3.21 ( 见 彩 图 ) 到 图 3.32 〈 见 彩 图 ) 所 示 。 

由 图 3. 22 ( 见 彩 图 ) 和 图 3. 24 可 以 看 到 ， 当 信 噪 比 高 于 30 :1 时 ， 含 
噪 时 间 域 响应 的 波 场 反 变 换 曲 线 与 期 望 曲线 吻合 得 还 可 以 ， 尤 其 是 单 峰 的 定 
位 很 准确 。 当 信 噪 比 进一步 降低 到 20 : 1 时 ， 波 场 反 变换 曲线 出 现 强烈 波 
动 ， 严 重 干 扰 了 单 峰 的 识别 。 从 图 3. 21 ( 见 彩 图 ) 和 图 3.23 可 以 看 到 信 噪 
比 为 20 : 1 的 时 间 域 响应 曲线 很 不 光滑 ,干扰 明显 ， 可 见 在 数据 处 理 中 光滑 
的 重要 性 。 
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图 3.21 正 向 单 脉冲 的 含 品 时 间 域 响应 
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图 3.22 正 向 单 脉冲 含 品 时 间 域 响应 反 变换 
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图 3. 26 正 向 双 脉 冲 含 噪 时 间 域 响应 反 变 换 


由 图 3.26 (WERD 可 以 看 到 ， 当 信 噪 比 为 30 : 1 时 ， 波 场 反 变换 曲线 的 
第 一 个 峰值 与 期 望 峰值 比较 吻合 ， 而 第 二 个 峰值 吻合 得 不 好 。 由 于 高 斯 白 噪 是 
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随机 给 的 ， 每 次 反 演 的 结果 可 能 会 有 些 不 同 。 从 图 3.28 ( 见 彩 图 ) 可 以 看 到 ， 
当 信 了 噪 比 为 30 : 1 时 ， 波 场 反 变换 曲线 的 两 个 峰值 吻合 得 都 不 错 。 但 是 由 图 
3. 26 和 图 3. 28 可 见 ， 当 信和 品 比 降低 到 20 : 1 时 ， 波 场 反 变换 曲线 出 现 较 多 波 
zj, 干扰 了 两 个 单 峰 的 识别 。 
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图 3.28 负 回 双 脉 冲 含 品 时 间 域 啊 应 及 变换 


由 图 3. 29 〈 见 彩 图 ) 和 图 3. 31 ( 见 彩 图 ) 可 见 ， 信 噪 比 为 30 : 1 的 时 间 
域 响应 曲线 已 不 太 光 滑 ， 由 图 3. 29 和 图 3. 31 可 见 当 信 品 比 降低 到 20 : 1 时 ， 
时 间 域 响应 曲线 “毛刺 ”明显 。 由 图 3. 30 ( 见 彩 图 ) 和 图 3.32 ( 见 彩 图 ) 的 
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波 场 反 变 换 图 上 可 以 看 到 ， 当 信 噪 比 不 低 于 30 : 1 时 ， 反 演 值 与 期 望 值 吻 合 得 
很 好 ， 一 正 一 负 两 个 波峰 定位 很 准确 ， 但 是 当 信 品 比 达 20 : 1 时， 波峰 定位 已 
不 是 很 准确 ， 而 且 波 场 反 变换 得 到 的 曲线 出 现 多 个 正 负 波峰 ， 造 成 了 假象 。 
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图 3.29 正 负 向 双 脉冲 的 含 噪 时 间 域 响应 
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图 3.30 正 负 向 双 脉 冲 含 噪 时 间 域 啊 应 反 变 换 


由 于 波 场 反 变 换 所 对 应 的 系数 矩阵 条 件数 很 大 ， 使 得 方程 组 的 求解 很 不 稳 
定 。 右 端 项 的 微小 扰动 就 可 能 严重 影响 求解 结果 。 由 噪声 实验 可 见 ， 当 信和 品 比 


. 71 ° 


瞬 变 电磁 法 拟 地 震 偏 移 成 像 研究 


不 低 于 30 : 1 时 ， 波 场 反 变换 曲线 与 波 场 期 望 曲线 可 以 吻合 得 很 好 。 
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图 3.31 负 正 向 双 脉 冲 的 含 品 时 间 域 啊 应 
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图 3.32 负 正 向 双 脉 冲 含 品 时 间 域 啊 应 反 变 换 


总 体 看 来 ， 波 场 反 变换 对 时 间 域 响应 的 测量 精度 要 求 是 很 高 的 ， 只 有 在 保 
证 测量 精度 的 前 提 下 进行 波 场 反 变换 才 有 意义 。 
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4.1 合成 孔径 成 像 算法 


瞬 变 电磁 波 场 变换 算法 的 提出 ， 实 现 了 由 具有 扩散 特征 的 瞬 变 场 癌 虚拟 波 
场 的 转变 ， 这 就 为 实现 航空 瞬 变 电磁 法 的 合成 孔径 成 像 创造 了 条 件 。 瞬 变 电 磁 
合成 孔径 成 像 技 术 是 借用 合成 孔径 雷达 成 像 的 思路 ， 利 用 机 载 具 实 孔径 发 射线 
央 与 目标 的 相对 运动 ， 把 尺寸 较 小 的 真实 天 线 孔径 用 数据 处 理 的 方法 合成 一 个 
较 大 的 等 效 孔 径 的 发 射线 圈 ， 使 其 分 辩 能 力 更 高 、 穿 透 能 力 更 蝇 。 对 于 航空 瞬 
变 电 磁 法 而 言 其 观测 方式 与 机 载 合成 孔径 雷达 十 分 相似 ， 完 全 可 以 借助 合成 筷 
径 雷 达 的 成 像 思 想 ， 实 现 虚 拟 波 场 条 件 下 的 瞬 变 电磁 合成 孔径 成 像 。 

瞬 变 电磁 场 经 过 波 场 变 换 ， 已 经 把 原来 的 感应 场 转换 成 了 波 场 ， 每 一 点 的 
数据 相当 于 变 成 了 自 激 目 收 的 波动 场 。 而 前 人 已 经 通过 实验 分 析 确 定 瞬 变 电 磁 
场 在 多 激励 源 情 况 下 存在 场 的 相关 去 加 性 ， 基 于 瞬 变 电磁 场 的 上 述 特点 ， 采 用 
相关 合 加 的 方法 来 进行 合成 孔径 。 其 合成 示意 图 如 图 4. 1 所 示 。 


2N+1 
Z=? 
















N// 
区 






图 4.1 合成 孔径 示意 图 


首先 选取 一 个 中 心 点 ， 取 为 第 i 点， 此 点 的 波 场 值 可 表示 为 U(r;，r)， 其 
Pr Ki 点 到 一 N，…，N 内 茶点 的 距离 ，r 为 相对 时 移 量 。 然 后 选 定 2N 十 1 
个 测 点 的 长 度 为 合成 孔径 的 长 度 ， 即 选取 i 点 左右 两 人 出 从 一 NN 到 NN 的 测 点 分 
别 导 中心 点 做 相关 ， 其 归 一 化 的 互相 关系 数 为 


è Ji s 
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>`UCr, sË; U(r sÉ; i r) 
=] 


Í P> [TO 61] > [Cr PY 
式 中 ，m 为 每 一 测 点 的 时 间 道 数 。 | 

互相 关系 数 olro r) 表示 了 两 列 波 场 的 相关 程度 ,但 它 与 相对 时 移 量 
rz 值 有 关 ， 通 过 改变 时 移 量 * 值 ， 找 出 最 大 相关 系数 o (r, r") 对 应 的 
时 移 量 值 ， 称 为 最 佳 延 时 ， 记 为 *。 由 此 可 以 得 到 2N+ 1 个 最 大 相关 系 
数 omx (ro C) 和 最 佳 延 时 区 (k= 二 一 N，…， 一 1，0，1，…，N)。 之 后 ， 
把 从 相关 计算 得 到 的 最 大 相关 系数 作为 权 系数 ， 用 权 系 数 分 别 乘 上 各 点 的 波 场 
值 ， 友 加 到 中 心 点 ， 最 终 可 得 到 中 心 点 的 合成 值 为 





p(riyr) = (4.1) 


N 
U Crati) — >, p (r cE Uir si — ri) (J mn 1 ,2 ,m) (4. 2) 
k=—N 


在 一 条 测 线 上 依次 移动 ， 则 可 得 到 第 itl, (2, i+3, PE Agr U sa 
的 合成 值 。 

以 上 是 为 了 说 明 原 理 ， 只 给 出 了 剂 面 上 的 一 维 合成 孔径 算法 ， 也 可 以 推广 
到 平面 上 的 二 维 合成 孔径 计算 ,进一步 还 可 以 发 展 日 聚焦 算法 等 。 


42 子 波 宽度 压缩 


虽然 变换 出 来 的 波 场 在 形式 上 与 地 震波 场 一 样 都 满足 波动 方程 ， 但 是 由 于 
两 种 波 场 的 物理 背景 不 同 ， 它 们 之 间 存 在 一 些 重 要 区 别 ， 前 者 是 与 感应 的 瞬 变 
电磁 衰减 曲线 相对 应 的 “反射 ” 子 波 ， 是 虚拟 的 ， 而 后 者 是 在 弹性 介质 中 传播 
的 客观 存在 的 地 震 子 波 . 另外 ， 虚 拟 波 场 在 每 一 种 介质 中 的 传播 速度 不 仅 与 本 
介质 的 电导 率 有 关 ， 而 且 还 受 相 邻 介 质 电导 率 的 影响 ， 这 与 地 震波 场 是 不 同 
的 ， 并 且 存 在 严重 的 波形 展 宽 效 应 ， 使 得 计算 得 到 的 虚拟 波 场 数值 分 辩 率 降 
低 ， 成 为 瞬 变 电磁 法 拟 地 震 成 像 解 释 的 主要 障碍 之 一 。 显 然 ， 这 一 现象 严重 影 
响 着 TEM 成 像 的 空间 分 辨 能 力 。 所 以 ， 有 必要 讨论 这 种 虚拟 的 波动 场 波 形 展 
宽 的 原因 及 对 策 。 

以 下 首先 给 出 了 波形 展 宽 的 现象 ， 然 后 讨论 了 虚拟 波形 展 宽 的 原因 ， 认 为 
不 是 由 于 波 在 介质 中 传播 时 能 量 损耗 所 致 ， 而 是 波 场 变 换 式 中 高 斯 分 布 的 核 函 
数 随 虚拟 时 间 的 增加 分 布 范围 增 大 的 分 布 特点 所 决定 的 。 最 后 对 虚拟 波 场 的 离 
散 数据 求 取 反 裙 积 ， 消 除了 波 场 变换 的 波形 展 宽 效 应 。 通 过 对 理论 模型 计算 处 
理 ， 证 明了 该 方法 可 以 压缩 子 波 的 宽度 ， 并 增强 瞬 变 电磁 法 识别 地 下 电 性 分 界 
面 的 能 力 。 
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4.2.1 波形 展 宽 的 现象 


以 三 层 地 电 断 面 为 例 进行 波 场 变换 ，K 型 断面 地 电 参 数 为 , pi = 二 5Q，m， 
hi=60m, o;=50Q * m, hs 二 60m， ps 二 5Q，m。A 型 断面 地 电 参 数 为 ，o = 
5Q * m, hi =60m, o;=50Q * m, h;=60m, o =500Q ° m, Q 型 断面 地 电 参 
数 为 ， 0 ¿250 * m, hi =60m, o =50O * m, h; =60m, o = 10 = m, 图 4.2 
为 以 上 三 个 模型 的 波 场 变换 结果 ， 对 三 个 模型 的 变换 结果 进行 分 析 ， 其 共同 的 
特点 就 是 波形 较 宽 ， 正 负 脉 冲 之 间 的 时 间 间 隔 较 大 ， 导 致 该 方法 在 拟 地 震 成 像 
处 理 中 对 电 性 界面 的 分 辩 能 力 下降 。 





Š 
Ë 
RE 
时 间 /sz 
(a) 
8 
4 
b 0 
Ë -4 
g 
-g| 





时 间 /s 


(b) 
图 4.2 地 电 模 型 的 波 场 变换 结果 
(a) K 型 断面 的 波 场 变换 结果 ; (b) Q 型 断面 的 波 场 变 换 结果 ; 
(c) A 型 断面 的 波 场 变换 结果 
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图 4.2 地 电 模 型 的 波 场 变换 结果 ( 续 ) 


4. 2.2 ”波形 展 宽 的 原因 


在 瞬 变 电磁 场 半 化 为 虚拟 波 场 的 过 程 中 ， 所 得 虚拟 波 场 波形 会 随 着 时 间 的 
增 大 而 不 可 避免 地 展 宽 ， 严 重 影 啊 成 像 的 分 辨 率 。 原 因 在 于 波 场 变换 式 (1.3) 
EHA% pR 


nt 





a(t,r) = (4. 3) 


的 特点 ， 即 高 斯 分 布 的 核 国 数 以 真实 的 物理 时 间 为 方差 ， 因 此 随 着 时 间 的 增 
加 ， 时 域 场 结束 的 时 间 越 滞后 ， 波 场 波 形 也 就 越 宽阔 。 尽 管 这 一 现象 是 由 数学 
变换 产生 的 ， 但 是 与 地 震 勘 探 中 的 “大 地 滤波 作用 ”非常 相似 。 从 滤波 角度 
看 ， 要 提高 分 辩 率 抵消 大 地 滤波 器 的 作用 ， 求 取 一 个 反 滤波 因子 ， 令 其 能 抵消 
大 地 滤波 的 作用 ， 则 分 辨 率 即 可 大 幅 提 高 。 而 这 可 以 通过 反 锐 积 运算 来 完成 。 
它 经 过 理论 和 实践 的 证 明 完 全 可 以 把 先前 得 到 展 宽 的 波 场 能 量 集中 化 ， 达 到 精 
细 分 层 的 效果 。 为 了 说 明 核 肾 数 的 分 布 宽度 随时 间 增 大 而 增 大 ， 按 照 式 
(2.22), 分别 计 算 t=80us, 325us, 800us, 2. lms、8. 7ms、27ms、8lms 时 
核果 数 值 随 虚拟 时 间 r 变化 时 的 分 布 形态 。 

由 图 4. 3 可 见 ，: 二 80ps 时 ， 核 肾 数 对 应 的 z+ 值 的 最 大 取 值 范围 为 0~ 
0.05s? ; t=325ps 时 ， 核 函数 对 应 的 r 值 的 最 大 取 值 范围 为 0 一 0. 1s? t= 
800us 时 ， 核 图 数 对 应 的 r 值 的 最 大 取 值 范围 为 0 一 0. 15sz ; t=2. lms Ff, 3 
函数 对 应 的 z 值 的 最 大 取 值 范围 为 0 一 0. 25sš ; t=8. 7ms Rf, RROTA r 
值 的 最 大 取 值 范围 为 0~0. 5st; 1=27ms 时 ， 核 函数 对 应 的 * 值 的 最 大 取 值 范 
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围 为 0 一 0.85sz; t 一 8lms 时 ， 核 函数 对 应 的 + 值 的 最 大 取 值 范围 为 0 一 
l. 5s, 随 着 时 间 t WIER, IAR a (t, +) 随 虚 拟 时 间 r 的 增加 分 布 范围 明显 
增 大 ， 这 就 导致 经 反 变换 获得 的 虚拟 波 场 出 现 了 随 着 时 间 的 推移 而 展 宽 的 现象 。 
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(a) (b) 
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; 0.4 0.8 1.2 
虚拟 时 间 /s 虚拟 时 间 /; 
(e) (D) 
图 4.3 核 函数 值 分布 性 态 
(a) t= 二 80ps 时 核 函数 值 分 布 ，(b) :=325us 时 核 函 数值 分 布 ; (ec) t 二 800ps 时 核 函 数值 分 布 ; 
(d) t= 二 2. 1ms 时 核 函 数值 分 布 ;(e) 上 一 8. 7ms 时 核 函 数值 分 布 ;(fT) t=27ms 时 核 函 数值 分 布 ; 
(g) =81ms 时 核 函 数值 分 布 
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虚拟 时 间 /st 
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图 4.3 核 函 数值 分 布 性 态 (#k) 


表 4.1 列 出 了 不 同时 间 的 核 消 数 分 布 情况 的 统计 ， 由 表 4. 1 可 以 看 出 ， 时 
加 越 大 ， 核 曙 数 分 布 的 范围 越 大 ， 说 明 随 时 间 的 推移 ， 由 于 核 男 数 分 布 的 范围 
加 大 ， 导 致 虚拟 波 的 宽度 加 大 。 


表 4.1 核 函 数值 分 布 性 态 (范围 ) 统计 


序号 时 间 / 虚拟 时 间 r 范 围 /sz ”虚拟 时 间 宽 度 /s? 
l SOps 0—0. 05 0. 05 
2 325ps 0 一 0. 1 0.1 
3 800ps 0~0. 15 0. 15 
4 2. Ims 0—0. 25 0. 25 
9 8. 7ms 0 一 (0. 3 0. 5 
6 2 Tms 0 一 0. 85 0. 85 
81ms 0—1.5 L5 


另外 ， 在 自由 空间 中 ， 电 磁 波 的 波 速 为 co 二 1/ Vjyweo， 而 导电 媒质 中 的 虚 
拟 波 速 为 v=1/ Vuo., MERZA 6 个 数量 级 。 由 于 导电 媒质 中 的 虚拟 波 速 非 
常 慢 ， 这 样 ， 传 播 同样 的 距离 差 ， 就 需要 更 多 的 时 间 。 所 以 ， 其 波形 的 展 宽 效 
应 就 变 得 很 显著 了 。 
4.2.3 利用 反 裙 积 技 术 压 缩 波 场 宽度 
波形 展 宽 这 一 现象 与 地 震 勘 探 中 的 “大 地 滤波 作用 ”非常 相似 。 在 地 震 勘 
探 中 ， 由 于 大 地 滤波 的 作用 ， 加 上 噪声 和 多 次 波 的 到 加 ， 使 本 来 可 以 清晰 呈现 
出 的 地 层 层 序 尖 脉冲 序列 “模糊 化 ”， 降 低 了 地 层 的 纵向 分 辨 率 。 反 裤 积 的 作 
用 就 是 压缩 子 波 ， 可 以 通过 求 取 一 个 反 滤 波 因 子 ， 令 其 能 抵消 大 地 滤波 的 作 
用 ， 可 以 把 先前 得 到 展 宽 的 波 场 能 量 集 中 化 ， 达 到 精细 分 层 的 效果 ， 这 样 分 辩 
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率 可 望 大 幅 提 高 。 

Peacock 和 Treitel 在 1969 年 提出 预测 反 窜 积 ， 此 项 技术 在 地 震 资 料 处 理 
中 得 到 了 非常 广泛 地 应 用 。 这 项 技术 同样 可 以 用 来 处 理 转换 出 来 的 虚拟 地 震 
子 波 。 

假设 y G) 为 实际 转换 出 来 的 拟 地 震 子 波 ，y GO = (y (O), y (1)，…， 
y(nm))， 我 们 希望 通过 一 个 滤波 因子 h(t) = (h (Cm), h (—m +1), 
h (>m +2), *, h (—moə+m)), h (G) 的 起 始 时 刻 为 一 ro， 其 延续 长 度 
为 (m+1), ， 进 行 如 下 计算 


yi(t+a) = ON (n—t) (4. 4) 


然后 ， z 后 的 子 流 》(t 十 a)， 并 且 使 该 时 刻 拟 地 震 子 波 与 
计算 误差 值 Q 最小， 


Ly(t +a) — y(t + aa) J (4.5) 
把 式 (4.5) RAA (4.6), 4 
Q= SLs +a) — DRO (pay (4. 6) 


要 使 Q 为 最 小 ， 数学 上 就 是 求 Q 的 极 值 问题 将 Q 对 滤波 因子 求 偏 导 ， 
化 简 ， 并 令 其 为 零 ， 即 


aQ _ 
Bi en 


Futa 

00 _ ` = is 

T T aa 之 znl? z: a(r)b(t — r) de) | 
-my Hrtan —m; +m 


=2 > | D abatre) — dt) bai) 


E= r= 


=2 Š h) >) yG—r)yG—D—2 > yt)y(t—l)=0 (4.8) 
t=—mç r= =m 


({=—mo, “mtl, °, —m tm) 


mg Hrn 
RP, 2 yG — oyt- D = rw(l — r) 是 虚拟 子 波 的 自 相关 项 ， 而 
t= 


-my tmtn 
5) Oyat- = ra (Q) 是 子 波 与 期 望 输出 的 互相 关 项 ， 所 以 方程 (4. 8) 
一 ID 
可 写成 
mg +m 
> h(t)re (Í — r) = ra(t) (Í =— m, — mM, + l, — m +m) 
r= ht 


(4.9) 
L 80 Š 
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将 式 (4.9) 成 矩阵 形式 为 


Teb (0) rpp CL) ... Yb (0) | h(— mo) rpd (— Mh ) 

rep CL) rs (0) à rpm — i) | hec m + 1) ra — m + 1) 

bm) rlmO— 1) ... rep (0) (一 m +m) Lra — m + m) 
(4. 10) 


HARDE (4.9), KH Ef h GO 后 ， 就 可 以 通过 式 (4.4), 
求 出 转换 虚拟 波 所 对 应 的 反 禄 积 处 理 计算 结果 。 

反 滤波 算 子 的 长 度 原则 上 可 以 任意 选择 。 地 震中 反 滤波 算 子 的 长 度 一 般 选 
择 80ms、120ms、160ms、200ms、240ms 等 ， 但 我 们 在 计算 时 要 将 它们 换算 
为 采样 点 数 。 要 根据 插值 后 数据 的 总 个 数 来 进行 选择 ， 当 达到 一 定 长 度 时 ， 反 
军 积 效果 趋 于 稳定 。 在 本 文 的 计算 中 ， 反 滤波 算 子 的 长 度 基本 选取 在 400 到 
600 之 间 。 


4.2.4 模型 计算 


图 4. 4 为 运用 了 反 宰 积 技术 的 结果 ， 图 4.4 (a) 是 反 裙 积 处 理 前 的 波形 ， 
在 反射 界面 处 虚拟 反射 波 波形 明显 较 宽 ， 分辨 力 下 降 ; 图 4.4 (b) 反 裙 积 处 
理 后 波形 ， 出 现 两 组 虚拟 反射 波 ， 波 形 相 对 较 罕 ， 由 图 可 见 ， 通 过 反 初 积 处 
理 ， 可 以 将 电 性 分 界面 处 产生 的 被 拉 长 的 虚拟 反射 波 锐 化 ， 使 瞬 变 电磁 法 的 分 
A EIA BERRIE. 





图 4.4 反 袜 积 计算 前 后 的 H 型 模型 子 波 对 比 
(a) 压缩 前 子 波 波 形 ;，(b) 压缩 后 子 波 波形 
H 型 模型 参数 ; m =p =1000 * m，m 一 50 * m, hi =hz=100m 
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为 了 验证 反 裙 积 的 处 理 效 果 ， 将 对 图 4. 4 中 地 电 模 型 的 反 锐 积 处 理 后 的 虚 
拟 波 场 进行 拟 地 震 成 像 ， 图 4. 5 WEED 给 出 了 H 型 模型 的 成 像 的 结果 ， 理 
论 上 ， 地 电 模 型 的 两 个 电 性 界面 分 别 位 于 地 下 100m 和 200m 处 。 由 图 4. 5 可 
知 ， 浅 层 电 性 界面 所 产生 的 波 场 能 量 集中 于 100m 左右 ， 而 较 深 界 面 能 量 高 峰 
产生 于 190m 到 200m 之 间 。 说 明 在 界面 深度 定位 方面 与 模型 预先 给 定 情况 大 
体 一 致 ， 效 果 已 令 人 满意 。 





图 4.5 H 型 模型 拟 地 震 成 像 图 
(a) H 模型 成 像 图 ; (b) H 模型 成 像 切 片 图 


距离 /m 
40 60 





(a) (b) 
图 4.6 MRR a — #Eba 8-1 XJ IE. 
(a) 三 维 地 电 模 型 ，(b) 计算 场 点 平面 分 布 ; (c) 压缩 前 子 波 波形 ; 
(d) 压缩 后 子 波 波形 ; (e) 拟 地 震 成 像 剖面 
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为 了 验证 反 裙 积 的 处 理 效 果 ， 我 们 设计 了 三 维 地 电 模 型 ， 取 均匀 半空 间 的 
电阻 率 为 pi 二 10Q0，m。 异 常 体 为 30mX 30m X 50m Wiki, 电阻 率 为 p = 
300Q，m， 埋 深 为 70m。 

从 图 4.6( 见 彩 图 ) 中 可 以 看 到 ， 经 反 滤 波 处 理 后 ， 波 形 得 到 锐 化 ， 上 部 
波形 能 量 集 中 在 70m 处 ， 下 部 波形 能 量 集中 在 120m 处 ， 与 模型 设置 吻合 得 很 
好 。 特 别 是 下 部 从 高 阻 到 低 阻 的 界面 ， 本 身 是 弱 反 射 面 ， 波 形 展 宽 很 严重 ， 经 
REE. HEE T KAE. 
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前 面 我 们 已 经 解决 了 瞬 变 电磁 场 的 波 场 变换 问题 ， 并 通过 实验 研究 证 明了 
多 孔径 瞬 变 电磁 场 确实 存在 相关 全 加 的 波 场 性 ， 这 就 为 瞬 变 电磁 场 的 偏 移 成 像 
创造 了 条 件 。 但 瞬 变 电磁 场 偏 移 成 像 问题 与 地 震 勘 探 中 的 弹性 波 偏 移 成 像 问题 
有 相似 之 处 但 也 不 完全 相同 ， 如 果 将 扩散 的 瞬 变 电磁 场 变 换 为 波 场 处 理 ， 将 波 
场 分 析 的 原理 ， 即 地 震波 场 从 地 面向 地 下 反 向 外 推进 行 偏 移 成 像 的 方法 用 于 对 
瞬 变 电磁 场 的 解释 ， 形 成 瞬 变 电磁 偏 移 方 法 。 也 就 是 用 基 尔 霍 夫 积 分 法 进行 电 
磁 波 场 偏 移 成 像 处理 ， 实 现 电磁 波 场 的 延 拓 计 算 ， 这 是 目前 瞬 变 电磁 测 深 方法 
研究 的 焦点 和 前 注 诛 题 。 本 章 对 该 方法 进行 了 探索 性 的 研究 ， 分 为 三 节 内 容 。 


5. 1 基 尔 替 夫 绕 射 积分 的 建立 


如 图 5. 1 所 示 ， 取 一 区 域 WW， 其 边界 为 Q, n 为 Q 的 外 法 线 ， 设 在 上 时 
Z|, M (zi, ys 2) 点 上 观测 到 的 势 郴 数 为 (tis ys zs t), MÆR — 
时 刻 所 一 上 一 二 ， 在 闭合 曲面 Q 上 的 位 移 势 是 p ays z, t——) = [gç], 
称 [o] 为 “推迟 势 ”， 而 > 为 M 点 至 Q 各 点 的 距离 








和 
图 5.1 区 域 与 边界 


为 了 计算 Q 上 各 位 移 势 Lp] 的 总 和 要 引用 格林 公式 
[u Vv—vVudw = (ue -v2 )aQ (5. 1) 
w Q 
AF, u 和 w HERRA, HX 
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s = @(z,y,z,t— T) = [ç J 
因为 需要 由 zt 的 前 一 时 刻 二 :一 二 的 推迟 势 [p] 来 确定 的 解 。 取 "一 
二 ,由 于 g 一 二 C1 (T), 二 是 gp 的 基本 解 ， 因 此 子 波 的 球面 扩散 取决 于 二 ， 
考虑 式 (5.1) 的 左 端 
u V‘u = Le] v>) 
由 于 


式 中 ， 二 是 拉 普 拉 斯 方程 的 解 ， 故 w Vv 一 0， 所 以 式 5. 1) EX 
-|4 view d an Pr jQ 


r d or . 
考虑 到 -一 一 全 。 Ar RS 
E- 1 Vz[o]dw = M. d EE S ee alel ag (5.2) 


dn r an dr 
ERZA REMV a= L Vš e] 而 
[Ve] = [V* Vo] = [V E] (5.3) 
式 中 ， 
E = Vo, p= 二 CiG 一 工 ) 
所 以 ， 有 
[Vo] = [LO 
At, [0] 表示 在 6 一 上 一 二 时 刻 的 激发 源 的 体积 相对 变化 ， 而 且 根 据 激发 源 
强度 函数 的 定义 ， 式 (5.2) 左边 的 体积 分 表示 的 就 是 激发 源 的 强度 。 格 林 公 
式 右边 曲面 积分 内 的 两 项 ， 是 由 于 激发 源 作 用 后 在 Q@ 面 上 的 解 . 不 过 此 时 该 QQ 
面 上 的 解 作 为 二 次 子 波 的 新 激发 源 ， 因 此 对 Q 面 作 面 积分 实际 上 是 求 这 些 作 
为 新 激发 源 的 二 次 子 波 在 Q 面 上 的 总 和 ， 于 是 式 (5.3) 所 描述 的 意义 是 ; 被 
Q 曲面 包围 的 W 区 域内 ， 在 时 刻 yman 一 次 震源 所 产生 的 能 量 应 该 等 于 
此 时 此 刻 Q 面 上 产生 的 二 次 子 波 的 总 和 。 因 此 求 某 一 点 上 由 激发 源 引 起 的 位 
BAHR, A (5.2) 的 意义 是 关键 。 下 面 建立 用 lo] 及 其 导数 表示 的 积分 
ya 
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dd I SE: kisa =- =. 
= d = idii v Dr v Jz It A U 










































2u d. arrap | 
[2] 一 二 了 [了 2 (5.4) 
同样 有 
9 9 9 | 
m = PS 1= E ErP Nan = Ə t) (5:9) 
2. kV? 
由 于 
w, = 9 4 二 gu 9u 
Ə x° 3 y’ J >° 
由 式 (5 
az _ 3 391 1araragj larap 
azz ` Pr U To ia U myra (S. 6) 
由 于 
a me rd _ l ear 2 IP wl dr dr dg 
rp Fa U k Fa bs š u 9x | 
9 [99] (9 — 12 OT Panl BF. Fan 
EFE or U Ds EL “= r: u dr Fš 
RAR (5. 6) 得 
S u 了 A tr e ) rag 1 azrrag] 
同 理 得 
9 šu 3’ 2 dr 1 3fr-3ğ 
3 y’ = 152 一 v EEE FEF + sp af- v ay It- 
Ə “wu k _— Ë dř 1 Ə ° 4 -ia drr 29% 
F 52 F- s pE 4 y [a = "i r Eya 
于 是 有 
2 3 
2 oe 
w LV g. 2 aW IAT 
1r r3? yr ar 1lragN 9°r 
Yara LT 591 s W: (5. 7) 


$h, W,=z, W,=y, W =r. WF = (x—=z, + ia Hie) 
我 们 得 到 
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dr r d y r dz r 
所 以 

As Ar r— It" — y. z— z 
e S s 十 (一气 ) 十 (一) =] 
i=] : r r 
J 

ary 3 ara 2 r, 2 ¿3r,— 2 
HGW o daos aa a r 


此 外 ， 考虑 到 波动 方程 
1 


2 
Ve = F LR +LFG,y zD 


将 以 上 各 式 代入 方程 (4.7) 合并 同类 项 ,得 
;2ralg) 2ragp 2N arr ai 
vam Ls P] == =i 2, w. e r (5. 8) 


1 Ë 





于 是 格林 公式 左 端 被 积 函数 一 Vzw 可 写成 

lyu = 22e 2r3a8]_ 2y W, —Wyi- ap J, F 

| wiar JŽ r“ (wot 
(5.9) 


为 了 把 上 式 写 成 一 个 有 意义 的 形式 ， 考 虑 如 下 运算 


a [znr ap) 9 (xz—a) .rap  z—m 3 rag 
arl r° [51] Ja 产 ) Arpa s; = Pipra 


_ lrəƏg,_2(z=—m) rap， 工 一 ZIr 日 _ (z—z)'- 32? 
Salad a La It laza a lazy 





aamke 222 [9]| 和 地 -和 各 (52)), ARAA 





š 可 W. — W, d _ ] dg 1 J? - W: — Wy 2° 
2 | z 5E) 2[ 1 s. lO, r* ELAT 


式 中 ,Wi 一 x，W, 一 y，Ws 一 z。 上 式 乘 以 一 之 再 加 上 全 就 等 于 式 (5.9)， 


于 是 
l vy, =— 2 2 f. 3h 393f 1 __ 2 v. f+ 上 (5. 10) 
r g x 可 y dz r V r 


若 邻 
| 1 d == n 1 ag = # 21 2 g == 
Tean ea Lare E 


A (5.10) 可 写 为 
l wz， —— 2 ¿2 f: 2b L 2 y E =— Z v. f+ 二 
r v gz d y dz r U F. 
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由 散 度 定理 
| v- faw = 中 da 
W Q 
可 以 得 到 
2 — F __2 I 0 他 人 — o“) 
| > yiudW = = V. faw + | 7 F aw aiai a 


由 于 f,=— E], BD 
[tv aw- É =— 21 [SF]dQ (5. 12) 


Rh, n 为 激发 源 到 各 参考 点 的 距离 。 — 
将 式 (5.12) MR 6.5) 代入 格林 公式 ， 并 考虑 到 w= [g]， 当 M 点 位 
于 Q 外 面 时 ， 格 林 公式 左 端 为 
viv vv dw =— J: Vk = Zra fJ Q+ Ë 
格林 公式 右 端 为 I 
ETET ain p| Ee sr 了 (一 ) 一 一 [eJ |dQ 
一 中 [由 (一 点 ) 9 一 了 [[352] 一 S52]])dQ 
根据 格林 公式 左 端 等 于 右 端 ， 有 
2p 1 Sr[2%Jaq+ Z = -去 CR UZAN 


vj v dn-at ro 
Q 


ERERRNTAN 
Ia ananas em 


RH, F=uð (1 一 0* )。 此 时 ， TI M 点 的 9 EAE, 
3. M 点 位 于 W 内 (W 为 扰动 区 域 ) 
当 产 =0 时 ，v 一 一 一 co， 格 林 公 式 不 成 立 ， 为 此 必须 去 掉 奇 点 M， 在 M 
点 周围 用 一 个 半径 为 e 的 小 球 包围 M 点 构成 5 域 ， 在 Q 面 和 小 球面 c 之 间 的 
区 域 应 用 格林 公式 ， 并 利用 式 (5. 13) 的 结果 


_lrap_1 9rr9pgp = 一 4) 
$ lz, r i piia Ur Ən gs i wa ay 
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考虑 
[aAa 3s3n]a 


= (方向 相反 )，do 二 edQ9，d0 为 立体 角 。 于 


由 于 e 很 小 ， 故 好 -= 


是 ， 上 式 变 为 
_r_ e ls iri fy l 2eyr 9 
| ra + + Cd 


a —— 


do = #ag|||[ g] +L 32]+ [32] dn 


局 











£ 一 > o| Colan = 4r| ç] 





代入 式 (5.14) 中 有 
| ady llr2Ə@. _ 1 2rr 2° = £. 
Jue: TA pm u [2 J) aQ + ng ro 
所 以 区 域 W AM 点 的 位 移 势 
| CESEN | Fels AF908 1 arrg E 
plam» ym» zmt) 一 一 | (ia sC — L[3$]— 1 2r 2e] j+ 
(5. 15) 


由 于 式 中 F = uó (t— o" ), 只 有 在 t= 0 处 ， F (0) 一 Moi 1&0 时 ， 
F (1)= 二 0。 在 一 般 情况 下 t 了 关 0， 故 不 需要 考虑 此 项 。 


5.2 基 尔 霍 夫 积 分 偏 移 (曲面 延 拓 ) 


我 们 知道 ， 波 场 在 地 下 传播 可 以 用 波动 方程 描述 ， 其 纵波 波动 方程 为 


k. É. diu __ r 
V‘ u gr F 


上 述 方程 的 基 尔 霍 夫 积分 解 为 
一 


人 一 上 让 L. 
u(r,y,z,t) = 4 an r rodan w dn dt 





(5. 16) 
式 中 ， 上 一 po6 (t—o*), Q=Q., +Q, 为 闭合 曲面 ， 如 图 5. 2 所 示 ，Q@Q 为 地 面 ， 
Q, 为 无 限 大 半球 面 ， 利 用 +r 一 co 时 ，u->0， si, 所 以 Q 的 曲面 积分 对 M 


(X, Y, Z) 点 波 场 函数 的 贡献 为 零 ， 故 
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dQ = 0 





[Zh +[5]— 1 SEL] 
Q 

因此 ， 式 (5.16) 可 写成 

ult ys zt) — 圳 (3 es se 二 [一 


Ən r 





图 5.2 区 域 与 边界 


虚拟 波 场 偏 移 处 理 是 获取 记录 的 逆 过 程 ， 现 在 我 们 已 知 的 是 地 面 上 (Q. 
E) 观测 点 的 电磁 响应 的 虚拟 波 场 记录 ， 需 要 确定 反射 面 上 作为 二 次 虚拟 辐射 
源 的 空间 位 置 。 对 于 xx (z, y, z, t) 而 言 ， 当 把 t 改变 成 一 :时 , 即 W (x<, 
ys z, t) =u (x, y, z, —t BF, W (z, x, z, t) 仍 可 满足 同样 的 波动 方 
程 。 对 于 W (xz, y, z, t) EWER ERTA, Fu (z, y, z, t) 就 
是 时 间 “ 人 倒退” 的 问题 。 可 以 把 反射 界面 的 各 点 看 为 同时 激发 上 行 波 的 源 点 。 
这 样 ， 我 们 可 以 把 地 面 上 的 接收 点 作为 二 次 发 射 源 ， 将 这 些 信 息 值 时 间 “ 倒 
退 ” 到 原来 状态 ， 寻 找 反 射 界面 的 波 场 函数 ， 以 确定 反射 界面 。 

设 自 激 自 收 的 上 行 波 为 G(x，y，zo，)， 它 是 地 下 反射 界面 作为 源 点 发 
射 的 波 场 g (w, y, z, t) 在 地 面 上 z=z 上 的 值 ， 由 式 (5.17) 可 得 


OSE.. 2 ily l 2 _1 LEAT T Lr 
g(TrYys Tt) = 可 | Tan = i =e 3: G(&,m b +Q 
pE (5. 18) 
ro 


RP, G (&, p h tHo Wito, EH TE BHEE HNE, 


事实 上 是 由 g (xw, y, z, t) 在 地 表面 上 的 值 G (ë, p G, t) 求 出 地 下 的 波 
场 值 ， 从 而 确定 反射 界面 ， 这 即 是 波 场 的 向 下 延 拓 。 
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5. 3” 波 场 延 拓 的 边界 单元 法 


5.3.1 基 尔 霍 夫 积 分 的 离散 化 


H n eh AX Q, 进行 前 分 ， 将 式 〈5. 18) 中 的 边界 积分 分 解 为 各 单 
元 积分 的 积分 之 和 ， 当 了 点 位 于 W q 


1 _ 1 3 工人 r r 

saa T ai C-L 212r EOGI 5 十 工 )dQ 
二 (5. 19) 
ro 


5.3.2 单元 分 析 


如 图 5. 3 所 示 ， 设 三 角形 行 单元 顶点 编号 为 ;)，k 上 ，m; 各 节点 坐标 分 别 为 
(Tj Yj Z;), CErs Yrs Zk) 和 (zm，yn，zm)， 和 有 先 我 们 做 如 下 坐标 变换 ， 
HERA E., &, & 表示 单元 中 任意 点 的 坐标 。 


Cy) 





(Ki) 


图 5.3 三 角 单 元 
T= Xj; T Tb Tm Lj Th Y Ty, [总 | T 
y = y;Ë; + Yik 十 Ymêm™> > 9 YV >| El H (5 20) 
= 一 -AF + zaê, - Esa Z; * Zk * Zm m- 


AP, £, 6, & 是 取 值 范围 为 0 一 1 的 函数 ， 且 
Ë + Š + ë, = Ls 
在 了 Ag = 1,& = &, = 0 
在 有 点 ,6 = 1,&, = š, -| 
TE m Pu 6; = 1,& = &, = 0 
可 见 ， Gr Bs bn JE: x, Ys z 的 线性 图 数 。 
由 于 单元 Q, 一 般 取 得 很 小 ， 故 可 假定 波 场 G 在 各 单元 上 是 线形 变化 的 ， 
即 单元 上 的 G 可 表示 为 


(5. 21) 
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G, 


G = £G; +&G, +$,G,, = [&,,& ,6,, IG (5. 22) 








la, 
AP, G;, Grs Gn RB H Q AER GH. 
考虑 单元 积分 


32 , 1 r . 
|z; ç GCE, 1» bo t +7) dl 
=j- x Felep „t+ 一 <ur 


_ || cos(r ° n) PASA 
= I — G (ë, 1)» Ço s. sz a = dr 


i 


KA, Arcos (r` n), 将 单元 上 的 波 场 G (ë, 1» ns pepeg 统一 用 G, 表 
R, WERS | 
_ || sos(r * n) Gien z 
| 了 GCE, £ „t+ )d]' 


r. 


ea . cos(r + n) jp. Ë =— y AG (5. 23) 
e= jk; 


e= jkm 


式 中 ， 
= fje, Rar (= jka 
该 式 可 用 高 斯 求 积 公式 计算 
r= le. colr: D ape DEl 1 ady, “A 


AP, W, 为 加 权 系数 ， A = 角形 面积。 








同 理 ， 单 元 积分 
| +26ar= > a, S: (5. 24) 
sr dn ej ksm an 
AP, 
d, = le. + tar = DEW =W, (5. 25) 
Jx 
单元 积分 
l 3r 2Gm_ V; 3G. > 
Er = i r (5. 26) 
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其 中 ， 
ca = ||, SE Teo Leay, „d (5. 27) 


R u = 


5.3.3 BEBE 


将 式 (5.23). = (5.24). Ñ (5.26) 代入 式 (5.19) 中 , i 点 处 的 波 
场 为 


m = 一 下 | FG- D, s =A 2a] 
H nn 个 市 点 得 到 一 个 方程 组 
(5. 28) 


式 中 ， 

h = [F; J; D= [ D; J; C= [Gy]. 其 中 ， F, j 节点 上 相 邻 单元 fj 之 
和 ; D; g j 节点 上 相 邻 单元 dy 之 和 ; C, 为 ) 节点 上 相 邻 单元 cy 之 和 。 若 已 知 
地 表 上 的 波 场 值 G — (g g "S gjr， 波 场 法 向 导数 32 一 

dg dg Ü. | nd [dm Š> sl", 

(F aa e Da) PEARES I Dr) ar) RT 
以 求 得 地 下 某 点 处 的 波 场 值 。 

5.3.4 单元 积分 的 计算 


1. cos (r。n) 的 求法 
因为 cos (r * n) 一 士 D/r， 其 中 D 是 表示 所 考虑 的 节点 i 到 某 个 面 元 素 
所 在 的 平面 的 距离 ，r 是 节点 i 到 面 元 素 上 任 一 点 的 距离 。 
1) D 的 求法 
设 面 元 如 所 在 平面 为 
arx tby tetli =o (5. 29) 

然后 把 面 元 素 上 三 个 节点 的 已 知 坐 标 代 入 上 式 

faz, Tby;À Tz; + 1 = 0 

[tt est = 0 (5. 30) 

ar, + by, + z,, + 1 = 0 


Tj Yj Zj a | 
Ex Zh | ， 一 一 H (5. 31) 
La yi Esa c 1 
XEMA RH as b, c 
. 04 。 
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设 i 扩 到 该 平面 上 任 一 点 (zx，y，z) 的 距离 为 rp， 则 


rp ° = (x — x;)° + (y — yD + (z — z)? (5. 32) 
由 平面 方程 得 
„= aty 5.39 


将 上 式 代 入 rp 中 ， 并 分 别 对 z，?y 求 导数 ， 再 令 其 为 零 ， 得 
Br ar) tAE |) 2 = o 
| (5. 34) 
Fm’ _ N E a a N, bo 
d y aii C "E 

KE x, y 分别 记 为 zp，yp， 然 后 把 他 们 分 别 代 入 平面 方程 (5. 29) rh, 

可 求 得 zp， 则 
D = i(z=p — xi) +j Oyo — yi) + k(zxp — zi) (5. 35) 





ID| = (zp — xi)? + (yp — y) + (zp — z)? (5. 36) 


cos(r ° n) = (2; (5. 37) 


(z>) +y) Hlaz), /£ (z, y, z) 是 元 素面 上 的 





2) £p» yp» zp 的 推导 
HK (5.34) 可 知 


(z—z HOYT S =o (5. 38) 
(ay) ET ae =o (5. 39) 


H< (5. 38) 一 人 (5. 39) 得 


一 一 一世 = 0 
所 以 
了 (5. 40) 
将 式 (5. 40) RAR (5.37) 中 ， 得 
a Ley EEL] y 


(ey) thei r) a isi. 


C 


所 以 ， 有 
e ÜG š 
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— (c° + a2)y; — bax; — bcz; — b 


a2 + b° + C (5. 41) 
将 式 (5. 41) RAJ (5.40) 中 ， 得 
NP S ( +a" dy — mz; — sg — b _ à 
r= +Í a+ y); he: | 
Fa = DEEN 
— W +h )ri — aby; — acz; —a e -a a a (5. 42) 
将 式 (5.41) 和 式 (5.42) RAR (5.33) 中 得 
i Tl = (a° + b )z=; —acx; — bey; — c (5. 43) 





c a? + b° + c° 
将 式 (5.41). 式 (5.42) 和 式 (5. 43) 中 的 Ty Ys z 分 别 用 Tps YD’ 
ZD 表示 ， 得 到 
kas; (b° + c° )=z; — ab F 
a° + b° +c’ 
c° — ! ESE 
ss == 和 (5.44) 
Xa to Je —acz; —Bcy; — £ 
a? + b° + c° 





2. 单元 积分 的 计算 

以 单元 积分 方 的 表达 式 为 例 说 明 其 计算 方法 ， 由 上 节 推 导 可 知 
=] sar . n) p — Ds ) Sat n) "E 

取 4 点 的 高 斯 积分 公式 ， 由 高 斯 积分 公式 表 5.1 知 


f; = cos(n en), gipa 3 cos(rs* M) (25/48)A 
2 


Ti coarsi* 到 (25148)A 十 革 cos(r, * n) (25 /48)A 
rs ra 
表 5.1 四 点 高 斯 积分 系数 表 
g g” &? g? Wy 
1 1/3 1/3 1/3 —9/16 
2 3/5 1/5 1/5 25/48 
3 1/5 3⁄5 1/5 25/48 
4 1/5 1/5 3/5 25/48 


1) 三 角形 面积 人 的 求法 
由 图 5. 3 知 


. On ° 
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Fu 一 Ela: ~a) jy — y; t K(z; — z; ) } (5.46) 

Fe = iCxm — L) + j(y,, — yr) + k(x,, — z,) K 
Fy Tm 一 | Fa | | rm | ° cos(x — 0) 

(5. 47) 


ee [rs | [Fim | cos 


IP, ry WAER j 节点 指向 大 W ra AE k 节点 指向 m 节点 ，6 


Kry r, 53248. 故 
(=, — =. )(m, — zx) + (3 yt ,) T (zm —z;)(z2, — z, ) 






















iii (z, — z)! F OQ — y Te) Van — m): + (y, A) + (z, — z): 
sin0 = V1 一 cos20 
三 角形 单元 的 高 
= w (ay —a y + (Qy — 3) + Can — zr“ sing (5. 48) 
三 角形 面积 
~ ll |ru | lr sing = $ lru | rl/1— TTT) 
=+ "Ñ 二 [二 
(5. 49) 
2) |r| 的 求法 





zü (5.45) PEI ri, rs rss r 的 表达 式 可 与 为 
F, = W (m, — m,y =F (y, — y.) -F Ça, — @;) 
AHP, q HARAT AR, ra 是 区 域 0 内 茶点 (zm, y> z) 至 单元 上 任 一 点 
(z, Yo z) 的 距离 。 
我 们 知道 三 角形 单元 上 的 坐标 可 表示 为 
x = LEP + rE 十 ones 
y = yP F mE + yném (5. 50) 
z = Zé H z£ F znén 
对 于 四 点 (q=1, 2, 3, 4) 高 斯 积分 公式 ， 当 q=1 时 有 








rn = V (m, ri) + (y, — y;)° +- (z, — z; (5.51) 
式 中 ， 
A (1) (1) (ly — 1 I j! 1 
Tı 一 j$ + rêi ga B T+ + 
_ e (1) a dl ] i 
yi = YE + y& T yném = ya s gY + 3 Ym 
z = zi6f + xD H znt = r * tatiz 


. 97 ° 
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4 q=2, 3, 4 时 有 








r = B =m Tti — 一 (5. 52) 
式 中 ， 
而 == Ez Hir 十 ==, 
3 1 1 
3 = EI Ty F-E Ym 
zz = Žr tien Hirn 
rs = y (zs — x) + Oa — y) + (zs — z) (5. 53) 
式 中 ， 
T; = ==; += +=, 
1 3 ] 
ya 一 5 2; kr TT Ym 
z = Fat n t en 
ra = V (m, — xi) + (y, — y;)° + (x, — z;)° (5. 54) 
式 中 ， 
Ia = ==; += 一 S=, 


1 


| 1 3 
Y 一 BY HEY FE Ym 


z = Sr tE Hen 


. OR 。 
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随 着 TEM 虚拟 波 场 深度 偏 移 方法 的 应 用 ， 人 们 开始 关注 直接 影响 仿 移 成 
像 质 量 的 偏 移 速度 建 模 技 术 。 在 地 震 勘 探 中 ， 速 度 建 模 技 术 可 分 为 射线 法 和 波 
动 方程 法 两 种 方法 ， 波 动 方程 法 又 分 为 深度 聚焦 分 析 和 剩余 曲率 分 析 ， 其 中 的 
射线 法 是 基于 射线 追踪 理论 。 按 所 基于 的 道 集 ， 分 为 共 中 心 点 道 集 、 共 成 像 点 
道 集 、 共 反射 点 道 集 、 共 聚焦 点 道 集 和 共 反 射 面 道 集 。 这 些 速度 分 析 方 法 在 
TEM 虚拟 波 场 深度 偏 移 中 均 可 借鉴 ， 但 需要 指出 ，TEM 虚拟 波 场 是 经 过 数 
学 变换 得 到 的 ， 不 是 客观 存在 ， 因 此 ， 还 有 其 特殊 性 。 本 文 从 瞬 变 电磁 场 本 喘 
` 的 特性 出 发 ， 基 于 等 效 导 电 平 面 的 理论 ， 提 出 瞬 变 电磁 虚拟 流 场 连续 速度 分 析 
方法 。 


6. 1 基于 等 效 导 电 平 面 法 的 速度 建 模 原理 


等 效 导电 平面 法 是 根据 视 纵向 电导 曲线 的 特征 值 直观 地 划分 地 层 的 一 种 近 
似 方法 ， 因 此 又 称 “ 视 纵向 电导 法 ”。 由 于 该 方法 可 以 形象 地 理解 为 ， 随 看 t 


的 增 减 ， 等 效 导电 平面 以 -速度 上 下 “浮动 >， 这 样 便 可 以 用 一 块 随时 间 + 的 


变化 而 “沉浮 ”的 “ 载 流 ”导电 平面 来 近似 代替 回 线 源 中 的 均匀 大 地 ， 以 便于 
求 出 地 表 任 一 点 的 异常 场 。 可 见 ， 导 电 平 面 上 下 浮动 的 速度 仅 与 介质 的 电导 率 
有 关 ， 这 一 点 与 TEM 虚拟 波 场 的 传播 速度 完全 一 致 。 可 以 设想 ， 用 等 效 导电 
平面 法 求 出 地 电 断 面 上 的 纵向 电导 分 布 ， 即 求 得 随 深 度 变 化 的 纵 同 电导 值 ， 由 
TEM 虚拟 波 场 的 速度 公式 计算 出 速度 值 ， 实 现 连续 速度 分 析 。 


6.1.1 等 效 导 电 平 面 法 基本 原理 


设 在 地 面 上 有 圆 形 回 线 ， 其 中 通 以 阶 跃 电流 ， 并 在 £= 0 时 刻 断 开 电 
源 ， 即 
I t<0 
0 £—ÜO 
则 在 上 过 0 时， 在 地 下 产生 涡 旋 电流 ， 在 地 表 任 一 点 均 可 观测 到 由 此 涡流 
产生 的 电磁 场 。 根 据 电磁 理论 ， 我 们 可 用 一 个 导电 平面 来 代 蔡 地 下 均匀 介质 ， 
然后 用 镜像 法 可 方便 地 求 出 空间 任 一 点 的 感应 电磁 场 。 如 图 6.1 所 示 ， 设 A 


. 909 . 


I(t) = 


8 S BOE R 585 pR: 


为 一 完全 导电 平面 ， 面 电导 率 为 as， 其 hh 处 有 一 条 半径 为 a 的 圆 形 回 线 ， 其 中 
EMREN I 





h 
| i 
| | | 4 平面 
| je 
Í r | 
图 6.1 等 效 导电 平面 图 


当 电源 断 开 时 ， 由 于 电磁 感应 ， 导 电 板 A 上 产生 涡流 ， 为 求 出 涡流 产生 
的 电磁 场 ， 在 A 前 方 h 处 ， 置 一 虚 源 代替 导电 平面 板 中 的 涡流 ， 因 此 只 需求 
出 虚 源 在 空间 任 一 点 产生 的 电磁 场 即 可 。 

根据 电磁 场 理 论 ， 引 入 矢 势 4， 则 导电 板 内 任 一 点 的 电场 强度 为 





4 
板 内 面 电 流 密度 为 


根据 板 内 涡流 的 分 布 可 知 ， 矢 势 A X A, DE, FAE r 的 函数 ， 而 且 ` 
WREED E r ARK. ABRE E MRE 1(7) 为 板 上 任 一 点 7 处 到 边界 上 
的 总 电流 ， 则 面 电流 密度 为 


. _ dI(r) 
J dr 
由 安培 环 路 定律 
tol(r)=4B.dl (6. 3) 


L 
积分 路 径 L 为 自 r 经 平板 上 方 绕 过 板 边缘 ， 再 由 下 方 返 回 > 处 。 由 式 
(6. 3) 得 | 





.— dl) dfr. 1 d pi —— 
poj = p SY = 48 dl 21] Roda =s— BIR 
故 
FE (6.4) 
pto 
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将 式 (6.4) RAR 6.2), A 
2Br -2A。 





Om to 3t (6.5) 
HT B=VXA, 将 B= — ÜA (6.5) 式 ， 有 
= = -S (o = om) (6. 6) 
该 方程 的 解 有 如 下 形式 
A, 一 了 (Z 十 7) (6.7) 


这 就 是 板 外 任 一 点 处 矢 势 4 的 A。 分量 所 满足 的 方程 。 

根据 镜像 法 原理 ， 像 源 线圈 在 地 面 上 一 点 产生 的 矢 势 分 量 A。( 即 异常 场 ) 
可 由 下 式 表 示 

A, = f (z+ -)_ = (2 + t | (6. 8) 

此 式 说 明 ， 在 t>0 时 ， 虚 源 以 速度 2 回 下 作 勾 速 运动 。 可 形象 地 理 
解 为 ， 随 t 的 增 减 ， 等 效 导电 平面 A 便 以 速度 1/jwo 上 下 “浮动 ?>。 当 上 增 大 
时 ， 它 下 沉 ， 当 1 减 小 时 ， 它 又 逐渐 “上 浮 ”。 这 样 便 可 以 用 一 块 随 时 间 上 的 
变化 而 “沉浮 ”的 “ 载 流 ”导电 平面 来 近似 代 蔡 大 地 的 均匀 介质 ， 从 而 便于 求 
出 地 表 任 一 点 的 异常 场 。 


6.1.2 水 平 层 状 介质 表面 的 瞬 变 电磁 场 的 近似 计算 


有 瞬 变 电 磁 法 中 ， 视 纵向 电导 指 的 是 在 断 电 后 的 某 一 时 刻 ， 电 磁场 的 有 效 分 
布 范围 内 地 电 断 面 能 被 具有 一 定 深度 的 均匀 导电 平面 等 效 代替 时 的 纵向 电导 。 
并 由 此 可 知 ， 该 等 效 导 电 平 面 和 其 纵向 电导 以 及 等 效 深 度 均 随时 间 而 改变 。 

下 面 导 出 均匀 大 地 表面 上 回 线 中 心 处 的 异常 场 表 达 式 。 

选取 柱 坐 标 系 ， 设 虚 源 中 的 电流 强度 为 I， 虚 源 平面 与 zoy FHES., A 
心 位 于 坐标 原点 。 任 取 一 线 电流 源 Ids, CE z 轴 上 任 一 点 P 处 产生 的 磁场 强 
度 垂直 分 量 为 


i= 二 S sing (6.9) 
则 整个 回 线 在 P 的 
j° 1 m: 
H, -ppi ° FEE (6. 10) 
考虑 到 等 效 导 电 平 面 的 “浮动 ”， 把 z= 二 2h 十 2t/joo RAER, 4 
.et (6.11) 


: [eraat] 
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于 是 
| 有 十 一 


APF, s= 为 等 效 导 电 平 面 的 电导 率 。 
式 (6.11) MA (6.12) 即 为 用 等 效 导电 平面 法 计算 均匀 大 地 表面 上 回 线 
中 心 处 瞬 变 场 的 近似 计算 公式 。 





A 
m = ES 
pos 
= ow í p 
HL — = 
a posa 
则 
人 | 
Qa (aiidm) 2a (1 + dm 3222 (6.13) 
ƏB,_ Sl eam  __6l m _ (6.14) 
dt sa? (a2 +4) sa? (1 + 4m°)5/ 9 
令 4221 EË a2<4m ， 即 得 到 晚期 瞬 变 场 的 近似 表达 式 
2 
H, = -2 (6.15) 
2 
S =a 2a (6. 16) 





对 于 水 平 层 状 地 电 断 面 ， 我 们 也 可 以 用 等 效 导 电 平 面 法 求 得 瞬 变 场 的 近似 . 
解 。 只 不 过 这 时 的 等 效 平面 随时 间 增 大 而 “下 沉 ” 的 速度 不 再 是 常数 ， 而 分 别 
与 各 层 介 质 的 电导 率 有 关 。 

对 于 水 平 层 状 介 质 ， 引 入 等 效 导 电 平 面 的 深度 hj4K， 它 等 于 组 成 断面 的 单 
位 薄 层 的 平均 厚度 。 


H r 
i c (z)zË dz 


hy. 一 





rH $ 
=: l anamda (6.17) 


Š 








[odz 
AP, o (z) 为 介质 电导 率 ; H 为 研究 深度 ， 即 某 一 时 刻 对 异常 场 作 实 质 贡 
献 的 深度 ; s 为 全 究 深度 以 上 介质 的 总 纵 癌 电导 ; g 为 在 计算 hx 过 程 中 ， 确 
定 上 、 下 地 层 相 对 权 的 参数 ，g 盖 0; O0<hx< H, H hy 等 于 等 效 导 电 板 的 中 
DORE., 


由 式 (6.14) 4, EPEHA D 


的 近似 公式 为 


. 102 ° 
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3 B-t) __ 6I m f 
at sa: (1+ 4m2 )5/2 (6. 18) 
将 该 式 写 成 如 下 形式 
ET (6. 19) 
£ S 
. sbl wa m 
式 中 ， p= z ° F (m) (Fim y” 
我 们 由 上 式 引入 视 纵向 电导 
BF) x 
w F) 3 B. (t) (6. 20) 
dt 


该 式 可 用 于 野外 实测 资料 的 整理 ， “为 实测 值 。 


6.1.3 元 参数 的 优化 提取 
由 式 (6. 20》 知 ， 在 获得 实测 “二 值 后 ， 要 求 得 S. (1) ， 必 须 得 到 F 
关 。 野外 


(m) 的 值 ， 而 (Cm) 与 等 效 导 电 平 面 的 深度 hh 有 关 ， 即 与 研究 深度 有 关 
一 3 全 曲线 进行 微分 的 广 


实际 工作 中 难以 直接 得 到 严 (元 ) 值 。 为 此 ， 利 用 对 了 3 (D 
法 间接 求 取 下 (元 ) 值 。 对 于 "qh 求 其 对 时 间 的 导数 

















2 ° B, = -一 (< 9B.\ 9 m) _ k _ ƏF(m) , at 
ar am\ 8! 3t s 3m 9 (m) 
Š, IF, (6.21) 
S 9 m posa ' 
5| A 388 BJpR38 olm) 
(6. 22) 
(6. 23) 
(6. 24) 





= — 
F (m) (1 + 4m 572 
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于 是 





pm)= (1 + dm 73 4—16] (6. 25) 
综 上 所 述 ， 求 取 纵向 电导 曲线 的 具体 步骤 为 ; 
(1) 首先 求 得 一世 实测 曲线 和 它 的 时 间 导数 曲线 ; 
9 £: X 





(2) 把 全 二 和 一 二 学 呈 代入 式 6.22)， 求 得 辅助 函数 ， 
(3) kinai t 的 gp(m) 值 对 应 的 值 m; 
(4) 把 m 代 入 式 (6.24)， 再 代入 式 (6. 20) 求 得 S. (i) 的 值 ， 并 代入 下 
式 求 H. (1)， 其 中 
H, (t)= (m—t/pos.) ($m =)" (6. 26) 
由 上 述 步 骤 可 知 ， 要 求 得 视 纵 向 电导 S. 和 视 深 度 H., WLR m, IQ 
m 值 可 由 式 (6.25) 中 的 gC(m) 求 得 ， 且 该 式 中 的 op(m) 的 大 小 只 与 mx 值 有 
关 ， 我 们 可 把 该 式 所 求 得 的 gp(m) 值 视 为 理论 的 olm), WE pr (m). H 


(1) 中 的 一 二 二 是 一 组 野外 实测 获得 的 值 ， 时 间 t 在 测量 过 程 中 给 定 ， 
二 和 可 根据 已 获得 的 2 已 与 值 求 出 ， 由 此 我 们 可 把 由 式 6. 22) 得 


出 的 辅助 困 数 视 为 实测 的 pm), IE ps (m), 
在 此 推导 过 程 中 ， 式 (6.22) 中 的 p(m) 与 式 (6.25) 中 的 gCm) 是 相等 
的 ， 也 就 是 说 ，qs (mn) 应 等 于 gr (m)， 这 就 为 我 们 用 最 优化 法 求解 m 值 创造 
了 条 件 。 反 过 来 说 ，m 值 的 大 小 直接 影响 到 了 等 效 导 电 平 面 解 释 法 的 最 终结 
果 ， 即 视 纵 向 电导 S, AREH.. WR m ARTERA. s AH, 也 能 得 到 
最 优 解 ， 反 演 成 像 结 果 也 最 能 接近 实际 目标 体 的 模型 。 为 此 建立 目标 函数 
a2B 
ra 
g B. (t) 1° 
k mi 
当 该 函数 取 极 小 值 时 m 值 取得 最 优 解 。 
但 是 ， 大 多 数 的 优化 方法 都 要 先 给 出 一 初始 数值 ， 根 据 给 定 条 件 判断 该 初 
始 值 是 否 是 最 优 解 ， 如 果 是 就 终止 算法 ， 否则 对 初始 值 进行 改变 再 进行 运算 ， 
直到 符合 终止 条 件 为 止 ， 这 些 方法 的 初始 数值 模型 对 最 后 的 解 有 很 大 的 影响 ， 
而 好 的 模型 的 选择 却 又 是 另 一 个 较 难 解决 的 问题 。 当 然 也 可 采用 一 种 无 需 给 定 
初始 模型 的 优化 方法 来 解决 对 m 值 的 求解 。 如 模拟 人 脑 机 理 和 功能 的 一 种 新 
型 计算 机 和 人 工 智 能 技术 ， 或 者 建立 在 生物 进化 论 基 础 上 的 遗传 算法 。 
。 104 。 
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6.2 早 个 测 点 的 速度 建 模 万 法 


图 6. 2 是 基于 等 效 导 电 平 面 法 提出 的 虚拟 波 场 速度 估算 方案 ， 对 于 自 激 自 
收 观测 方式 ， 我 们 从 最 顶层 的 估算 开始 ， 给 出 初始 参数 m”， 便 可 确定 ss” 和 
hr ， 给 原始 模型 参数 mx 一 个 小 的 扰动 bmx， 然 后 再 利用 牛顿 迭代 法 计算 模型 
参数 更 新 


w ent di (6. 27) 






最 顶层 的 原始 模型 
x \ 下 一 次 迭代 


= 
计算 参数 更 新 向 量 | 更 新 模型 
ëm. ôv | w)i0+ëv 





1 拾取 下 一 个 同 相 轴 


估算 的 速度 模型 
图 6.2 速度 模型 选 代 算法 


是 
下 一 个 后 


得 到 新 的 模型 后 ， 再 由 目标 肾 数 B(m) 的 极 小 值 确 定 m 参数 的 一 次 估计 值 
m, HA (6.20) 得 到 s Ah. 的 一 次 信 计 值 s:” 和 h. ， 重 复 以 上 步骤 便 可 


得 到 一 组 随 视 深 度 变化 的 视 纵向 电导 分 布 曲线 。 求 得 随 深 度 变化 的 纵向 电导 值 
后 ， 由 总 纵向 电导 公式 


N 
S (A) = > ai (6. 28) 
1 三 ] 
可 以 导出 


_ S.(h,) ~S. Chi) 
i 


即 可 求 出 电导 率 值 c， 然 后 由 TEM 虚拟 波 场 的 速度 公式 


(6. 29) 
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v = 1/Vuo (6. 30) 
计算 出 速度 值 ， 实 现 速度 建 模 。 


6.3 连续 速度 分 析 


对 于 从 实际 测量 的 数据 中 得 到 的 虚拟 速度 值 w， 因 为 实际 测量 的 数据 量 不 
可 能 很 大 ， 所 以 计算 得 到 的 速度 值 在 空间 分 布 上 也 会 表现 得 较为 稀少 ， 这 样 就 
会 给 进一步 的 波 场 延 拓 成 像 工 作 带 来 不 便 。 为 了 解决 虚拟 速度 数据 不 够 丰富 的 
问题 ， 我 们 将 使 用 三 维 空间 插值 计算 的 办 法 ， 来 使 该 数据 量 在 保证 正确 的 基础 
上 得 以 迅速 扩大 ， 而 在 多 种 空间 插值 方法 中 ， 我 们 应 该 以 如 下 的 过 程 来 选择 适 
合 我 们 使 用 的 方法 (图 6. 3) 。 


内 插 方 法 选择 





| 根据 交叉 验证 统计 指标 


| 对 空间 数据 进行 估计 


内 插 方 法 评价 | 


重 
新 
选 
择 
内 
插 
方 
法 





图 6. 3 空间 内 插 方法 选择 流程 图 
最 终 ， 根 据 在 正确 性 和 速度 上 的 双重 标准 ， 我 们 选择 以 如 下 的 两 步 插值 方 
法 来 进行 计算 。 
6.3.1 全 局 距离 加 权 插 值 


及 用 全 局 距离 加 权 插 但 的 办 法 ， 将 计算 得 到 的 原始 虚拟 速度 值 进 行 适度 数 
据 量 的 插值 。 所 谓 全 局 距离 加 权 插 值 ， 是 指针 对 散乱 数据 的 全 局 按 距 离 成 反比 
的 加 权 插 值 计算 ， 其 速度 受 原始 数据 点 数目 和 插值 点 数目 的 影响 。 
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— 


假设 平面 上 分 布 一 系列 离散 点 ， 已 知 其 位 置 坐标 o y) 和 属性 值 P 
(i=1, 2, n), p (z, y) 为 任 一 格 网 点 ， 全 局 距离 加 权 插 值 法 假设 P 点 
的 属性 值 是 在 局 部 邻 域内 中 所 有 数据 点 的 距离 加 权 平 均值 ， 可 以 进行 确切 的 或 
者 平滑 的 方式 插值 。 周 围 点 与 P 点 因 分 布 位 置 的 差异 ， 对 P (=) 影响 不 同 ， 
我 们 把 这 种 影响 称 为 权 肾 数 W, (zx，y)， 其 主要 与 距离 有 关 ， 男 外 也 与 方向 
有 关 ， 这 时 


| _ i 
i=] [d; (x,y) J 





P(z) = I (6.31) 
2 [d,(z=,y) ]” 
式 中 ， 
d.(z=,y) = / (zw — x) + (y — y) (6. 32) 


它 表 示 由 离散 点 (x;，y;) 2 P (=, y) 点 的 距离 。P (z) 为 要 求 的 待 
捅 点 的 但 。 

RAR W, (z, y) 与 预测 点 和 已 知 样 点 之 间 的 距离 的 zx KERRE. A 
此 ， 随 大 距离 的 增加 ， 权 重 迅 速 地 减 小 ， 权 重 减 小 的 速度 取决 于 指数 的 大 
小 。 对 于 一 个 较 大 的 指数 ， 较 近 的 数据 点 被 给 定 一 个 较 高 的 权重 份额 ;对 于 一 
个 较 小 的 指数 ， 权 重 比较 均匀 地 分 配给 各 数据 点 。 

全 局 距离 加 权 插 值 方法 创建 一 个 粗略 的 数据 表面 ， 该 数据 表面 受 奇异 值 
(特别 是 边缘 的 极 高 、 极 低 值 ) 的 影响 较 大 。 同 时 ， 网 格 化 的 朴 密 程度 也 对 结 
果 有 影响 ， 如 图 6.4 ( 见 彩 图 ) 所 示 ， 图 6.4 (a) AX, Z 方向 网 格 30 的 插 
值 结果 ,图 6.4 (b) 为 X，Z 方 向 网 格 100 的 插值 结果 。 可 以 明显 看 出 下 图 


的 图 像 非 常 圆 请 ， 设 有 过 多 的 锯齿 图 案 。 
从 Zim 





ys — === 


09 s. 





(b) 
图 6.4 全 局 距离 加 权 插 值 的 不 同 网 格 化 程度 的 对 比 
(a) X，Z 方 向 网 格 30 的 插值 结果 (b) X，Z 方 向 网 格 100 的 插值 结果 
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采用 这 个 方法 做 第 一 步 插值 ， 主 要 是 为 了 保证 在 扩大 数据 量 的 同时 ， 还 保 
证 在 插值 出 来 的 点 上 ， 虚 拟 速度 值 与 周围 原始 数据 过 渡 圆 消 ， 其 什 与 只 实 但 从 
可 能 接近 。 而 只 做 小 范围 内 的 计算 ， 是 为 了 不 让 插值 点 数 过 于 超过 已 知 点 数 
目 ， 导 致 假 值 的 过 多 出 现 ， 并 且 不 会 使 计算 速度 过 慢 ， 时 间 过 长 。 


6.3.2 近 点 线性 插值 


对 由 全 局 距离 加 权 插 值 出 来 的 数据 ， 进 一 步 做 近 点 线性 择 什 计算， 根据 我 
们 的 实际 需要 ， 在 我 们 需要 的 空间 点 上 ， 插 值 出 虚拟 速度 值 。 

首先 需要 对 现 有 的 源 数据 进行 网 格 化 ， 为 了 配合 延 拓 成 像 工 作 ， 我 们 将 采 
取 曲 面 延 拓 一 样 的 网 格 划分 方法 ， 使 用 三 角 网 格 。 如 图 6. 5 所 示 ， 图 中 普通 数 
字 为 网 格 节点 编号 ， 圈 型 数字 是 单元 编号 。 
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图 6.5 ” 近 点 线性 插 伍 网 格 划分 方法 示意 图 
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由 此 可 以 得 到 原始 文件 格式 ， 分 别 给 出 了 各 个 已 知 点 在 地 面 上 的 坐标 和 在 
地 下 的 省 度 以 及 该 点 上 的 电 性 参数 数值 。 

紧 接 者 ， 我 们 要 确定 捅 值 点 的 邻 域 。 空 间 位 置 相 近 的 两 个 样品 点 间 存 在 痢 
一 定 的 相关 性 ， 但 随 着 距离 的 增加 这 种 相关 性 将 变 弱 。 在 一 定 的 距离 外 ， 采 样 
点 与 预测 点 不 相关 。 因 此 在 实际 应 用 中 ， 需 要 通过 定义 邻 域 的 形状 与 邻 域内 已 
知 样品 点 的 数目 。 

一 个 好 的 预测 方案 要 求 参 与 预测 的 观测 点 来 自 多 个 方 申 上 。 如 图 6.6 所 
示 ， 一 个 圆 形 邻 域 被 分 成 四 个 局 区 并 定义 每 个 山区 内 参与 预测 的 观测 点 的 数目 
为 3 一 5 个 。 


A a Zk waska 
e@ q IDTVTpeoeoeooeose 
AN | esen. 
图 6.6 邻 域 中 参与 预测 的 观测 点 数目 的 定义 
取样 点 的 位 置 在 空间 插值 有 着 重要 的 影响 。 当 数据 点 相关 且 均 习 分 布 时 ， 能 
更 好 地 反映 研究 要 素 在 空间 的 分 布 特征 ， 对 于 源 数据 的 分 析 和 插值 方法 的 选择 有 很 
大 的 帮助 ; 当 数 据点 成 群 或 不 均匀 分 布 时 ， 难 以 很 好 地 反映 研究 要 素 在 整个 空间 的 
分 布 特征 ， 可 能 会 造成 插值 结果 的 不 合理 。 取 样 时 ， 尽 量 使 样本 点 均匀 分 布 。 
如 图 6.7 〈 见 彩 图 ) 所 示 ， 对 于 要 捅 值 的 每 一 点 ， 在 其 上 方 的 所 有 已 知 点 





| ” Fo 
Fa | 
图 6.7 近 点 线性 插值 计算 原理 示意 图 
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都 要 对 这 一 点 的 速度 值 产 生 影响 ， 因 此 针对 任意 一 个 插值 点 ， 地 表 上 每 一 个 已 
车 氮 都 将 和 它 构成 一 个 选择 预测 点 的 邻 域 ， 为 了 确定 每 一 次 计算 时 的 邻 域 大 
小 ， 我 们 将 采用 通过 两 点 来 确定 一 条 直线 的 办 法 ， 把 地 表 已 知 点 与 插值 点 的 连 
线 计算 出 来 ， 则 以 该 连 线 作为 对 角 线 的 立方 体 即 为 所 需要 的 邻 域 ， 该 区 域 的 每 
E igp ag 
近 的 一 点 数据 作为 插值 点 的 临时 值 记录 下 来 ; 然后 换 一 个 地 表 已 知 点 重复 这 个 
步骤 ， 直 到 地 表 的 已 知 点 数据 全 部 做 完 ; 最 后 ， 对 所 有 的 临时 值 取 平 均 ， 所 得 
到 的 数据 就 作为 插值 点 的 最 终结 果 保存 进 结果 文 件 ， 这 样 ， 就 完成 了 一 次 近 点 
线性 插值 计算 。 

在 计算 过 程 中 ， 根 据 格林 公式 的 原理 ， 为 了 保证 延 拓 时 每 一 次 单元 积分 都 
泊 逆 时 针 方向 ， 在 捅 全 过程 中 ， 地 表 已 知 点 的 循环 顺序 必须 严格 按照 节点 编号 
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图 6.8 近 点 线性 插值 计算 循环 顺序 图 
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在 上 面 的 计算 过 程 中 ， 关 键 的 步骤 就 是 确定 邻 域 的 对 角 线 ， 即 算出 地 表 已 
MASER LEEREN E. 

如 图 6.9 BUR, a, b, m 三 点 的 空间 坐标 为 已 知 ， 那 么 ， 可 以 写 出 6，m 
或 者 a，m 之 间 的 直线 点 回 式 方程 

T— To _ yy _ z 
m n p 
Ms Ns b 为 方 回 数 ， MSX To» NS Y2 Yos P=: z, Zo. 
由 以 上 方程 可 以 得 到 直线 的 参数 方程 
-+ 





(6, 33) 


y = yo + nt (6. 34) 


z = zo + bt 









a E A A A E 
EERE oa 
£ R L Z 





b(x;,y EE 0) 


>, A A 


A m(xe.yo, zo) 
图 6. 9 近 点 线性 插值 计算 方法 示意 图 


对 于 波 场 延 拓 需要 的 数据 点 来 说 ， 可 以 事先 规定 好 每 次 向 下 延 拓 的 距离 ， 
则 由 此 就 可 以 知道 z 值 ， 将 该 值 代入 直线 参数 方程 中 的 第 三 个 方程 ， 即 可 求 得 
t 值 ， 表 把 t 值 代入 前 两 个 方程 中 ， 就 可 以 把 需要 点 的 7，y 值 求解 出 来 。 

然后 我 们 就 可 以 在 以 5 点 与 m 点 所 确定 的 直线 为 对 角 线 的 立方 体 范 围 内 ， 
寻找 与 所 求 点 空间 距离 最 近 的 一 个 已 知 点 ， 用 该 点 的 虚拟 速度 值 代替 需求 点 的 
速度 值 。 

最 终 ， 可 以 得 到 如 图 6. 10 的 近 点 线性 插值 的 流程 图 。 

这 种 三 维 空间 插值 方法 ， 它 的 精度 较 高 ， 主 要 受 原 始 数 据 的 数据 量 多 少 
影响 比较 大 。 当 然 ， 这 种 插值 方法 本 质 上 仍然 属于 一 种 线性 插值 ， 相 比 于 前 
面 介绍 的 样 条 函数 插值 法 以 及 克 里 金 〈Kriging) 插值 法 等 方法 ， 显 然 精 度 上 
要 低 一 些 , 但 其 优点 是 插值 所 耗费 的 时 间 非 常 少 ， 在 保持 较 高 精准 度 的 同 
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时 ， 还 保证 了 很 快 的 计算 速度 ， 从 而 大 大 提高 了 这 种 方法 应 用 于 实际 的 可 
能 性 ， 
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更 换 插 信 点 | 





_ 给 出 连续 速度 分 析 _》 
图 6. 10 近 点 线性 插值 算法 流程 图 

从 图 6.11 〈 见 彩 图 ) 可 以 看 出 ， 下 边 的 插值 结果 很 好 地 保持 了 上 边 的 原 
始 数据 的 基本 形态 ， 并 且 对 细节 的 反映 更 加 清晰 一 些 。 

为 了 更 深入 地 了 解 该 插值 方法 的 准确 性 ， 我 们 再 对 原始 数据 体 和 插值 后 的 
数据 体 进 行 切片 对 比 ， 如 图 6. 12 所 示 ( 见 彩 图 )。 


° 112 ° 


第 6 章 TEM BE 13 BF yr 


. Zim 








P o S 
. € %@ ° 


Ym ` 
(b) 
图 6.11 近 点 线性 插值 结果 对 比 图 
(a) Fins: (b) MARE 
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Y/m ` 
图 6.12 近 操 线性 插值 结果 切片 对 比 图 ( 续 ) 


从 图 6.12 可 以 看 出 ,在 X, Z 方向 的 切片 图 中 ， 插 值 数据 的 形态 基本 相 
同 ， 数值 上 也 和 原始 数据 几乎 完全 一 致 ， 而 在 Y 方向 的 切片 图 中 ,插值 数据 
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的 形态 仍然 与 原始 数据 相近 ， 但 在 数值 上 有 了 一 定 的 偏差 .这 是 需要 进一步 改 
进 的 地 方 。 

最 后 ， 我 们 来 对 这 种 插值 方法 做 一 个 有 效 性 评价 ， 本 次 评价 将 采用 平均 估 
计 误 差 百 分 比 PAEE 作为 交叉 验证 统计 指标 。 

PAEE = 100% pA Z. Ge y (6. 35) 

AP, Z 为 位 置 zx 随机 变量 Z 估计 值 ，Z; (z) AME zs EWS, Z, 
为 样 点 平均 值 ，n 为 样 点 数 。 

如 果 PAEE 趋向 于 0， 则 认为 估计 是 无 偏 的 。 

经 过 计算 ， 这 种 方法 的 平均 估计 误差 百分比 PAEE 为 1.07% ， 该 误差 值 
是 很 小 的 ， 因 此 可 以 认为 该 方法 是 切实 有 效 的 。 
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为 了 更 好 理解 瞬 变 电磁 法 拟 地 震 解释 技术 及 偏 移 成 像 效 果 ， 本 章 给 出 了 对 
瞬 变 电磁 理论 模型 进行 偏 移 成 像 的 数值 模拟 结果 ， 以 及 作者 近年 来 的 工作 
实例 。 


7. 1 模型 计算 
7.1.1 BRAE 


在 地 表 布 置 面积 工作 如 图 7.1 所 示 ， 共 布置 11 条 剖面 ， 每 条 剖面 布置 
11 个 测 点 ， 取 发 射 回 线 边 长 为 200m， 时 窗 采 样 宽度 为 32.5us—8.7ms, X} 
全 部 121 个 测 点 进行 正 演 计算 ， 对 所 得 电磁 响应 进行 波 场 变换 ， 得 到 各 测 点 
的 虚拟 波 场 值 ， 以 数据 文件 的 形式 输入 计算 机 ， 计 算 地 电 断 面 上 的 波 场 延 
拓 值 。 





”地 下 电 性 界面 


图 7.1 计算 模型 示意 图 


朋 变 电磁 中 心 回 线装 置 下 的 一 个 三 层 A 型 地 电 模型 的 参数 为 : p 一 500 + m, 
二 5000。m，m 三 50000。m， 太 二 100m， 和 和 200m。 在 这 个 模型 条 件 下 ， 
对 每 个 测 点 进行 正 演 计算 ， 得 到 电磁 响应 值 ， 之 后 进行 波 场 变换 ， 为 了 清楚 起 
见 ， 只 对 变换 后 的 波峰 进行 显示 ， 如 图 7. 2 所 示 〈 见 彩 图 )。 从 图 中 可 以 看 出 ， 
随 着 虚拟 时 间 z 的 增 大， 虚拟 波 场 出 现 了 两 个 正 的 峰值 变化 ， 这 两 个 峰值 就 标 
示 出 了 这 个 A 型 地 电 模 型 的 两 个 地 下 分 界面 的 位 置 。 利 用 每 个 测 点 的 电磁 响 
应 值 ， 用 = -三 一 计算 虚拟 速度 值 。 速 度 计算 的 结果 如 图 7. 3 所 示 ( 见 彩 

v Ho 
图 )。 可 以 看 到 ， 虚 拟 速度 的 分 布 明显 呈现 三 层 分 布 ， 与 模型 的 电阻 率 值 分 布 
e 116 = 
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是 对 应 的 ， 这 是 因为 虚拟 波 场 的 速度 值 取决 于 地 下 电导 的 分 布 ， 因 此 ， 虚 拟 速 
度 表 现 出 与 电阻 率 一 致 的 三 层 分 布 规律 。 





图 7.3 A 型 地 电 模 型 波 场 速度 分 析 结 果 


最 后 ， 将 三 维 边界 元 拉 术 应 用 于 基 尔 和 霍 夫 积分 计算 中 ， 利 用 速度 分 析 结 采 
实现 三 维 延 拓 成 像 ， 得 到 图 7.4( 见 彩 图 ) 所 示 的 结果 ， 两 个 峰值 明显 表示 出 
了 A 型 地 电 模 型 的 两 个 地 下 分 界面 的 位 置 。 这 两 个 位 置 所 对 应 的 座 度 也 恰好 
在 100m 与 300m 的 附近 ， 与 模型 的 参数 基本 吻合 。 
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图 7.4 A 型 地 电 模型 波 场 延 拓 成 像 结 果 


7.1.2 三 维 模型 


为 了 验证 本 文 提 出 方法 的 有 效 性 ， 设 计 了 半空 间 中 高 阻 块 状 异常 体 和 低 阻 块 
状 异常 体 的 三 维 模 型 ， 采 用 中 心 回 线 激发 ， 对 所 得 的 电磁 响应 进行 波 场 变换 ， 经 合 
成 孔径 计算 得 到 合成 波 场 值 ， 最 终 用 三 维 边界 元 法 实现 基 尔 霍 夫 偏 移 成 像 计算 。 

其 模型 三 维 示 意图 如 图 7.5 (a) 所 示 ( 见 彩 图 )， 异 常 体 为 30mX 30m X 
50m 的 块 体 ， 发 射 边 长 为 100m。 地 表 测 点 分 布 图 如 图 7.5 (b) 所 示 ， 图 中 所 
示 主 剖面 上 共 布 置 11 个 测 点 ， 点 距 为 10m, 





(a) 
图 7.5 模型 三 维 示 意图 
(a) 三 维 模型 立体 图 ; (b) 主 剖 面 平 面 图 
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距离 /m 





(b) 
图 7.5 模型 三 维 示 意图 ( 续 ) 


1. 高 阻 异常 体 

取 均 匀 半 空间 的 电阻 率 为 m 二 10Q0，。m， 高 阻 块 状 异 常 体 电阻 率 为 o= 
3000。m， 顶 板 埋 深 为 70m。 首 先 ， 用 正 演 方法 获得 了 主 放 面 上 的 视 电 阻 识 
断面 图 ， 如 图 7.6 所 示 ( 见 彩 图 )， 可 以 看 出 高 阻 异常 的 中 心 与 模型 的 位 置 并 


| 视 电 阻 率 /0 m) 
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出 现 了 向 上 偏 移 的 现象 ， 从 电阻 率 断 面 图 中 已 经 很 难 确定 异常 体 的 位 


不 对 应 ， 
置 和 大 小 。 


利用 瞬 变 电磁 合成 孔径 算法 ， 对 该 模型 计算 了 11 个 测 点 的 数据 ， 图 7. 7 


彩 


结果 , 图 7.7 (b) ( 见 


(a)( 见 彩 图 ) 为 合成 孔径 前 瞬 变 电磁 波 场 变换 计算 


AHH 
d BRDN ODA DOR OARS D O T 


g: FLLIL L LI1ITI- 


= s. 





p(w) i PIERE 


图 7.7 模型 计算 结果 


(b) 合成 孔径 后 波 场 变换 结果 


(a) 合成 孔径 成 像 前 该 场 变 换 绩 采 ; 
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图 ) 为 合成 孔径 后 瞬 变 电磁 波 场 变 换 计 算 结果 ， 对 比 这 两 个 图 ， 可 以 看 出 ， 在 
相关 性 较 好 的 中 心 部 分 合成 后 的 波形 得 到 了 一 定 的 增强 。 在 相关 性 不 好 的 边缘 
部 分 ， 波 形 得 到 罚 化 ， 合 成 孔径 的 有 效 性 得 到 验证 。 

最 后 ， 经 过 速度 分 析 、 波 场 基 尔 霍 夫 偶 移 成 像 ， 得 到 主 放 面 上 的 成 像 结 
果 ， 如 图 7.8 所 示 《〈 见 彩 图 )。 从 图 中 可 以 看 到 ， 上 下 两 组 波峰 分 别 对 应 在 
70m 和 120m 处 ， 这 和 模型 设置 吻合 很 好 ， 纵 向 分 辩 率 得 到 了 提高 。 第 二 界面 
波形 幅 值 相对 第 一 界面 减 小 了 很 多 ， 应 该 是 从 高 阻 到 低 阻 本 号 是 弱 反 射 面 上 所 
致 。 与 以 上 电阻 率 断 面 图 相 比 ， 它 具有 更 高 的 分 辨 率 。 
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WaS 
图 7.8 主 剖 面 偏 移 成 像 断 面 图 


2. 低 阻 异常 体 

均匀 半空 间 的 电阻 率 取 为 250。m。 低 阻 异 常 体 电阻 率 为 50，m， 模 型 大 
小 尺寸 不 变 ， 将 其 顶板 埋 深 加 大 到 120m。 主 剖面 上 的 视 电 阻 率 断面 图 如 图 
7.9 所 示 〈 见 彩 图 )， 可 以 看 出 低 阻 异常 的 中 心 与 模型 的 位 置 也 不 对 应 ， 出 现 
了 向 下 偏 移 的 现象 ， 从 电阻 率 断 面 图 确定 异常 体 的 位 置 和 大 小 也 会 出 现 仿 差 。 
经 波 场 变换 、 相 关 合 成 等 一 系列 处 理 ， 得 到 合成 前 后 的 主 剖 面 对 比 图 ， 如 图 
7.10 所 示 〈 见 彩 图 )， 可 以 看 出 ， 在 相关 性 较 好 的 中 心 部 分 合成 后 的 波形 得 到 
了 增强 。 在 相关 性 不 好 的 边缘 部 分 ， 波 形 得 到 弱化 ， 同 样 也 证 明了 合成 孔径 的 
有 效 性 。 
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图 7.10 模型 计算 
(a) 合成 孔径 成 像 前 波 场 变换 结果 ; (b) 合成 孔径 后 波 场 变换 结果 
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虚拟 时 间 /(ms) 





(b) 
图 7. 10 ”模型 计算 结果 (E) 


图 7. 11 (WEED 给 出 了 仿 移 成 像 结 采 ， 我 们 可 以 看 出 两 个 界面 的 存在 ， 
EE BATRA E E NAR: MEE FEA MWE 
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| 
图 7.11 主 剖 面 偏 移 成 像 断 面 图 
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7.2 应 用 实例 


7.2.1 隧道 超前 预报 实例 


锦 屏 水 电 枢 纽 工程 辅助 洞 位 于 四 川 省 凉山 青 族 自治 州 的 木 里 、 盐 源 、 晃 
宁 三 县 交界 处 的 雅 奢 江干 流 锦 屏 大 河 湾 上 ， 辅 助 洞 工程 地 处 青藏 高 原 向 四 川 
盆地 过 渡 地 貌 斜 坡地 带 ， 沿线 地 形 起 伏 ， 高 程 均 在 3000m 以 上 ， 最 高 山峰 
达 4125m， 埋 深 大 于 1500m 的 洞 段 长 度 约 12875m， 占 全 洞 长 度 的 73. 1%， 
地 质 条 件 非常 复杂 。 为 了 探 明 辅助 洞 的 水 体 不 良 地 质 体 ， 在 辅助 洞 不 同 的 掌 
+I E, 布置 了 四 条 剖面 [A 洞 掌 子 面 (S AK11 十 843)、A 洞 掌 子 面 
(里 程 AK11 十 447)、A 洞 掌 子 面 (里 程 AK11 十 598) 、 横 通道 (7 一 2#) ği 
面 」 通 过 对 锦 屏 水 电 枢纽 工程 辅助 洞 的 探测 实验 ， 验 证 本 项 目 提出 方法 的 有 
效 性 。 

瞬 变 电磁 超前 地 质 预 报 测 线 布置 如 图 7. 12 所 示 ( 见 彩 图 )， 在 一 个 框 内 
测量 5 条 测 线 ， 每 条 测 线 上 测量 3 一 5 个 点 ， 然 后 移动 线 框 ， 重 复 上 面 的 步 
又 ， 直 到 将 整个 掌 子 面 测 完 。 多 条 测 线 测 量 不 仅 可 以 保证 资料 的 充分 ， 同 时 
对 连续 速度 分 析 以 及 延 拓 成 像 来 说 也 是 确保 源 数据 可 以 构成 空间 目标 体 的 必 
要 条 件 。 





图 7.12 瞬 变 电磁 测 线 布置 图 
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首先 我 们 根据 实测 的 电磁 啊 应 人 计算 出 一 套 初 始 的 虚拟 速度 伍 ， 然 后 进 
行 两 次 三 维 空间 择 什 计算， 完成 速度 分 析 。 速 度 分 析 的 绪 末 如 图 7. 13 所 示 
( 见 彩 图 )。 可 以 看 到 ， 虚 拟 速度 的 分 布 大 致 呈现 为 三 层 分 布 ， 分 别 在 7m 与 
30m 左右 的 地 方 出 现 大 的 数值 变化 ， 其 规律 和 电阻 率 分 布 关系 是 相对 应 的 。 
这 是 因为 虚拟 波 场 的 速度 值 取决 于 地 下 电导 的 分 布 ， 而 电导 又 与 纵向 电导 密 
切 关 系 ， 纵 向 电导 与 电阻 率 为 倒数 关系 ， 因 此 ， 虚 拟 速度 表现 出 与 电阻 率 一 
致 的 三 层 分 布 规律 。 
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图 7. 13 ”虚拟 波 场 速度 分 析 三 维 结果 


接 下 来 对 采集 到 的 实际 数据 做 波 场 变换 ， 波 场 变 换 的 结果 如 图 7. 14 所 示 
( 见 彩 图 )， 可 以 看 到 ， 曲 线 充 分 显示 除了 波 的 特性 ， 出 现 了 多 个 波峰 和 波 合 。 
~ yim + pim 





-一 一 xm 


瞬 变 电磁 法 拟 地 震 偏 移 成 像 研究 





Ñ ~” <i yr 


图 7.14 实测 资料 虚拟 波 场 三 维 切 片 结果 ( 续 ) 


把 波 场 变换 的 结果 用 三 维 铁 空 技术 显示 ， 只 保留 与 水 体 不 良 地质 体 相对 应 
的 低能 部 分 ， 如 图 7. 15 所 示 〈 见 彩 图 )， 可 以 看 到 它 清 晰 的 反映 出 了 几 扶 务 
面 ， 且 其 形态 有 一 定 的 显示 。 





AFE i Sanka 
图 7.15 不 同方 位 水 体 不 良 地 质 体 波 场 变换 三 维 结果 


同样 ， 可 以 把 延 拓 的 结果 用 三 维 铁 空 方法 显示 ， 如 图 7. 16 〈 见 彩 图 ) 清晰 
的 在 7m 左右 位 置 与 20 一 35m 范围 内 反映 出 了 水 体 不 良 地 质 体 的 三 维 空间 分 布 。 

最 终 与 该 实际 测量 面前 方 已 知 的 情况 对 比 ， 可 以 确定 ， 这 个 经 过 一 系列 波 
场 变换 ， 速 度 分 析 以 及 延 拓 成 像 后 得 到 的 结果 与 实际 情况 吻合 得 很 好 。 


7.2.2 煤矿 采 空 区 探测 实例 


1. 实例 1 
在 煤矿 开采 过 程 中 ， 积 水 、 采 空 区 、 陷 落 柱 等 水 体 病害 一 直 困 扰 着 生产 。 
在 山西 某 煤 矿区 开展 采 空 区 调查 ， 目 的 是 要 搞 清 煤田 地 下 充 水 采 空 区 的 采 空 深 
度 、 采 空 高 度 、 分 布 范 围 ， 避 免 在 采掘 过 程 中 发 生 突 水 灾害 。 
勘查 区 表层 主要 被 第 四 系 黄土 覆盖 ， 厚 度 一 般 不 超过 20m， 下 伏地 层 主要 
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图 7.16 水 体 不 良 地 质 体 延 拓 成 像 三 维 结果 


为 第 三 系 砂砾 岩层 、 二 从 系 砂 泥岩 夹 煤 系 、 石 误 系 砂 泥岩 夹 煤 系 地 层 及 奥 陶 系 
谈 宕 地 层 ， 该 区 分 布 的 主要 煤层 有 两 屋 ， 分 别 位 于 地 下 150m 和 260m 人 处。 各 
地 层 岩 矿石 电阻 率 差异 较 小 ， 难 以 引起 明显 的 电阻 率 异 常 ， 当 地 下 煤 体 局 部 被 
来 出 后 ， 在 岩 体 内 形成 一 个 有 一 定 规模 的 空间 ， 使 周围 的 应 力 平 稳 状 态 遭 受 破 
坏 ， 产 生 局 部 的 应 力 集中 ， 采 空 区 顶板 在 上 有 覆 宕 层 压 力 的 作用 下 ， 发 生变 形 、 
断 型、 位移、 冒 落 ， 形 成 的 冒 落 市 、 断 裂 市 、 变 形 弯 曲 禹 ， 在 地 下 水 的 充填 太 
地 表 水 沿 裂缝 向 采 空 区 渗 漏 ， 其 电阻 率 将 明显 发 生变 化 ， 形 成 一 个 低 阻 电 性 体 ， 
也 与 围 岩 电 性 形成 较 明显 的 差异 ， 这 样 就 为 采 空 区 的 勘查 提供 了 前 提 条 件 。 

根据 该 工区 附近 的 采 空 区 的 埋藏 深度 和 厚度 情况 ， 结 合 有 瞬 变 电磁 特点 ， 选 
择 大 定 源 回 线 装置 ， 发 射 回 线 边 长 为 200m， 测 线 间 距 为 20m， 测 点 间距 为 
20m， 布 置 了 面积 工作 。 

(1) 波 场 变换 结果 

图 7.17 ( 见 彩 图 ) 是 经 波 场 变 换 ， 并 做 均衡 处 理 、 反 袜 积 、 多 点 平滑 后 最 
终 的 三 维 成 像 全 图 ， 波 场 值 颜色 深 的 部 分 为 虚拟 波 能 量 高 的 区 域 ， 也 是 反射 界 
面 所 在 处 。 
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(a) (b) 

图 7.17 波 场 三 维 成 像 效 果 图 

(a) 于 方向 切片 ; (b) Y 方向 切片 
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(2) 波 场 延 拓 成 像 结果 

为 了 更 好 地 展示 来 空 区 的 空间 分 布 ， 将 波 场 延 拓 成 像 的 结果 采用 两 种 显示 
方式 ,一 种 是 去 掉 低能 反射 (低速 或 者 低 阻 分 布 〉 部 分 ， 如 图 7. 18 (a) 所 示 
( 见 彩 图 )， 可 以 清楚 地 显示 采 空 区 在 上 下 围 岩 中 的 分 布 位 置 。 男 一 种 是 去 掉 高 
能 反射 (高 速 或 者 高 阻 分 布 ) 部 分 ， 如 图 7.18 (b) 所 示 ( 见 彩 图 )， 展 示 出 
了 采 空 区 的 空间 分 布 形态 。 通 过 对 成 像 结 果 进 行 解释 ， 认为 有 两 层 采 空 区 分 
布 ， 上 层 采 空 区 位 于 地 下 140m 处 ， 充 水 较 少 ， 冒 落 带 部 分 范围 小 ; FERE 
区 位 于 地 下 260m 处 ， 充 水 较 多 ， 冒 落 带 部 分 范围 大 ， 由 于 充 水 的 影响 ， 导 致 
低 阻 异常 放大 ， 采 空 区 的 分 布 范围 也 会 有 一 定 的 放大 。 经 验证 解释 结果 与 实际 
地 质 情 况 吻合 较 好 ， 这 也 证 明了 本 书 算 法 的 有 效 性 。 





(a) (b) 
图 7.18 三 维 延 拓 成 像 效果 
(a) 去 掉 低 能 反射 部 分 ! (b) 去 掉 高 能 反射 部 分 


2. 实例 2 

对 于 山西 很 多 地 方 煤矿 ， 在 煤矿 开采 过 程 中 ， 煤 矿 积 水 、 陷 落 柱 、 采 空 区 
等 矿井 地 质 灾害 长 期 以 来 一 直 困 扰 着 生产 ， 严 重 影响 着 矿井 的 生产 建设 。 需 要 
采用 瞬 变 电磁 法 对 煤矿 的 采 空 区 做 系统 的 研究 ， 探 明 其 范围 和 埋 深 .。 

工作 中 使 用 的 仪器 为 GDP 一 32II 电 法 工作 站 。 考 虑 到 工区 的 实际 地 质 情 
况 及 所 需 探 测 的 地 质 体 ， 采 用 瞬 变 电磁 法 进行 探测 。 工 作 中 根据 实际 试验 情 
况 ， 结 合 瞬 变 电 磁 特点 ， 选 择 采 用 大 定 源 回 线装 置 。 其 发 射 边框 采用 边 长 较 大 
的 和 矩形 回 线 ( 边 长 300m) ， 而 接收 采用 专用 的 瞬 变 电磁 探头 接收 〈 接 收 面 积 
2000m: ) 。 发 送 频率 为 25Hz， 时 间 范 围 为 0.087 一 8ms， 其 中 测 线 间 上 距 为 
20m， 测 点 间距 为 20m。 

图 7.19 a) 〈 见 彩 图 ) 是 视 电 阻 率 平 面 图 ， 由 图 可 知 ， 在 测 区 的 西北 部 、 
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东南 部 出 现 两 个 比较 明显 的 低 阻 异常 区 ， 推 测 是 由 于 采 空 区 充 水 引起 。 参 考 全 
区 视 电 阻 率 断面 图 成 果 ， 以 650 - m 等 值 线圈 定 两 个 低 阻 区 ， 即 A # B 区。 
图 中 红色 的 正方 形 框 为 三 维 成 像 数 据 体 采集 区 域 。 图 7. 19 (b) ( 见 彩 图 ) 为 
该 区 域 的 主 剖面 视 电 阻 率 等 值 线 断 面 图 ， 红 线 左 侧 为 三 维 成 像 数 据 体 采集 区 
域 ， 从 断面 图 中 可 以 看 出 ， 在 浅 层 ， 电 阻 率 较 高 ， 达 到 1000 ° m 以 上 。 从 
120m 左右 开始 ， 电 阻 率 快 速 下 降 至 150m 处 达到 极 小 ， 约 为 55Q，m。 之 后 
随 深度 增加 ， 电 阻 率 又 开始 缓慢 增加 ， 至 500m 又 恢复 到 1000 ° m 左右。 从 
图 中 可 以 隐约 看 出 两 个 电 性 分 界面 ， 第 一 层 位 于 130m 到 160m 深度 附近 ， 第 
二 层 界面 由 于 变化 缓慢 不 易 确定 ， 初 步 估 计 在 250m 到 350m 深 附 近 。 


图 例 
医 肥 推断 采 空 区 





测 点 /m 
(a) 
图 7.19 工区 成 果 图 
(a) 视 电阻 率 平 面 图 ; (b) 和 典型 剖面 视 电 阻 率 等 值 线 断面 图 
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00 ; 
300 340 380 420 460 500 540 580 620 660 700 740 780 820 860 900 
测 点 号 
(b) 
图 7.19 工区 成 果 图 ( 续 ) 


将 工区 1300 一 1500 号 线 的 300 一 500 点 的 数据 组 成 一 个 三 维 数据 体 ， 经 波 
场 变 换 和 适当 均衡 处 理 后 ， 进 行 合成 孔径 计算 ， 得 到 整个 区 域 的 合成 和 孔径 波 场 . 
三 维 像 图 ， 其 匀 空 效果 图 如 图 7. 20 所 示 〈 见 彩 图 ) ， 可 以 看 出 在 经 过 相关 委 加 
后 ， 突 出 了 有 用 信和 号， 增强 了 信 噪 比 ， 成 像 效 果 明 显 。 





图 7. 20 典型 地 外 成 像 三 维 链 空 图 
将 合成 后 的 数据 进行 偏 移 成 像 ， 得 到 以 深度 为 纵 轴 的 三 维 成 像 图 ， 如 图 
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7.21 所 示 〈 见 彩 图 )。 图 中 蓝 色 代表 波 场 负 值 区 ， 蓝 色 越 次 则 表示 负 值 越 大 ; 
红色 代表 波 场 正 值 区 ; 而 黄色 、 绿 色 为 过 渡 色 ， 表 示 波 动 幅 值 较 小 的 地 方 。 从 
图 中 可 以 看 到 ， 和 电阻 率 的 断面 图 相 比 ， 此 时 的 分 层 已 经 非常 明显 ， 通 过 对 成 
像 结 果 进 行 解释 ， 认 为 有 两 层 采 空 区 分 布 ， 上 层 采 空 区 顶板 位 于 地 下 140m 
处 ; 下 层 采 空 区 顶板 位 于 地 下 280m 处 ， 由 于 充 水 的 影响 ， 导 致 低 阻 异常 放 
大 ， 采 空 区 的 分 布 范围 也 会 有 一 定 的 放大 。 经 验证 解释 结果 与 实际 地 质 情况 吻 
合 较 好 ， 充 分 证 明了 波 场 变换 后 合成 孔径 方法 的 有 效 性 。 





图 7.21 三 维 波 场 成 像 效 果 图 
(a) X PJU]; (b) 立方 向 切片 
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图 7.3 A 型 地 电 模 型 波 场 速度 分 析 结 果 





图 7.4 A 型 地 电 模 型 流 场 延 拓 成 像 结 果 
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(a) | E 
u (b) 
图 7.5 模型 三 维 示 意图 
(a) 三 维 模 型 立体 图 ; (b, 主 齐 面 平面 图 





图 7.6 视 电 阻 率 剖 面 
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图 7.7 模型 计算 结果 
(a) 合成 孔径 成 像 前 波 场 变换 结果 ; 


(b) 合成 孔径 后 波 场 变换 结果 


= XVI > 


9 


10 11 





偏 移 成 像 断 面 图 
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图 7.8 
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计算 结果 


7.10 ”模型 
变换 结果 ; 


(a) 合成 孔径 成 像 前 波 场 


s 


波 场 变换 结果 


(b) 合成 孔径 后 
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图 7.11 





图 7.12 BEER REMER i m 
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图 7. 14 
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图 7.13 虚拟 波 场 速度 分 析 三 维 结 来 
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实测 资料 虚拟 波 场 三 维 切片 结果 
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掌 子 面 zim I zim 
图 7.15 不 同方 位 水 体 不 良 地 质 体 波 场 变换 三 维 结果 
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图 7.16 KRR BERKER RREN 





虚拟 时 间 /st 


图 7.17 波 场 三 维 成 像 效 果 图 
(a) X isi] I (b) 工 方向 切片 
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图 7.18 三 维 延 拓 成 像 效果 
(a) 去 掉 低 能 反射 部 分 ; (b, 去 掉 高 能 反射 部 分 
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测 点 号 
(b) 
图 7.19 工区 成 果 图 
(a) 视 电 阻 率 平面 图 :(b) 典型 剖面 视 电 阻 率 等 值 线 断面 图 
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图 7.21 三 维 波 场 成 像 效 果 图 
(a) 时 方向 切片 (b) 立方 向 切片 
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